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ABSTRACT 



The Naval Postgraduate School possesses software designed to perform thermal 
analysis of electronic components. At the core of this package is a model builder whose 
purpose is to generate a thermal model for use in steady state and transient thermal 
analyzers. The current version of the model builder requires excessive amounts of time 
for data input and model verification. This thesis describes the development of a model 
builder specifically designed to reduce the time required to model a printed circuit board 
containing up to four copper layers. 
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I. INTRODUCTION 



The Naval Postgraduate School possesses software designed to perform thermal 
analysis of electronic components. The software package contains two routines de- 
signed to generate a thermal model (in the form of an ASCII input data file) to be read 
by the thermal analyzer program. The first routine is a generalized model builder used 
in developmental stages, as well as an editor used to modify existing models. The second 
routine was designed to work with models that have a specific geometric configuration. 

In order to accurately model electrical components, the structure must be subdivided 
into small, equal sized subvolumes. The centroid of each subvolume is referred to as a 
node, and due to the assumed isothermal nature of each subvolume, a node is considered 
as representative of the total subvolume. 

Producing an accurate thermal model, requires the design engineer to deal with 
enormous amounts of equations and temperatures, each describing an indiviual node. 
Additionally, each node is connected to adjacent nodes by thermal conductances. As the 
desired accuracy increases, the number of nodes and equations also increases. Modeling 
the electronic component without the aid of a computerized model builder is a task re- 
quiring inordinate amounts of time. 

Current versions of the thermal analyzer contain provisions for the generation of the 
node equations and node interrelationships; however, data input is still a manual, time- 
consuming process. At present there exists a need for a program that will generate a 
data-input file for the thermal analyzer that is both generated through a menu-driven 
routine, and allows the design engineer enough flexibility to model the electronic com- 
ponent to suit his or her needs. 

This thesis describes the development of a model builder designed specifically to re- 
duce the time required to model the copper and epoxy layers of a printed circuit board. 
A typical printed circuit board configuration is shown in Figure 1. The printed circuit 
board may contain up to four copper layers. Additional features of the model builder 
covered are: 

1. The capability of working in SI or English units. 

2. The choice of a total of sixteen aspect ratios. 

3. The provision for up to 740 nodes. 

4. The ability to specify the percent of copper coverage for each layer. 
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5. The ability to specify dimensional parameters for each individual copper layer. 

6. The ability to input heat dissipation using several methods. 

7. The provision for six ambient temperatures. 

8. The automatic calculation of conductance values based on user input. 
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Figure 1. Typical printed circuit board configuration. 
Source: Reference 2. 



II. A JUSTIFICATION FOR THERMAL ANALYSIS 



A printed circuit board is a conglomeration of organic and inorganic materials with 
external and internal wiring, which allows electronic components to be physically sup- 
ported, and electrically connected. [Ref. 1J 

Over the past several decades, printed circuit board technology has developed sub- 
stantially. Early printed circuits were fabricated by printing a pattern of polymer resist 
on a copper plane and then chemically etching. Holes drilled in the laminate held the 
component leads that were soldered to the copper-printed patterns. The technology has 
progressed in developing the sophistication and uses of the printed circuit board inter- 
connections. 

Today, the basic functions of the printed circuit board are the same; the intercon- 
necting copper signal lines join two I/O leads from two different components. The 
components may be resistors, inductors, capacitors, or semiconductor chips. When 
applying multichip technologies, there may be hundreds of components attached to the 
printed circuit board. The ever-increasing level of complexity of printed circuit boards 
has forced primitive boards using a few yards of printed wiring in the 1960s, to evolve 
into sophisticated, multi-layered structures requiring kilometers of printed wiring. This 
increase in the circuit board level of sophistication can be attibuted to the integration 
of semiconductors and an increased need for I/O capabilities. [Ref. 2] 

The most complex printed circuit boards contain kilometers of copper intercon- 
nection, roughly 50 to 100 microns wide, and half as thick. These boards distribute 
KIVIm 2 of power internally in very densely packed layers of copper. The drive for higher 
performance means that there is a greater requirement for power-handling and cooling 
capabilities. [Ref. 2] 

It is the responsibility of the designers to ensure that cooling on the printed circuit 
board is adequate under all possible load conditions in order to allow proper perform- 
ance of the individual components and the board as a whole. Therefore, it is imperative 
that the designers understand and are able to predict the temperature distribution on 
multilayered structures prior to prototype production. The overriding reasons for per- 
forming a precise thermal analysis are to increase component reliability, ensure proper 
material selection, reduce the possibility of catastrophic thermal failure, and guarantee 
electrical performance. 
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A. RELIABILITY 



Reliability is defined as the probability that a component is functioning as designed, 
while failure is defined as the probability that a component is not functioning as de- 
signed. There is a predictable relationship between the operating temperature of elec- 
tronic components and reliability. The materials used in the fabrication of these 
components have thermal limitations, and should these thermal limitations be exceeded, 
the physical and chemical properties of the material are affected, and the device fails. 
[Ref. 1] 

For a large number of components, a typical plot of failure rate as a function of time 
is shown in Figure 2. Failures at short times are called early fails, or infant mortality, 
while failures at long times are called wearout fails because they result from usage. At 
all times, failures can occur from intrinsic mechanisms, or from random overstress. [Ref. 

3] 

Provided the device has been adequately designed, early failures can arise as a result 
of manufacturing defects. Defects that occur early on, or the "burn-in" period, are 
considered to be the result of poor or inadequate quality control mechanisms in the 
manufacturing process. [Ref. 2J 

Of greatest concern are the failures that occur during the useful life of the device 
because the probability of failure during this period should be nearly zero. Should a 
device fail during its useful life, the. probable cause of failure would be due to a variety 
of external factors, and are unpredictable. [Ref. 2] 

As time and usage progress, the terminal period of wearout is encountered, usually 
well past the the end of system life (EOL), in which the failure rate increases. Con- 
ductor electromigration is a typical example of a wearout mechanism, in which the 
electron flow itself causes irreversible mass flow, which causes the formation of voids and 
consequent conductor failure. Printed circuit boards seldom have significant failure rates 
when they are produced with sufficient characterization and control. When failures do 
occur, they can normally be attributed to either manufacturing defects or lack of integ- 
rity to the design specifications. The most notable concerns are the resin and laminate 
effects that can give rise to insulation integrity. As printed circuit boards become more 
densely packaged and continue to spread into more unconventional environments such 
as homes, automobiles or marine engines, there is a need to establish improved materials 
and process controls for even better reliability. Thermal analysis of the printed circuit 
board becomes one of the fundamental facets of the design process in order to better 
characterize the printed circuit board. [Ref. 2] 
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The infant mortality, - or early fail 
portion of the curve is the result 
of failures during the "burn-in" 
period. The high failure rate shown 
on the right side of the curve is due 
to usage, or "wear-out." 




Figure 2. Failure rate as a function of time for typical component. 
Source: Reference 3. 
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B. MATERIAL 



The fabrication of printed circuit boards results in the joining of different materials. 
The materials selected can have a significant impact in thermal properties of the printed 
circuit board. In the largest circuit boards hundreds of amps may be switched at once. 
As packaging densities increase thermal, mechanical, electrical, and chemical coupling 
becomes very strong. In view of the complexity of electrical structures today, computer 
modeling of total thermal response of a printed circuit board is a requirement conducive 
to understanding one of many interactions. [Ref. 2] 

C. THERMAL FATIGUE AND CATASTROPHIC THERMAL FAILURE 

Printed circuit boards are comprised of dissimilar materials that expand at different 
rates of heating. Table 1 shows the thermal coefficients of expansion of materials com- 
monly used in the fabrication of printed circuit boards. The differential expansion of 
mismatch must be accommodated by the various elements on the board. The increasing 
levels of packaging densities and board complexity dictate the need for designing a 
thermal environment that can accommodate the diverse components that are in close 
proximity to each other. [Ref. 4] 



Table I. MATERIALS USED IN PCB PRODUCTION 



Chemical Name 


Coefficient of Thermal Expansion 1 0 -7 / °C 


Polyethylene 


650 


Bismaleimide 


500 


Polyamide imide 


360 


Polyarylate 


625 


Peek 


400 


Polvether amide 


500 


Polyimide 


500 


Polytetrafluorethylene 


700-1200 


Epoxy glass cloth 


170 (x,y); 600 (z) 


PBZT 


90 (x,y); 200 (z) 



Source: Reference 2. 
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As previously mentioned it is very important that the board designers have an 
understanding of the operating environment in which the board will be operating in or- 
der to incorporate into the design the tolerances that will allow the product to operate 
reliably. 

Catastrophic thermal failure is defined as an immediate, thermally induced, total loss 
of electronic function in a specified component. This type of failure comes as a result 
of excessive temperature, or a thermal fracture. Catastrophic failure comes about as a 
result of many factors including the operating environment, equipment history, me- 
chanical loading, and operational modes of the component. Although it is difficult to 
predict the temperature at which thermal failure may occur, it is possible to establish 
with the aid of thermal analysis the boundaries at which the board can be expected to 
operate reliably and within its useful operating life cycle. [Ref. 1,2] 
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III. HEAT TRANSFER 



The degradation of the heat flow capabilities of a printed circuit board can lead to 
reliability problems due to excessive operating temperatures. It is imperative that de- 
signers incorporate into the board the capability to maintain temperatures within upper 
operational limits while operating in all possible environments in which the board will 
be exposed. [Ref. 2] 

Heat transfer is defined as all energy flows that arise as a result of temperature dif- 
ferences. Because the components mounted on printed circuit boards and, indeed, the 
printed circuit boards themselves are not one hundred percent efficient, heat is gener- 
ated. The primary modes of heat transfer are conduction and convection. Conductive 
modes include mechanical thermal contact and solid thermal interfaces between materi- 
als, such as copper, solder, or epoxy layers. Convective modes include natural and air- 
forced cooling, and forced liquid cooling. Radiation is also a factor; however, it is not 
as significant as conduction and convection at the temperatures in which printed circuit 
boards operate. [Ref. 5] 

A. CONDUCTION 

Conduction is the transfer by molecular motion of heat between one part of a body 
to another part of the same body or one body and another in physical contact. [Ref. 1] 

For the case of conduction, the heat flow equation is the basis for understanding this 
behavior 



Q = -kA 



AT 

AX 



( 1 ) 



where 

Q = heat flow along the thermal gradient, 

A = the area through which the heat is flowing, m 2 

W 

k = the thermal conductivity of the material, — _ 
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A T 

—~-rr = change in temperature per unit length, 
IXa. 



°_C 

m 



Rearranging Equation (1) leads to 

R-th = ^q~ ( 2 ) 

which is the thermal resistance 



R-th ~ 



AX 

kA 



in 



IQ. 

w 



( 3 ) 



Equation (3) demonstrates that thermal resistance will increase with an increase in 
path length for heat flow, AX, with a decrease in area of heat flow, A, or change in 
conductivity k to a low'er value. 

1. General Equation of Heat Conduction 

The general equation of heat conduction is 



J- (k lL ) + JL 

dx K dx ’ oy 



(k J-L ) + JL( k JLL 

{k dy }+ dz [k cz 



s . r cT 

) + pe- 



rn 



w'here 



p = density, — — 
m 

C = specific heat, 

T = temperature, °C 

x,y, and z = cartesian coordinates, m 

t = time, sec 



k = thermal conductivity, 



W 

m — °C 



q = internal heat generation, — — 

m 
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Assuming k , C and p are independent of temperature, direction, and time, the resulting 
equation is 



£t_ £t_ £t_ ± _]_ dT_ 

dx 2 dy 2 dz 2 k « 



where 



a = thermal diflusivity, 



, 2 
k m 

pC ’ sec 



(5) 



Several variations of the general conduction equation exist. Fourier's equation, 
which contains no heat sources is 



o 2 r e 2 r a 2 r _ 1 bt 

ox 2 dy 2 dz 2 a ot 



(6) 



Another variation, known as Poisson's equation, solves a system in which tem- 
perature is not time dependent 



c 2 T 

~ 2 
OX 



oy 





(7) 



The last variant of the general equation of conduction is intended for a system 
operating in a steady-state condition, and does not contain any heat sources. This 
equation, known as La Place's equation is 



c 2 r d2j_ £t 

- 2 + ,2 + 2 

ox dy oz 



= V 2 T= 0 



( 8 ) 



2. Single Plane Slab Modeling 

For illustrative purposes, a single plane slab, with 7, and T 2 being the face 
temperatures, is considered. If Equation (8) is limited to only one coordinate then 

-7T = 0 (9) 

dx 

If Equation (9) is integrated twice, and boundaries are established, then the temperature 
distribution across the slab can be expressed as follows 
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( 10 ) 



T=T x -j-(T,-T 2 ) 

If Equation (10) is substituted into Equation (1), then a solution for the heat flow across 
the single plane slab can be obtained 

q = ~kA\_ 1= J j-(T,-T 2 ) (11) 

It is known that Ohm's law relates the voltage drop across a resistor with the 
current flow through the resistor, V=IR. Ohm's law is analogous to Equation (1), 
where the current flowing through a resistor is equivalent to heat flow’, and the electrical 
resistance is equivalent to the thermal resistance expressed in Equation (2). 

Ohm's law r can be expressed as 

'-T < l2 > 

w'here the analogy between Ohm's law’ and Equation (1) is 
Current 7<=>Heat flow’ q 

Potential E<^Temperature difference A T 
Resistance /?<=> Thermal Resistance R 

It can be observed that for the heat flow across a simple plane slab described 
by Equation (11), the thermal resistance is described by Equations (2) and (3). 

B. CONVECTION 

Convection is defined as the process by which thermal energy is transferred to or 
from a solid by a fluid flowing past it. If the movement of fluid is a result of a temper- 
ature differential, then the process is called free or natural convection. When natural 
convection is present, the movement of fluid can be accelerated by increasing the tem- 
perature differential. When a pressure differential is introduced to force the movement 
of the fluid by using a pump or a fan, the process is called forced convection. [Ref. 
1,5] 

Newton's law’ of cooling states that the heat flow' through a body is proportional to 
the normal area and the temperature difference between the body and the surrounding 
fluid [Ref. 2). To make it an equality, a proportionality constant, h , is introduced. The 
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proportionality constant is the surface heat transfer coefficient. Newton's law of cool- 
ing can be expressed as 

q = hA{T 0 -T f ) (13) 



where 

h is a proportionality factor known as the surface heat transfer coefficient 



Comparing Newton's law of cooling with Fourier's law leads to an expression re- 
lating the surface heat transfer coefficient to thermal conductivity, the surface fluid 
temperature differences, and the wall temperature gradient of the fluid 



h = 



q 

AAT 



-(1L) 
* dy [ 

AT 



(14) 



Consequently, correlating heat transfer coefficients must be based on the depend- 
ence of h on the thermal conductivity of the fluid and on the ratio of the wall temper- 
ature gradient to the temperature difference. [Ref. 1] 

1. Electrothermal Analog 

For the case of convective heat transfer the thermal resistance is represented by 



R = 



hA 



(15) 



Consequently, the total thermal resistance is now defined for a single slab and convective 
heat transfer on its two faces as follows 



R — ~7~r + ~rj + ~i~~7 — ~~r C "T h ~r~ + ~r~ 1 

h x A kA h 2 A A h ] k h 2 



The heat transfer equation can now be represented as 
(7) - T 2 ) 



q = 



AT_ 

R 






(16) 



(17) 



C. RADIATION 

The third mechanism for heat transfer is radiation in the form of electromagnetic 
waves. The rate at which a body radiates thermal energy is proportional to the area 
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of the body and to the fourth power of the absolute temperature [Ref. 5] 
found empirically by Josef Stefan in 1879, is written as 

q = eoAT 4 

where 

q = power radiated, W 

2 

A = area, m 

e = emissivity of the body, a value between 0 and 1 

o [{/ 

a = Stefan-Boltzmann constant, 5.6703 x 10 — r 

m 2 .h A 

When radiation falls on an opaque body, part of the radiation is reflected and part 
is absorbed. Light-colored bodies reflect most of the radiation, whereas dark bodies 
absorb most of it. [Ref. 5] 

Materials employed in the manufacture of electronic components are classified as 
gray. Gray bodies are diffusely-reflecting opaque surfaces. These surfaces reflect equal 
amounts of energy over the thermal radiation spectrum in all directions. [Ref. 2] 

1. Transformation of the General Radiation Equation 

The use of the thermal radiation equation in analytical studies is a difficult task 
due to the fourth power relationship with temperatures. The complexity of the calcu- 
lations involved with radiation dictate the need for the aid of a computer [Ref. 2]. The 
general equation for radiation interchange is 

q = oF a F e A{T A s - t T ) ( 19 ) 



. This result, 
(18) 



where, as before, o is the Stefan-Boltzmann constant and 

F a = shape factor accounting for source and receiver arrangement. 

2 

A = area, m 

F e = emissivity factor accounting for properties of the source and receiver. 
T s = temperature of the source, °A" 
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T r = temperature of the receiver, °K 



It must be observed that for the radiation case absolute temperature is the norm. 

Radiation based heat transfer is represented by transforming the general radi- 
ation equation, Equation (19), into a form compatible with Fourier’s law'. Linearization 
of the general radiation equation is the method used to produce the desired result. The 
difference between the tw r o fourth powers can be reduced to: 



(7? - Tf) - (7? + 7?X7? - 7?) 

= (7? + Tj)(T s + T r )(T s - T r ) (20) 

Inserting this into Equation (19) results in 

? = <r W(7? + 7?)( 7) + T r )( T s -T,) (21) 

A radiative heat transfer coefficient may be defined as 

K - oF a F e (Ts + T?XT, + T r ) (22 a) 

or 

K - oF.Fjf, +T]T r + TJ? + 7?) (22*) 



Substitution of h r into Equation (21) shows that radiative heat transfer can now 7 be 
treated similarly to convection at the boundary [Ref. 6). Thermal resistance in the case 
of radiation heat transfer can now be denoted as 



Now that all three methods of heat flow have been discussed, a graph depicting 
the electrothermal equivalent is shown in Figure 3. Here the plane slab has both radi- 
ation and convection on its tw r o faces, and the nonlinearity of the radiation coefficient, 
h r , leads to a more detailed analysis procedure. 

When considering the calculation of heat transfer by radiation, it is usually 
necessary to approximate real material behavior with gray-body idealization. There is 
also a great deal of data required in order to model radiation for a real body not only 
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because the surfaces are non-black, but also because the non-uniformity of the body's 
surfaces creates greater complexities in the collection and processing of data. [Ref. 7] 

Basic problems, such as the one shown in Figure 3, defy a simple solution be- 
cause of h r . Therefore, a computer-aided approach becomes necessary in solving even 
elementary problems. Models exist for aiding the user in solving radiative heat transfer 
problems; however, their development is not within the boundaries of this thesis. 
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IV. FINITE DIFFERENCE ANALYSIS 



Finite difTerence methods represent the most appropriate approach in determining 
the temperature distribution within an electronic component. By using a finite differ- 
ence approach, complex problems involving difficult geometries, non linearities, or 
complicated boundary' conditions can be successfully approached. The purpose of this 
section is to explain the fundamental concepts behind finite difference analysis and its 
use in solving differential equations. [Ref. 6] 

A. FUNDAMENTALS OF FINITE DIFFERENCE ANALYSIS 

The general equation for heat transfer by conduction, Equation (4), can be ex- 
pressed in a three-dimensional variant as 



_d_ 

dx 



(k 



dT 

dx 




dT 

dy 




cT 

dz 



) + q t = pc 



dT 

dt 



(24) 



where all the equation components have been previously defined. In order to formulate 
a problem using Equation (24), it is first necessary to approximate the first and second 
derivatives. 

1. First and Second Derivative Approximation 

The derivative of a function at a point can be formulated by a finite difference 
approximation. It is first necessary to determine the slope of a line tangent to the point 
of interest on the plot of temperature as a function of x, y, or z. (Figure 4). In choosing 
x as the coordinate system, it is observed that the slope of the line (m) is defined as 



lim 

Ajc-»0 



J\xq + Ax) -TI-Vp) 
Ax 



(25) 



The term of Equation (25), [j^x 0 + Ax) — y(x 0 )]/Ax is referred to as the difTerence quo- 
tient, and is the ratio of the change in the value of the function at x 0 and x 0 + At to the 
change in x. The limit of the difference quotient is called the derivative of the function 
at x 0 . [Ref. 8] 



/'(x 0 ) = lim 
Ajc->0 



Axq + Ax) -Axp) 
Ajc 



( 26 ) 
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T 




Figure 4. First and second derivative approximation. 
Source: Reference 8. 
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By letting y[x) = 7~(x), the derivative for this particular curve at the point x 0 can be de- 
rived 



/'(* o) = lim n 

Ajc-»0 



T(x 0 + Ax) - 7Yx 0 ) 

Ax 



As Ax -+• 0 Equation (27) becomes 
, dT T(x o + Ax) - T(x 0 ) 

= S 

Thus, for Ax finite, but very’ small 






(27) 



(28) 



(29) 



The second derivative, using the definition of a derivative, will be the ratio of the 
change in value of the first derivative to the change in Ax [Ref. 8]. Again, looking at 
Figure 4, the second derivative is 



/''(•%) = lim 



/'(*o + A*) -/'(xp ) 
Ax 



(30) 



f"(x o) 



d 7 T _ T(x, + 2 Ax) - 277x 0 + Ax) + 1\x Q ) 
dx 2 Ax 2 



( 31 ) 



For both Equations (28) and (31), the approximations are valid as long as Ax is small. 

Taylor's theorem states that a function can be approximated by a polynomial 
of the form 

Ax) =A a ) +/'(«)(* - a) + ~ 2^ ( x ~ a ) 2 + - + ~ ^ ( x ~ a T ( 32 ) 



where the polynomial is for the function expanded about x = a [Ref. 8]. The Taylor se- 
ries for the function T(x) at T(x 0 + 2Ax) is 



r(x 0 -I- 2Ax) = T(x 0 + Ax) + 



dT(x 0 + Ax) 

dx 



Ax + 



1 d 2 T(x 0 + Ax) 

2 dx 2 



Ax 2 + 



(33) 



1 d 3 T(x 0 + Ax) 
6 dx 3 



Ax 3 + ... + 



j d n T(x 0 + Ax) 
n\ dx n 



Ax n 
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and for the function T(x ) at T(x 0 ) : 



7Y \ T( . a , dT{x Q + Ax) x d T(x 0 + Ax) 2 

/ (x 0 ) = T(x 0 + Ax) 7 Ax + — 7 Ax 



dx 



dx 



1 d 2 T(x 0 4- Ax) 3 i d n T(x 0 + Ax) 

—Ax +... + — Ax 



dx 



n\ dx n 



(34) 



If Ax is kept small, then the terms above second order become negligible. When 
Equations (33) and (34) are added, an expression for the second derivative can be ob- 
tained 

d 2 T{x 0 + Ax) 2 

T(x 0 + 2 Ax) + 7(x 0 ) = 2 7(x 0 + Ax) + — 2 Ax 2 (35) 

dx 



d 2 T(x 0 + Ax) 

dx 1 



T(x 0 4- 2Ax) — 2 7(x 0 + Ax) + T(x 0 ) 
Ax 2 



=/"W 



(36) 



As previously stated, the configuration to be analyzed can be subdivided into 
small, finite subvolumes considered to be isothermal. The centroid of each subvolume 
is called a node, and the node is representative of the entire subvolume. Nodes may are 
connected to their adjacent nodes through thermal resistances and the nodal analysis for 
the node temperatures can best be solved by means of a computer aided model [Ref. 8]. 

The first law of thermodynamics states that energy can neither be created nor 
destroyed, but can be transformed from one form to another. An energy balance may, 
therefore, be formed on the typical node shown in Figure 5. For the sake of simplicity, 
the environment is composed by a single node, node 5 (Figure 6). Close observation of 
node 5 shows that it is connected to nodes 2, 4, 6, 8, 14, and 100. If any energy is di- 
rectly applied to or removed from this node, it would become part of the equation. The 
node equation for node 5, with rate of heat input q„ becomes 

#2 + q$ + + q% + <7i4 + <7ioo — qi = 0 (37) 

where each one of the q values with a numerical subscript represents the rate of heat 
flow from node 5 to the node indicated by the subscript. [Ref. 1,6,8] 
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Figure 5. Graphical representation of a single node. 
Source: Reference 6. 
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Figure 6. Node arrangement on typical printed circuit board. 
Source: Reference 8. 
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Equation (37) does not display the node temperatures; however, q = AT/R and 
Equation (37) can be combined to produce 



?2 — 




kAxAz 

Ay 



(7s - T 2 ) 



?4 = 



M 

L 



AT = 



kAyAz 

Ax 



(T 5 ~T a ) 



<76 = 



M 

L 



AT = 



kAyAz 

Ax 



( T 5 -T 6 ) 



<7s = 




kAxAz 

A y 



(T 5 ~ T & ) 



(38a) 

(386) 

(38c) 

(384) 



<7l4 — 



M 

L 



AT = 



kAxAy 



7 *. 4 ) 



(38c) 



<7100 — 



M 

L 



AT = 



/:A.xAy’ 

Ac 



(7*5 - T'iqo) 



2 



(3S/) 



where r i00 is on the top face of node 5 and, Ax, Ay, A z = the distance between nodes in 
the x,y, and z directions respectively. Assuming Ajc = A^ = Az = 1, then the node 
equation becomes 

k(T 5 - T 2 ) + k(T 5 - r 4 ) + k(T s - T 6 ) + (39) 



«T S - r 8 ) + k(T 5 - r 14 ) + 2 k(T 5 - r 100 ) = q t 



or 



-r 2 -r 4 + 7r 5 -r 6 -r 8 -r 14 = -|- + 2r 100 (40) 

where the terms on the right hand side of Equations (39) and (40) are knowm values. 
[Ref. 1,6,8] 

A node equation is written for each node in the configuration, and in Figure 
8, which contains 18 nodes, 18 node equations would be required to obtain the 18 un- 
known node temperatures. A larger problem having several hundred nodes would re- 
quire prohibitive amounts of time to solve or even enter data for computer-aided thermal 
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analysis. That is why a computer-aided model builder capable of constructing the 
needed input data for thermal analysis, is so desirable. The thrust of this thesis is not 
to solve the equations, but to allow the user to model a printed circuit board to his or 
her specifications with minimal effort, and then produce a data file compatible for use 
in the thermal analyzer. 
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V. THE MODEL BUILDER 



As previously indicated, the model builder currently in use in the thermal analysis 
software package requires time-consuming manual data entry. This thesis proposes that 
a model builder, PCB, can be incorporated into the software package. PCB is a menu- 
driven, user-friendly program which will assist the user in developing a thermal model 
of a printed circuit board with pre-defined geometries. PCB will generate a properly 
formatted ASCII output data file for use with the thermal analyzer. The current version 
of the thermal analyzer is designed to handle a maximum of 300 nodes [Ref. 9], An 
updated version capable of handling the amount of nodes generated by the model builder 
is under development; however, abridged versions of the model builder have been suc- 
cessfully tested with the current version of the model builder. This chapter addresses 
many aspects of the model builder, including a terminal session demonstrating the design 
process used in modeling a circuit board, and the interface between the output data file 
and the thermal analyzer. Also addressed is the methodology used in validating the PCB 
model builder output data file. Finally, limitations and possible improvements to the 
model builder are discussed. 

A. FEATURES 

PCB is a menu driven interactive program providing the user with a variety of sys- 
tematic choices progressing in a logical sequence that eventually will lead to the devel- 
opment of the model. Upon entering the program, the user will be provided with the 
option of viewing an optional overview. First-time users should consider reading the 
overview. The program is designed to operate with upper-case letters, therefore, it is 
suggested that the caps-lock key is pressed before starting the program. 

Following the overview option the user will be asked to choose either SI or English 
units. For this choice as well all others, the program provides the user with the op- 
portunity to correct any erroneous entries. 

The third section queries the user for the number of copper layers desired in the 
design. The allowable number ranges from one to four layers. For every choice of 
number of copper layers, there are four possible board configurations the user can select. 
If one copper layer is selected, the program will display the following choices of aspect 
ratios: 

1. 10 BY 36 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 
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2. 10 BY 35 PROVIDES A 1:3.5 RATIO WITH 700 NODES. 

3. 12 BY 25 PROVIDES A 1. 2:2.5 RATIO WITH 600 NODES. 

4. 15 BY 20 PROVIDES A 1.5.2.0 RATIO WITH 600 NODES. 

When two copper layers are selected, the user will be provided with the following 
choices: 

1. 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES. 

2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES. 

3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES. 

4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES. 

If three copper layers are selected, the user will have the following choices: 

1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 

2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 

3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 

4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 

By selecting four copper layers, the program 'will provide the user the following 
options: 

1. 5 BY 18 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 

2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES. 

3. 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES. 

4. 13 BY 6 PROVIDES A 1:1 RATIO WITH 624 NODES. 

After selecting the desired aspect ratio, the user must specify the epoxy character- 
istics. The program will prompt the user to enter the epoxy layers length, width, 
thickness and thermal conductivity. An option available to the user will be the choice 
of a default epoxy thickness value, 0.0625 in., (0.15875 cm.). 

Once the epoxy characteristics have been entered, the program will then display the 
corresponding length and thickness for the copper layers selected. The length and width 
for the copper layers will be the same as those of the epoxy layers. The user will then 
be provided with the option of entering the thickness of each selected copper layer either 
by length or by weight. If input by length is selected, the units should be consistent with 
the system of measurement selected; however, if entry of thickness by weight is selected, 
then the entries should be in ounces. A copper layer weighing one oz. is equivalent to 
a copper layer 0.0014 in. thick with a surface area of 144 in 2 . The program adjusts 
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thelayer thickness to the surface area provided automatically. Entry of thickness by 
weight is a standard practice used in industry. [Ref. 2] 

After entering the copper layer thickness the program will query the user for the 
thermal conductivity of copper. After this entry, the program will redisplay all infor- 
mation pertaining to the copper layers and will require the user to confirm all entries in 
order to proceed. 

The following section requests the input of initial and ambient temperatures. PCB 
has six ambient temperatures. 

After the temperatures have been entered the program requires the user to input the 
external heat sources. PCB provides for heat input into the upper copper layer. There 
are four methods of external heat input from which to choose. The first alternative al- 
lows for a total rate of heat applied to the upper surface. An entry for this choice would 
be divided by the number of nodes and distributed uniformly. The second choice pro- 
vides for the entry of average heat per unit area. The third alternative gives the user the 
ability to enter heat in specifically designated nodes. The last option provides for no 
heat input and was developed to test compatibility with the thermal analyzer. 

The following selection is the percent copper coverage for each copper layer. The 
program will prompt the user to enter the percentage of copper coverage for each layer 
in the printed circuit board. 

Finally, the program requires a name and a title for the file that will be created. This 
file will be in the proper format for use with the thermal analyzer. 

B. THE THERMAL ANALYZER INPUT DATA FILE 

The model builder generates an ASCII data file from the physical characteristics of 
the printed circuit board provided by the user. In order for the input data file to be ac- 
ceptable to the thermal analyzer, it must be in a specific format which is compatible to 
the thermal analyzer [Ref. 9]. Because each value and position of the output data file 
has a meaning to the thermal analyzer, and is not readily identifiable to the user, it is 
beneficial to describe each line and data set of the output data file and their relationship 
to PCB. Figure 7 shows a partial output data file. 

Line one is the title line. It may be left blank or may contain up to 79 alphanumeric 
characters. The user-selected title appears at the top of the data file. [Ref. 6] 

Line two is the problem data line. It has nine entries of which two are under user 
control, the number of nodes under consideration and the unit type. One entry, the 
number of constant temperatures, is preset at six for this specific model. The remaining 
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Figure 7. Sample PCB partial output data file. 
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entries have applications to models associated with heaters, unique exponents, secondary 
heat input, temperature coefficients and curves, and nodes controlling fast heat. These 
entries are not applicable to this model and are preset to zero. [Ref. 6] 

Line three places a zero at three points and is beyond the user's control. Therefore, 
no further discussion is required. [Ref. 6] 

Line four is the problem capability line. This line defines the maximum values for 
the entries in line two. The first entry is 750, which is the number of nodes for which 
the analysis is dimensioned. The number 750 is significant because the first constant 
temperature will be assigned to node number 751. The second entry is 50 which re- 
presents the largest possible number of constant temperatures in accordance with the 
analyzer dimension statement. The third entry is set to 6 and does not change. This 
entry is related to heaters and is not applicable to the model. The balance of the entries 
in line four represent a listing of data sets that are required for the particular analysis 
at hand. PCB uses three data sets that will be discussed in what follows. [Ref. 6] 

The fifth line contains five values that relate to the accuracy level that the thermal 
analyzer will achieve. These entries are preset. The first value provides the level of ac- 
curacy between iterations. The accuracy level number is critical because too small a 
tolerance will cause the computer to run for excessive amounts of time, and too large a 
number will provide inaccurate results. The second value is the damping factor used 
between iterations in order to prevent temperature oscillations between iterations. The 
third number provides the maximum number of iterations. If erroneous data is entered, 
the computer will not run for excessive amounts of time. The fourth value is the con- 
vergence factor which adjusts the damping in order to close the critical value. The fifth 
entry is the initial temperature at which the iterative process begins. This value is sup- 
plied by the user. [Ref. 6] 

Line six contains the temperature dependent coefficients, and is not used in this 
model. Line seven contains up to 50 constant temperature inputs. This model has six 
ambient temperatures entered by the user. [Ref. 6] 

The following lines contain all pertinent information concerning the n-node 
equations. Each node requires two lines of data. Even-numbered lines are used for 
specifying the nodes that interact with the node in question and the modes by which this 
interaction takes place. For example, a line of the form: 

6 7551 21 7521 111 7511 2411 
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is for node number one of aspect ratio selection one. The first entry indicates the num- 
ber of connections to that node. The second entry is indicating that node 755, an am- 
bient temperature, is connected to the node in question and the 1 indicates that the 
connection is conductive. The same procedure applies to the rest of the values on the 
line. If the entry is 9991 then the heat input is external. [Ref. 6] 

Other data sets that appear relate to unique exponents, secondary heat, and 
temperature-dependent heat input curves. These data sets are not used in this model. 
[Ref. 6] 

C. PCB MODEL BUILDER SAMPLE PROBLEM 

The following terminal session is a typical example of the printed circuit board 
model builder, PCB, and its capabilities. Figure 8 depicts the printed circuit board 
modeled in the terminal session. 

1. Printed Circuit Board Specifications 

1. Unit system: British. 

2. Copper layers in printed circuit board: three. 

3. Aspect ratio desired: 12 by 10, with 720 nodes. 

4. Epoxy layer length: 7.2 in. 

5. Epoxy layer width: 4.2 in. 

6. Epoxy layer thickness to default value: Yes (.0625 in.) 

7. Epoxy layer thermal conductivity: 0.087 Btu/hr/°F. 

8. Specify copper layer thickness by weight or length: length. 

9. Thickness for the three copper layers, respectively: .20, .10, .11 in. 

10. Copper thermal conductivity: 243.000 Btu/hr/ 0 / 7 . 

11. Initial board temperature: 87 °F. 

12. Upper surface ambient temperature: 86 °F. 

13. Lower surface ambient temperature: 85 °F. 

14. Right surface ambient temperature: 84 °F. 

15. Left surface ambient temperature: 86 °F. 

16. Front surface ambient temperature: 86 °F. 

17. Rear surface ambient temperature: 86 °F. 

18. Heat input to the upper copper layer: 300.23 Btu/hr/°F. 

19. Percent copper coverage for copper layers, respectively: 87.3, 76.2, 65.7 %. 

20. Name of output data file: LAYERS. 
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21. Header of the output data file: THIS IS THE OUTPUT OF THE THERMAL 
ANALYZER. 

2. Terminal Session 

Using the board specifications and Figure 8, the terminal session for the PCB 
program is as follows: 



THIS PROGRAM WAS WRITTEN TO INTEGRATE WITH EXISTING 
THERMAL ANALYSIS SOFTWARE AND TO REDUCE THE AMOUNT 
OF TIME REQUIRED FOR DATA ENTRY. 

WOULD YOU LIKE AND OVERVIEW OF THE PROGRAM PRIOR TO 
BEGINNING? ENTER Y FOR YES AND N FOR NO: N 



************************************************* 

************************************************* 



** ** 

** PLEASE SELECT UPPER CASE LETTERS ** 

** PRIOR TO BEGINNING ** 

** ** 



************************************************* 

************************************************* 



PRIOR TO ENTERING DATA INTO THIS PROGRAM ENSURE THAT 
YOU HAVE A DRAWING OF YOUR DESIGN AND ALL PERTINENT 
DATA. 



PRESS <ENTER> TO CONTINUE 
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Figure 8. Printed circuit board modeled in the terminal session. 
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THIS PROGRAM IS CAPABLE OF OPERATIONS IN EITHER SI OR 
ENGLISH UNITS. AFTER THE SELECTION OF THE UNITS, ALL 
ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION. 

S FOR SI NOTATION 
E FOR ENGLISH NOTATION: E 



YOU HAVE SELECTED ENGLISH NOTATION. 



IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 



THE PCB CAN HAVE UP TO FOUR COPPER LAYERS 

HOW MANY COPPER LAYERS DO YOU DESIRE? 

PICK A NUMBER 1 TO 4 : 3 



YOU SELECTED 3 COPPER LAYER (S) FOR THE PCB 

IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 
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YOU SELECTED 3 COPPER LAYERS, GIVING YOU THE FOLLOWING ALTERNATIVES 



1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 

2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 

3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 

4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 

PLEASE SELECT A NUMBER 1 THROUGH 4: 3 



YOU SELECTED NUMBER 3 OF THE FOLLOWING 4 ALTERNATIVES: 



1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 

2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. 

3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. 

4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. 



IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 
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******************* ************** ************************* ************** 

************************* * *EPOXY CHARACTERISTICS ************************ 
************************************************************************ 

ALL ENTRIES ARE IN ENGLISH NOTATION 

ENTER EPOXY LAYER LENGTH (in): 7.2 
ENTER EPOXY LAYER WIDTH (in): 4.2 

THE EPOXY LAYER THICKNESS WILL DEFAULT TO 0.0625 in (0.15875 cm) 

DO YOU WANT TO CHANGE THE EPOXY THICKNESS? (Y OR N) : N 

ENTER EPOXY LAYER THERMAL CONDUCTIVITY (Btu/hr/F ): 0.087 



YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE EPOXY LAYER: 



1.) 


LENGTH : 


7.2000 


in 


2.) 


WIDTH: 


4.2000 


in 


3.) 


THICKNESS: 


.0625 


in 


4.) 


k: 


.0870 


Btu/hr/F 



DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES AND N FOR NO: N 
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*********************************************************************** 

**************** COPPER LAYER CHARACTERISTICS *************************** 
*********************************************************************** 

UjL ENTRIES ARE IN ENGLISH NOTATION 



THE COPPER LAYER LENGTH IS THE SAME AS THE EPOXY LAYER: 7.2000 in 

THE COPPER LAYER WIDTH IS THE SAME AS THE EPOXY LAYER: 4.2000 in 

(OU SELECTED 3 COPPER LAYER (S) FOR THE PCB. 

fOU NOW WILL BE ASKED TO ENTER THE COPPER LAYER THICKNESS 

)0 YOU WANT TO SPECIFY THICKNESS BY LENGTH OR WEIGHT? 

•INTER l if you want to enter length, w if you want to enter WEIGHT: l 



•INTER THE THICKNESS FOR LAYER 1 (in): .2 

•INTER THE THICKNESS FOR LAYER 2 (in): .1 

•INTER THE THICKNESS FOR LAYER 3 (in): .11 

NTER COPPER LAYER THERMAL CONDUCTIVITY (Btu/hr/F ) : 243 



YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPER LAYER (S) . 



1. ) LENGTH: 

2. ) WIDTH: 

3. ) THICKNESS LAYER 1: 

THICKNESS LAYER 2: 
THICKNESS LAYER 3: 

4. ) k: 



7.2000 

4.2000 
.2000 
.1000 
.1100 

243.0000 



in 

in 

in 

in 

in 

Btu/hr/F 
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************************************************************************ 
***************** * AMBIENT TEMPERATURE INPUTS**************************** 
************************************************************************ 



ENTER THE 
ENTER THE 
ENTER THE 
ENTER THE 
ENTER THE 
ENTER THE 
ENTER THE 



INITIAL BOARD TEMPERATURE (F) : 87 
UPPER SURFACE AMBIENT TEMPERATURE (F) : 86 
LOWER SURFACE AMBIENT TEMPERATURE (F) : 85 
RIGHT SURFACE AMBIENT TEMPERATURE (F) : 84 
LEFT SURFACE AMBIENT TEMPERATURE (F) : 86 
FRONT SURFACE AMBIENT TEMPERATURE (F) : 86 
REAR SURFACE AMBIENT TEMPERATURE (F) : 86 



YOU HAVE MADE THE FOLLOWING AMBIENT TEMPERATURE ENTRIES: 



1. INITIAL BOARD TEMPERATURE: 

2. UPPER AMBIENT TEMPERATURE: 

3. LOWER AMBIENT TEMPERATURE: 

4. RIGHT AMBIENT TEMPERATURE: 

5. LEFT AMBIENT TEMPERATURE: 

6. FRONT AMBIENT TEMPERATURE: 

7. REAR AMBIENT TEMPERATURE: 



87.000 F 

86.000 F 

85.000 F 

84.000 F 

86.000 F 
86.000 F 
86.000 F 



DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES AND N FOR NO: N 



DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES AND N FOR NO: N 
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************* ****************heat inputs******************************* 

********************************************************** 4 ************ 



HEAT INPUT TO THE PCB OCCURS ONLY ON THE UPPER COPPER 
LAYER. HEAT INPUT IS ACCOMPLISHED BY ONE OF THE 
FOLLOWING METHODS : 

1. ENTER AS TOTAL HEAT APPLIED TO THE PCB 

2. ENTER AS AVERAGE HEAT PER UNIT AREA 

3. ENTER HEAT NODE BY NODE 

4. NO HEAT INPUT 

PLEASE SELECT NUMBER 1 THROUGH 4: 1 



YOU HAVE SELECTED NUMER 1 OF FOUR ALTERNATIVES. 
IS THIS THE DESIRED SELECTION? (Y OR N) : Y 



YOU HAVE SELECTED TO INPUT HEAT AS A TOTAL HEAT 
APPLIED TO THE SURFACE. 

ENTER TOTAL HEAT APPLIED TO THE SURFACE (Btu/hr) : 300.23 
IS THIS THE CORRECT ENTRY? (Y OR H): Y 

TOTAL HEAT PER NODE IS: 2.5019 Btu/hr 

<PRESS ENTER TO CONTINUE> 



ENTER 


THE 


PERCENT 


COPPER 


COVERAGE 


(i.e. 


67) 


FOR 


LAYER 


1: 


87.3 


ENTER 


THE 


PERCENT 


COPPER 


COVERAGE 


(i.e. 


67) 


FOR 


LAYER 


2: 


76.2 


ENTER 


THE 


PERCENT 


COPPER 


COVERAGE 


(i.e. 


67) 


FOR 


LAYER 


3: 


65.7 
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YOU HAVE SELECTED: 87.30% COVERAGE, LAYER 1 

76.20% COVERAGE, LAYER 2 

65.70% COVERAGE, LAYER 3. 

IS THIS YOUR DESIRED ENTRY? (Y OR N) : y 



YOU SELECTED LAYERS FOR YOUR DATA FILE NAME 

IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 



THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY INTO THE 
EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES 
NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE 
THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS 
PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND 
SHOULD NOT HAVE ANY SPACES. 



PLEASE ENTER THE DESIRED DATA FILE NAME: LAYERS 



YOU SELECTED LAYERS FOR YOUR DATA FILE NAME 



IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FOR NO: Y 
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ENTER THE DESIRED TITLE TO BE PLACED ON LINE 
NUMBER ONE OF THE OUTPUT DATA FILE: 



THIS IS THE OUTPUT OF THE THERMAL ANALYZER 



DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPUT DATA FILE? 
ENTER Y FOR YES AND N FOR NO: N 



THE OUTPUT DATA HAS BEEN PLACED IN A FILE NAMED LAYERS 



<PRESS ENTER TO CONTINUE> 
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3. The output data file 

Figure 9 shows a partial listing of the model builder output computed in the 
terminal session. A detailed description of the output data file is provided in section B 
of this chapter. 

4. Thermal analyzer output 

The file obtained in the model builder produces an output for 720 nodes. In 
order to use this data file with the thermal analyzer the data file must first be shortened. 
In this particular case, the output file, LAYERS, shortened to eight nodes was entered 
to the thermal analysis program, TASS. Three files are produced by the thermal analysis 
program, TASS. The first two output files consist of detailed and summarized thermal 
analysis results, respectively. The third file is an error message summary' in the event 
there is an error in the input data file. Figure 10 shows the detailed output file of the 
eight node thermal analysis. 

D. PCB MODEL BUILDER VALIDATION 

The output of the model builder was validated by following a three-step process. 
The first step in the validation process concerned format. In order to proceed with any 
other validation steps it was first necessary' to ensure that all data fields were in accord- 
ance with the specifications set forth in the user's manual of the thermal analyzer [Ref. 
9). The thermal analyzer input data file consists of five lines and as many as seven input 
data sets. The thermal analyzer's manual model builder program, THANSS, was in- 
strumental in verifying the accuracy of the model builder output data file format. 

The second step entailed content. The output of the model builder was checked for 
accuracy and completeness. Thermal resistances in the output data files were checked 
and compared with manually calculated benchmark models. A benchmark model was 
developed for each copper layer configuration. Other information such as node numbers 
and node relationships were also compared with the benchmark model and verified for 
accuracy. To aid in the verification process, the manual model builder program, 
THANSS, part of the thermal analyzer software was used to verify node relationships 
[Ref. 9J. Other items checked for correctness were specification codes relating to units 
(British or SI), temperature inputs, and thermal analyzer default values. 

The third and final step in the validation process required running the thermal ana- 
lyzer program with the PCB output data file as the program input. Abridged versions 
of the model builder were successfully executed with the current version of the thermal 
analyzer (maximum node capability is 300 nodes). The results of the thermal analysis 
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THIS IS THE 


OUTPUT OF THE 


THERMAL ANALYZER 








720 6 


0 


0 0 


0 


0 0 


1 




0 0 


0 












750 50 


6 


2 4 


6 


0 0 


0 




.0500000 . 


6GG6C70 


12 


.8000000 


87.00000 






66.000 


66.000 


66.000 64. 


.000 8G. 


000 85. 


000 


7 7551 


21 


7521 


131 


7511 


1211 




242.445 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 11 


31 


7521 


141 


7511 


1221 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 21 


41 


7521 


151 


7511 


1231 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 31 


51 


7521 


161 


7511 


1241 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 41 


61 


7521 


171 


7511 


1251 




67.280 


67.280 


41.249 


20.625 


.701 


.701 


2.50 


7 51 


71 


7521 


181 


7511 


12 61 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 61 


81 


7521 


191 


7511 


1271 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 71 


91 


7521 


201 


7511 


1281 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 81 


101 


7521 


211 


7511 


1291 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 91 


111 


7521 


221 


7511 


1301 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 101 


121 


7521 


231 


7511 


1311 




87.280 


87.280 


41.249 


20.625 


.701 


.701 


2.50 


7 111 


7541 


7521 


241 


7511 


1321 




87.280 


242.445 


41.249 


20.625 


.701 


.701 


2.50 


7 121 


141 


11 


251 


7511 


1331 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 131 


151 


21 


261 


7511 


1341 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 141 


161 


31 


271 


7511 


1351 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 151 


171 


41 


281 


7511 


1361 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 161 


181 


51 


291 


7511 


1371 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 171 


191 


61 


301 


7511 


1381 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 181 


201 


71 


311 


7511 


1391 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 191 


211 


81 


321 


7511 


1401 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 201 


221 


91 


331 


7511 


1411 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 211 


231 


101 


341 


7511 


1421 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 221 


241 


111 


351 


7511 


1431 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 231 


251 


121 


361 


7511 


1441 




67.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 241 


261 


131 


371 


7511 


1451 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 251 


271 


141 


381 


7511 


1461 




87.280 


87.280 


20.625 


20.625 


.701 


.701 


2.50 


7 261 


281 


151 


391 


7511 


1471 





Figure 9. Partial output data file of PCB terminal session. 



were compared with the results obtained using the analyzer's manual model builder, 
THANSS, and checked for discrepancies. 
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PlM He 



1 



this is TX OUTPUT OF TX WlH AHAlYTEP 
Ttaeer «tw r «« . In 



Node 

let 


1 

• s.oc 


• s.oc 


2 


s 

•s.oc 


4 

• S.OC 


1 

• S .00 


4 

• S.OC 


Hoe* 
let Tap 


7 

• S.OC 


•s.oc 


• 











hoot 
hi» Tap 
ht*-oid 


1 

• S.OC 
.7994 


2 

•s.e: 

.?»»» 


i 

os.oc 

.7MB 


4 

• S.OC 

.?m 


s 

•s.oc 

.7777 


4 

• S.OC 
.7 00 7 


hoot 
hew Tap 
hew-01 0 


7 

• S.OC 
.7007 


0 

•s.oc 

.7m 











hoot 
hew Tap 
hew-Clo 


1 

os.u 
• 1 4 CC 


2 

OS. 94 
. 14CC 


I 

• S .94 

.uc: 


4 

• S 94 

. 1 4 c: 


5 

OS. 94 

. 1 4 C C 


4 

• S . 9 4 
.1400 


Node 
he* Tap 
he w-c: fl 


7 

• 5.94 

. 14CC 


0 

• S . 9 4 

• 1 4 CC 










hede 
hew Tat 
hew-ClC 


1 

• 4 . 0 C 
.4C34E-0J 


2 

• 4.00 
.4C34E-C1 


I 

• 4 .0C 
.4 0 34 E-C 1 


4 

• 4 . ec 
.4041E-01 


s 

• 4.0C 
.4043E-C1 


4 

04. 0C 
.40A4E-01 


hefle 
hew Tap 

htw-Cifl 


7 

04. c: 
.4C45E-01 


O 

04 .e: 
.404AE-C1 











P»?f Nc. I 

TH10 IS THr OUTPUT OF The THERMAL ASAL V TER 
TEMPERATURES. CEGF 

1 ire. 78 : 10U.C6 3 1C5.30 A lCE.^A 5 10 6.06 6 lCt>.C3 

7 m.bQ 8 1 Oa . J 7 



Figure 10 . Output data file of TASS thermal analyzer. 
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E. MODEL BUILDER LIMITATIONS AND POSSIBLE IMPROVEMENTS 

Although the model builder provides the user with a significant advantage in terms 
of time and effort necessary to model a printed circuit board, there are at present certain 
limitations to the program. There is also significant potential for improvements that 
could enhance the capabilities of the current version of the model builder. 

1. PCB Model Builder Limitations 

1. The model builder is restricted to a maximum of four copper layers due to con- 
ventional personal computer operating system limitations. Current versions of 
MS-DOS ® do not allow the user to access extended memory. Access to extended 
memory would allow for a larger copper layer capability for the model. 

2. The model builder will not run when using MS-DOS ®, version 4.01. When run- 
ning PCB on a computer operating with MS-DOS ®, version 4.01, the computer 
will display an "insufficient memory" message. Writing more compact code for the 
model builder could overcome this problem. The program will run on any other 
version of MS-DOS ®. 

3. The model builder does not account for convective and radiative heat transfer. 

4. The model builder only accepts heat inputs from the top surface of the printed 
circuit board. 

2. Possible Improvements to the Thermal Model Builder. 

1. The present version of the model builder requires extensive amounts of TORTRAN 
code. An updated model builder could be designed using a more compact and ef- 
ficient higher level language, such as C++. Smaller, more efficient source code 
would allow for incorporating more copper layers into the model builder. 

2. The menu driven queries used in the thermal analyzer could be improved by in- 
corporating assembly language routines enabling the use of a mouse. Attempts 
were made to introduce the use of mouse driven menus; however, memory handling 
limitations in the Microsoft ® linker, version 3.55, prevented the use of a mouse. 

3. Allowing for the injection of heat not only from the top surface of the printed cir- 
cuit board, but also from other external surfaces (to include the heat being trans- 
ferred through the wedges supporting the printed circuit board), is another 
potential improvement. 

4. Adding graphics capabilities to the model builder would enable the user to model 
the printed circuit board without the need for preliminary drawings or sketches; 
however, adding this feature would cause the model builder executable code to be- 
come excessively large. 
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VI. CONCLUSIONS 



1 he purpose of this thesis is to develop a printed circuit board thermal analysis 
model builder that will effectively interact with thermal analysis software used by the 
Naval Postgraduate School. [Ref. 9] 

The current version of the model builder enables the user to choose from sixteen 
possible printed circuit board configurations, ranging from one to four copper layers. 
Other features of the model builder include the ability to set the physical characteristics 
of the board including dimensions, thermal conductivities, percent copper coverage, 
board temperatures, and heat input modes. 

The main goal of the model builder is to enable the user to model a printed circuit 
board with minimum effort, and produce a thermal analysis input data file. Using the 
model builder, PCB, relieves the user of the tedious, time-consuming manual data entry 
required of the thermal analyzer's current model builder, THANSS. [Ref. 9] 

The expansion possibilities for the model builder are significant. Adding more cop- 
per layers and nodes would enable the user to model more sophisticated printed circuit 
boards. Other potential features include the addition of graphics capabilities and the 
introduction of user-friendly peripherals, such as mouse drivers or digitazing pads. As 
the level of complexity of future model builders increases, it will become necessary to 
import the source code from FORTRAN compilers limited to 640 Kilobytes of memory, 
to 32 bit FORTRAN compilers allowing use of all RAM. The source code for PCB 
contains over 5000 lines of code, effectively reaching the memory limits of both the 
compiler and personal computer. 

The possible applications for the model builder are substantial. As the level of 
complexity of printed circuit boards increases, there is a real need for a tool that allows 
the designer to efficiently perform thermal analysis of printed circuit boards during the 
design process. 1 he model builder employed must allow the designer the necessary 
flexibility needed to model a board that meets the desired design criteria. The PCB 
model builder is the first step in creating a tool that enables the designer to effectively 
build a model that meets specified requirements. 
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APPENDIX 



PROGRAM LISTING FOR PCB MODEL BUILDER 



$ LARGE 



TITLE: MODEL BUILDER -- MAIN PROGRAM 

AUTHOR: LT STEVE GLASER 

DATE: 09 JUL 1991 

COMPILER: MICROSOFT VERSION 4. 01 

LINKER: MICROSOFT VERSION 3. 55 

DEFINE REAL VARIABLES 



REAL EL,EW,DELE,UKE,CL,CW,T1L1,T2L1,T2L2,T3L1,T3L2,T3L3,T4L1,T4L2, 
+T4L3 , T4L4 , T 1W 1 , T2W 1 , T2W2 , T3W 1 , T3W2 , T3W3 , T4W 1 , T4W2 , T4W3 , T4W4 , ARE A , U 
+N , VOL 1 , V0L2 , V0L3 , V0L4 , UKC , IBT , UPRT , LWRT , RT , LT , FT , BT , THEAT , THPN , AHE 
+AT , NHEAT , DELX , DELY 

COMMON I BT, UPRT, LWRT, RT,LT, FT, BT,UT, ANSWER 
INTEGER NWIDE,NDEEP,NPL,NUM,NUMA,CH,H 



DEFINE ALL ONE CHARACTER VARIABLES 



CHARACTER* 1 SELECT , ANS , ANSL , OVR , ANSWER , ANSN , ANSA , LAY 1 , LAY2 , LAY3 , LA 
+Y4 , ANSE , ANSEL , ANSEW , ANSET , ANSEK , ANSQ , WORL , ANSE 1 , ANST1 , ANSK 1 , ANSQ 1 , 
+SPEVAL , ANSTEM , ANSTEB , ANSTEU , ANSTEL , ANSTER , ANSLFT, ANSFRT , ANSBCK 



CHARACTER VARIABLES OF MORE THAN ONE POSITION 



CHARACTER UK*11,UT*1,UH*13,UAH*6,UL*2,UW*2 
DEFINE MATRICES 

REAL HEAT( 100, 100) ,COEF( 1000, 12) 

INTEGER IH( 1000) , JH( 1000) ,NC0N( 1000 , 12) 
PROVIDE THE USER WITH AN INTRODUCTORY STATEMENT. 
88888 CALL CLS 
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WRITE(* , 7001) 

7001 F0RMAT( ///////,' ' ,/, 

+' THIS PROGRAM WAS WRITTEN TO INTEGRATE WITH EXISTING ' ,/, 

+' THERMAL ANALYSIS SOFTWARE AND TO REDUCE THE AMOUNT ' ,/, 

+ ' OF TIME REQUIRED FOR DATA ENTRY. ' ,//, 

+' WOULD YOU LIKE AND OVERVIEW OF THE PROGRAM PRIOR TO 

+' BEGINNING? ENTER Y FOR YES AND N FOR NO: ' ,2X, ) 

READ(* , 7002) OVR 

7002 FORMAT(Al) 



VARIABLE, CONSTANT, AND STRING DEFINITION 
PHYSICAL CHARACTERISTICS 

CL, EL - EPOXY AND COPPER LENGTHS 
CW.EW - EPOXY AND COPPER WIDTHS 

UKE.UKC - EPOXY AND COPPER THERMAL CONDUCTIVITIES 

DELE, T1L1 THROUGH T4L4 - THICKNESS OF EPOXY AND COPPER LAYERS 

DELX - CL/NDEEP 

DELY - CW/NWIDE 

NPL - NUMBER OF NODES PER LAYER 
NWIDE - NUMBER OF NODES WIDE 
NDEEP - NUMBER OF NODES DEEP 
UL - UNITS OF LENGTH (SI OR ENGLISH) 

UK - UNITS OF THERMAL CONDUCTIVITY (SI OR ENGLISH) 

INITIAL AND AMBIENT TEMPERATURES 

IT - INITIAL BOARD TEMPERATURE 
LT - LEFT SIDE AMBIENT TEMPERATURE 
RT - RIGHT SIDE AMBIENT TEMPERATURE 
FT - FRONT AMBIENT TEMPERATURE 
BT - BACK AMBIENT TEMPERATURE 
UPRT - UPPER AMBIENT TEMPERATURE 
LWRT - LOWER AMBIENT TEMPERATURE 

UT - UNITS OF TEMPERATURE (CENTIGRADE OR FAHRENHEIT) 

HEAT INPUT 

THEAT - TOTAL INJECTED HEAT 

THPN - TOTAL HEAT PER NODE 

AHEAT - AVERAGE HEAT OVER A GIVEN SURFACE 

NHEAT - NODE PER HEAT INJECTED NODE PER NODE 

UH, UAH - UNITS OF HEAT (SI OR ENGLISH) 

NUM, NUMA, CH, H - DUMMY VARIABLES 

TOTNOD, NN, NC, N - VARIABLES USED TO ALLOW NODAL HEAT INPUT 
IH, JH - VECTORS USED TO CORRELATE NODE NUMBER WITH 
MATRIX POSITION 

HEAT - MATRIX USED TO HOLD HEAT INPUTS 



COEFFICIENT DEFINITIONS Y IMPLIES WIDTH, AND X IMPLIES DEPTH 

Z IMPLIES HEIGHT 

C IMPLIES COPPER, E IMPLIES EPOXY 
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EYLR - PROVIDES COEFFICIENT FOR LEFT OR RIGHT EDGE NODES 
TO THE EXTERNAL NODE, EYLR IMPLIES EPOXY. 

CYLR - SAME AS ABOVE, COPPER LAYER 

EXFB - PROVIDES COEFFICIENT FOR FRONT AND BACK EDGE NODES 
TO THE EXTERNAL NODE. 

CXFB - SAME AS ABOVE FOR COPPER LAYER 

EYY - INTERNAL COEFFICIENT IN THE Y DIRECTION, EPOXY 

CYY - INTERNAL COEFFICIENT IN THE Y DIRECTION, COPPER 

EXX - INTERNAL COEFFICIENT IN THE X DIRECTION, EPOXY 

CXX - INTERNAL COEFFICIENT IN THE X DIRECTION, COPPER 

SZE - COEFFICIENT FOR EPOXY TO UPPER EXTERNAL USE 
EZC - COEFFICIENT FOR COPPER TO EPOXY 
PCC - PERCENT COPPER COVERAGE 

ALL ONE CHARACTER STRINGS ARE SIMPLE 'YES' OR 'NO' TYPE 
VARIABLES. 

CALL INTRODUCTION 



IF (OVR. EQ. 'Y' ) THEN 
CALL INTRO 

ELSEIF (OVR. EQ. ’N’) THEN 

GOTO 899 

ELSE 

GOTO 88888 
END IF 

SUPPLY A LIST OF ACCEPTABLE LAYER LEVELS 
BEGIN PROGRAM 



899 

800 



C 

801 



CALL CLS 
WRITE (*,800) 
FORMAT( /// , ' 



+**** * / 
+ ' 

+ ’ 

+ ' 

+ ' 

+ ' 

+ ' 

+ ’ 



************************************************* ' , / , 



** ** ' j / } 

** PLEASE SELECT UPPER CASE LETTERS **',/, 
** PRIOR TO BEGINNING **',/, 
** **',/, 



A************************************************ ' , / , 



WRITE(*,801) 
FORMAT( // , 1 
+' PRIOR TO 

+' YOU HAVE 

+' DATA. 

+ ' 



ENTERING DATA INTO THIS PROGRAM ENSURE THAT ',/, 
A DRAWING OF YOUR DESIGN AND ALL PERTINENT ' ,/, 

’,///, 

PRESS <ENTER> TO CONTINUE ' ,2X, ) 
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READ(*, 802) ANSWER 
802 FORMAT(Al) 

C 

C 

C 

CALL CLS 

357 WRITE(*,367) 

367 FORMAK ///////, ' THIS PROGRAM IS CAPABLE OF OPERATIONS IN El 

+THER SI OR ’ ,/, 

+' ENGLISH UNITS. AFTER THE SELECTION OF THE UNITS, ALL 

+ ' ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION. ’/ 

+/,’ S FOR SI NOTATION’,/, 

+' E FOR ENGLISH NOTATION: ’,2X, ) 

READ(*,368) ANSN 

368 FORMAT(Al) 

C CHECK FOR CORRECT UNIT SELECTION 
C 

3333 IF(ANSN.EQ. 'S' ) THEN 
WRITE(*, 369) 

369 FORMAT( ////,' YOU HAVE SELECTED SI NOTATION.’) 

ELSEIF (ANSN.EQ. ’E’) THEN 

WRITE(*,370) 

370 FORMAT( //// , ’ YOU HAVE SELECTED ENGLISH NOTATION. ') 

ELSE 

CALL CLS 
GOTO 357 
END IF 
C 

3711 WRITE(*,3710) 

3710 FORMAT( // , ’ IS THIS THE DESIRED SELECTION? ENTER Y FOR YES A 

+ND N FOR NO: ' , ) 

READ(*,372) ANSA 
372 FORMAT(Al) 

C 

IF( ANSA. EQ. 'Y') THEN 
GOTO 378 

ELSEIF( ANSA. EQ. 'N') THEN 
CALL CLS 
GOTO 357 
ELSE 

CALL CLS 
GOTO 3333 
378 ENDIF 
C 
C 
C 
C 

7 CALL CLS 
SPEVAL='A' 

WRITE(*,301) 

301 FORMAT( //////,' THE PCB CAN HAVE UP TO FOUR COPPER LAYERS') 
WRITE(*,302) 

302 FORMAT( // , ' HOW MANY COPPER LAYERS DO YOU DESIRE?') 

WRITE(* , 303) 

303 FORMATC/,' PICK A NUMBER 1 TO 4: ' ,2X, ) 
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RE AD(*, 304) SELECT 

304 FORMAT(Al) 

C 

C 

CALL CLS 

IF ( SELECT. NE. ' 1 AND. SELECT. NE. ' 2 ' . AND. SELECT. NE. ’ 3 ' . AND. 

+SELECT. NE. '4' ) THEN 

GOTO 7 

ELSE 

10 WRITE(*, 305) SELECT 

305 FORMAT( //////, ' YOU SELECTED ' ,A1,' COPPER LAYER(S) FOR THE PCB 

+’,//) 

C 

ENDIF 

C 

WRITE(*,306) 

306 FORMAT ( ' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 

+R NO: ' ,2X, ) 

READ(*,307)ANS 

307 FORMAT(Al) 

C 

C 

IF( ANS. EQ. ’N') THEN 
GOTO 7 
ELSE 

CONTINUE 

ENDIF 

IF(ANS. EQ. 'Y') THEN 
GOTO 9 
ELSE 

CALL CLS 
GOTO 10 
9 ENDIF 
C 

c 

CALL CLS 
C 
C 
C 

IF( SELECT. EQ. ? 1* ) THEN 
GOTO 123 
ELSE 

CONTINUE 

ENDIF 

IF( SELECT. EQ. ’ 2 ' )THEN 
GOTO 124 
ELSE 

CONTINUE 

ENDIF 

IF( SELECT. EQ. ' 3 ' )THEN 
GOTO 125 
ELSE 

CONTINUE 

ENDIF 
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IF( SELECT. EQ. ' 4 1 )THEN 
GOTO 126 
ELSE 

CONTINUE 
END IF 
GOTO 7 
C 
C 

11 CALL CLS 

123 WRITE(*,308) 

308 FORMAT( //////, ' YOU SELECTED 1 COPPER LAYER, GIVING YOU THE FOLLO 
+WING ALTERNATIVES: ' ) 

WRITE(*,309) 

309 FORMAT(//,’ 1. 10 BY 36 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 

2. 10 BY 35 PROVIDES A 1:3.5 RATIO WITH 700 NODES. ',/, 

+' 3. 12 BY 25 PROVIDES A 1. 2: 2. 5 RATIO WITH 600 NODES. ',/, 

+' 4. 15 BY 20 PROVIDES A 1. 5: 2. 0 RATIO WITH 600 NODES.',//, 

+' PLEASE SELECT A NUMBER 1 THROUGH 4: ',2X, ) 

READ(*,310)LAY1 

310 FORMAT(Al) 

IF (LAY1.NE. 'l'.AND. LAY1.NE. ' 2 ' . AND. LAY1. NE. 

+ ’3 '.AND. LAY1.NE. '4') THEN 
GOTO 11 
ELSE 

3334 CALL CLS 

WRITE ( *,311) LAY1 

311 FORMAT( //////, ’ YOU SELECTED NUMBER ' A1 ' OF THE FOLLOWING 4 ALTER 
+NATIVES: ’,///, 

+ ' 1. 10 BY 36 PROVIDES A 1:3.6 RATIO WITH 720 NODES. ',/, 

+' 2. 10 BY 35 PROVIDES A 1:3.5 RATIO WITH 700 NODES. ',/, 

+ ' 3. 12 BY 25 PROVIDES A 1. 2: 2. 5 RATIO WITH 600 NODES.',/, 

+ ' 4. 15 BY 20 PROVIDES A 1. 5:2.0 RATIO WITH 600 NODES. ',///) 

ENDIF 

13 WRITE(* , 312) 

312 FORMAT(' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 

+R NO: ' ,2X, ) 

READ(*,313)ANSL 

313 FORMAT(Al) 

C 

C 

IF(ANSL. EQ. 'N') THEN 
GOTO 11 
ELSE 

CONTINUE 

ENDIF 

IF(ANSL. EQ. ' Y' ) THEN 
GOTO 12 
ELSE 

GOTO 3334 

12 ENDIF 
GOTO 1234 

14 CALL CLS 

124 WRITE(*,314) 

314 FORMAT( //////, ' YOU SELECTED 2 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ') 
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WRITE (*,315) 

315 FORMATC//,' 1. 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES. 

+ 2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES.',/, 

+' 3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES. ',/, 

+' 4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES. ',//, 

+' PLEASE SELECT A NUMBER 1 THROUGH 4: ' ,2X, ) 

READ(*,316)LAY2 

316 FORMATC Al) 

IF (LAY2.NE. ' 1 ' . AND. LAY2. NE. ' 2 ' . AND. LAY2. NE. 

+'3'.AND. LAY2.NE. '4' ) THEN 
GOTO 14 
ELSE 

3335 CALL CLS 

WRITE (*, 317 )LAY2 

317 FORMATC//////,' YOU SELECTED NUMBER 'Al' OF THE FOLLOWING 4 ALTER 
+NATIVES: ',///, 

+ ' 1. 10 BY 18 PROVIDES A 1:1.8 RATIO WITH 720 NODES.',/, 

+ ' 2. 9 BY 20 PROVIDES A 1:2.11 RATIO WITH 720 NODES. ',/, 

+' 3. 8 BY 20 PROVIDES A 2:5 RATIO WITH 640 NODES. ',/, 

+ ' 4. 10 BY 15 PROVIDES A 1:1.5 RATIO WITH 600 NODES. ',///) 

ENDIF 

15 WRITE(*,318) 

318 FORMATC' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 

+R NO: ' ,2X, ) 

READ(*,319)ANSL 

319 FORMATC Al) 

C 

C 

IF(ANSL.EQ. ' N' ) THEN 
GOTO 14 
ELSE 

CONTINUE 

ENDIF 

IFCANSL.EQ. 'Y' ) THEN 
GOTO 16 
ELSE 

GOTO 3335 

16 ENDIF 
GOTO 1235 

17 CALL CLS 

125 WRITE(*, 320) 

320 FORMATC//////,' YOU SELECTED 3 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ' ) 

WRITE(*, 321) 

321 FORMATC//,' 1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 

+ ',/,' 2. 10 BY 10 PROVIDES A 1:1 RATIO WITH 600 NODES. ',/, 

+' 3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES.',/, 

+' 4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES.',//, 

+' PLEASE SELECT A NUMBER 1 THROUGH 4: ' ,2X, ) 

READ(*,322)LAY3 

322 FORMATC Al) 

IF (LAY3.NE. ' 1 ' . AND. LAY3. NE. ' 2 ' . AND. LAY3. NE. 

+ '3'. AND.LAY3.NE. '4') THEN 
GOTO 17 
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ELSE 

CALL CLS 

WRITE (*,323)LAY3 

323 FORMATC//////,' YOU SELECTED NUMBER ' A1 ’ OF THE FOLLOWING 4 ALTER 
+NATIVES: ',///, 

+ ' 1. 5 BY 24 PROVIDES A 1:4.8 RATIO WITH 720 NODES. 

+' 2. 10 BY 10 PROVIDES A 1: 1 RATIO WITH 600 NODES.',/, 

+' 3. 12 BY 10 PROVIDES A 1.2:1 RATIO WITH 720 NODES. ' /, 

+' 4. 8 BY 15 PROVIDES A 1:1.875 RATIO WITH 720 NODES. ,///) 

ENDIF 

19 WRITE(* 324) 

324 FORMAT ( ' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 

+R NO: ' ,2X, ) 

READ(*,325)ANSL 

325 FORMAT(Al) 

C 

C 

IF( ANSL. EQ. 'N') THEN 
GOTO 17 
ELSE 

CONTINUE 

ENDIF 

IF(ANSL. EQ. ' Y' ) THEN 
GOTO 18 
ELSE 

GOTO 3336 
18 ENDIF 

GOTO 1236 

20 CALL CLS 

126 WRITE(*,326) 

326 FORMAT( //////,' YOU SELECTED 4 COPPER LAYERS, GIVING YOU THE FOLL 
+OWING ALTERNATIVES: ' ) 

WRITE(*,327) 

327 FORMATC//,' 1. 5 BY 18 PROVIDES A 1:3.6 RATIO WITH 720 NODES. 

+',/,' 2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES.',/, 

+ ' 3. 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES.',/, 

+' 4. 13 BY 6 PROVIDES A 1: 1 RATIO WITH 624 NODES. ’ ,//, 

+' PLEASE SELECT A NUMBER 1 THROUGH 4: ',2X, ) 

READ(*,328)LAY4 

328 FORMAT( Al) 

IF (LAY4.NE. ' 1 ' . AND. LAY4. NE. ' 2 ' . AND. LAY4. NE. 

+ '3 '.AND. LAY4.NE. '4') THEN 
GOTO 20 
ELSE 

3337 CALL CLS 

WRITE(*,329)LAY4 

329 FORMATC//////,' YOU SELECTED NUMBER 'Al' OF THE FOLLOWING 4 ALTER 

+' 1. 5 BY 18 'PROVIDES A 1:3.6 RATIO WITH 720 NODES.' /, 

+' 2. 9 BY 10 PROVIDES A 1:1.11 RATIO WITH 720 NODES. ,/, 

+ ' 3. 8 BY 10 PROVIDES A 1:1.25 RATIO WITH 640 NODES.',/, 

+' 4. 13 BY 6 PROVIDES A 1: 1 RATIO WITH 624 NODES. ',///) 

ENDIF 

21 WRITEC*,330) 

330 FORMATC IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N FO 

+R NO: ' ,2X, ) 
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READ(*,331)ANSL 
331 FORMAT(Al) 

C 

C 

IF(ANSL. EQ. 'N') THEN 
GOTO 20 

ELSE 

CONTINUE 

ENDIF 

IF( ANSL. EQ. ’ Y ' ) THEN 
GOTO 22 

ELSE 

GOTO 3337 
22 ENDIF 

GOTO 1237 

1234 CONTINUE 

IF(LAY1.EQ. ' 1 ' )THEN 

NWIDE=10 

NDEEP=36 

NPL=360 

ELSEIFCLAY1. EQ. ' 2' )THEN 
NWIDE=10 
NDEEP=35 
NPL=350 

ELSEIF(LAY1. EQ. ' 3* )THEN 
NVIDE=12 
NDEEP=25 
NPL=300 

ELSEIF(LAY1. EQ. '4’)THEN 
NVIDE=15 
NDEEP=20 
NPL=300 

ENDIF 

GOTO 1241 

1235 CONTINUE 

IF(LAY2. EQ. ' 1 1 )THEN 
NWIDE=10 
NDEEP=18 
NPL=180 

ELSEIF(LAY2.EQ. ' 2 1 )THEN 
NVIDE=9 
NDEEP=20 
NPL=180 

ELSEIF(LAY2. EQ. ' 3' )THEN 
NWIDE=8 
NDEEP=20 
NPL=160 

ELSEIF( LAY2. EQ. ' 4 ' )THEN 
NWIDE=10 
NDEEP=15 
NPL=150 

ENDIF 

GOTO 1241 

1236 CONTINUE 

IF(LAY3.EQ. 'l')THEN 
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NWIDE=5 

NDEEP=24 

NPL=120 

ELSEIF(LAY3. EQ. 1 2 1 )THEN 
NWIDE=10 
NDEEP=10 
NPL=100 

ELSEIF( LAY3. EQ. ' 3 ' )THEN 
NWIDE=12 
NDEEP=10 
NPL=120 

ELSEIF(LAY3. EQ. '4')THEN 
NWIDE=8 
NDEEP=15 
NPL=120 

ENDIF 

GOTO 1241 
1237 CONTINUE 

IF(LAY4. EQ. ' 1 ' )THEN 
NWIDE=5 
NDEEP=18 
NPL=90 

ELSEIF( LAY4. EQ. ' 2 ' )THEN 
NWIDE=9 
NDEEP=10 
NPL=90 

ELSEIF( LAY4. EQ. ' 3 1 )THEN 
NWIDE=8 
NDEEP=10 
NPL=80 

ELSEIF(LAY4. EQ. '4')THEN 
NWIDE=13 
NDEEP=6 
NPL=78 

ENDIF 

GOTO 1241 
C 

C EPOXY CHARACTERISTICS 

C 

1241 CALL CLS 

WRITE(*, 8000) 

8000 FORMAT ( / //, 1 Vf***yc******************yc**yc , *******Vr*** , ****** , ** , sV*yc* , ** , * , 3V 

+******************* f 3 / , 

+ 1 *****************Vc**********ePOXY CHARACTERI STICS******^***^^* 

+******** r a j 

+******** f , / 3 ) 
c 

C PROVIDE CORRECT UNIT ABBREVIATIONS 

C 

IF(ANSN. EQ. 'S') THEN 
WRITE(*,5400) 

5400 FORMATC ALL ENTRIES ARE IN SI NOTATION.',/) 

UL=' cm' 

UK=' Watts /cm/C ' 

UT=’C' 
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UW= ’ gm ' 

ELSEIF(ANSN. EQ. * E * ) THEN 
WRITE (*,5401) 

5401 FORMATC’ ALL ENTRIES ARE IN ENGLISH NOTATION’,/) 
UL=’ in’ 

UK=’ Btu/hr/F 1 
UT= 1 F ’ 

UW=’oz’ 

ENDIF 

C 

C 

c 



WRITE(*,5402) UL 

5402 FORMATC/,’ ENTER EPOXY LAYER LENGTH (’, A2 ,’) : ’ ,2X, ) 

READ * ,EL 

C 

WRITE(*,5403) UL 

5403 FORMATC/,’ ENTER EPOXY LAYER WIDTH (’ ,A2 ,’ ): ’,2X, ) 

READ *,EW 

4337 WRITE(*,6403) 

6403 FORMATC//,’ THE EPOXY LAYER THICKNESS WILL DEFAULT TO 0.0625 in ( 
+0.15875 cm )’,/,’ DO YOU WANT TO CHANGE THE EPOXY THICKNESS? (Y OR 
+ N): ’ ,2X, ) 

READ(*,5799) ANSQ 
5799 FORMATC Al) 

IFCANSQ. EQ. ’Y’ ) THEN 
GOTO 5798 

ELSEIF(ANSQ.EQ. ’ N ’ . AND. ANSN. EQ. ’ S ’ ) THEN 
DELE=0. 15875 
GOTO 4338 

ELSEIFCANSQ. EQ. ’ N’ . AND. ANSN. EQ. ’ E ’ ) THEN 
DELE=0. 0625 
GOTO 4338 
ELSE 

GOTO 4337 
ENDIF 

5798 WRITE(* ,5404) UL 

5404 FORMATC/,’ ENTER EPOXY LAYER THICKNESS (’ ,A2, ’ ): ’ ,2X, ) 

READ *,DELE 

C 

4338 WRITE(*,5405) UK 

5405 FORMATC/,’ ENTER EPOXY LAYER THERMAL CONDUCTIVITY C’.AIO,’): ’,2 

+X, ) 

READ *,UKE 
5459 CALL CLS 
C 

C MAKE CHANGES OR CORRECTIONS TO EPOXY ENTRIES 
C 



WRITE(*,5406) 

5406 FORMATC///,’ 

+ LAYER: ’ , / , ) 

5407 WRITEC*, 89343) 

89343 FORMATC/,’ 
WRITEC*, 89344) 

89344 FORMATC/,’ 



YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE EPOXY 
EL,UL 

1. ) LENGTH: ’ ,F9.4,1X,A2, ) 

EW.UL 

2. ) WIDTH: ’ ,F9. 4,1X,A2, ) 
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WRITEC*, 89345) DELE ,UL 

89345 FORMATC/,' 3.) THICKNESS: ' ,F9. 4, IX, A2, ) 

WRITEC*, 89346) UKE,UK 

89346 FORMATC/,' 4. ) k: ' ,F9. 4, IX, All , /// , ) 

C 

WRITEC*, 5408) 

5408 FORMAT(/,' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES 

+ AND N FOR NO: ' ,2X, ) 

READC* , 5409 ) ANSE 

5409 FORMATC Al) 

C 

C 

IFCANSE.EQ. 'Y') THEN 

5441 CALL CLS 
WRITEC*, 5410) EL,UL 

5410 FORMATC////,' THE CURRENT ENTRY FOR LENGTH IS: ' ,F9. 4 , IX, A2) 
WRITEC* , 5411) 

5411 FORMATC/,' WOULD YOU LIKE TO CHANGE THE LENGTH? CY OR N): 

+* 2X, ) 

READC *,54 12) ANSEL 

5412 FORMATC Al) 

PRINT * 

IFC ANSEL. EQ. ' Y r ) THEN 
WRITEC* ,5422) UL 

5422 FORMATC/,' ENTER THE EPOXY LENGTH C ' ,A2, ' ): ',2X, ) 

READ * , EL 

ELSEIFC ANSEL. EQ. 'N' ) THEN 

GOTO 5440 

ELSE 

GOTO 5441 
END IF 
C 

5440 CALL CLS 

WRITEC*, 5442) EW,UL 

5442 FORMATC////,' THE CURRENT ENTRY FOR WIDTH IS: ' ,F9. 4, IX, A2) 

WRITEC*, 5443) 

5443 FORMATC/,' WOULD YOU LIKE TO CHANGE THE WIDTH? CY OR N): 

+ ' , 2X , ) 

READ(* ,5444)ANSEW 

5444 FORMATC Al) 

PRINT * 

IFCANSEW. EQ. ' Y’ ) THEN 
WRITEC*, 5445) UL 

5445 FORMATC/,' ENTER THE EPOXY WIDTH C',A2,'): ' ,2X, ) 

READ * ,EW 

ELSEIFC ANSEW. EQ. 'N' ) THEN 

GOTO 5446 

ELSE 

GOTO 5440 
ENDIF 
C 

5446 CALL CLS 

WRITEC*, 5447) DELE ,UL 

5447 FORMATC////,' THE CURRENT ENTRY FOR THICKNESS IS: ' ,F9. 4, IX, A2) 

WRITEC*, 5448) 

5448 FORMATC/,' WOULD YOU LIKE TO CHANGE THE THICKNESS? CY OR N 



58 



+) : 1 2X , ) 

READ(*,5449)ANSET 

5449 FORMAT( Al) 

PRINT * 

IF(ANSET. EQ. 'Y' ) THEN 
WRITE(* ,5450) UL 

5450 FORMAT(/,' ENTER THE EPOXY THICKNESS (' ,A2, ' ): ' ,2X, ) 

READ *,DELE 

ELSEIF(ANSET.EQ. 'N') THEN 

GOTO 5451 

ELSE 

GOTO 5446 
ENDIF 
C 

5451 CALL CLS 

WRITE (*,545 2) UKE,UK 

5452 F0RMAT( ////,' THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS: ’, 
+F9. 4, 1X,A10) 

WRITE (*,5453) 

5453 FORMAT( / , ' WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCTIVIT 

+Y? (Y OR N): ' ,2X, ) 

READ(* ,5454)ANSEK 

5454 FORMAT(Al) 

PRINT * 

IF( ANSEK. EQ. ' Y* ) THEN 
WRITE(*,5455) UK 

5455 FORMAT( / , ' ENTER THE THERMAL CONDUCTIVITY (' ,A10 ,' ): ' ,2X, ) 

READ * ,UKE 

ELSEIF( ANSEK. EQ. 1 N ’ ) THEN 

GOTO 5456 

ELSE 

GOTO 5451 
ENDIF 

5456 CALL CLS 
WRITE(*,5457) 

5457 FORMAT( /// , ’ YOU HAVE MADE THE FOLLOWING CORRECTIONS TO THE 

+ EPOXY ENTRIES: ’ ,//) 

GOTO 5407 
C 
C 

ELSEIF(ANSE.EQ. ’N’) THEN 

GOTO 5458 

ELSE 

GOTO 5459 
ENDIF 
C 
C 
C 

C CALL SUBROUTINE COPPER 
C 

c 

c 

5458 CALL COPPER( EW , EL , ANSN , UK , SELECT , T1L1 , T2L 1 , T2L2 , T3L1 , T3L2 , T3L3 , T4L 
+ 1 , T4L2 , T4L3 , T4L4 , S PE VAL , UKC , CL , C W ) 

C 

C 
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oooono oooooo nnnnnnnnnnnnnnonn nnnno ooonoooooooon 



CALL SUBROUTINE PCBS1 



1180 CALL PCBS1 



CALL SUBROUTINE PCBS2 



CALL PCBS2 ( THE AT , THPN , AHEAT , NHEAT , NPL , NW IDE , NDEEP , HEAT , I H , JH , ANSN , 
+EL,EW) 



CALL SUBROUTINE PCBS3 



CALL PCBS3(EW ,EL,CW ,CL,UKE ,DELE ,UKC , SELECT, T1L1 ,T2L1 ,T2L2 ,T3L1 ,T3L 
+2 , T3L3 , T4L1 , T4L2 ,T4L3 , T4L4 , NWIDE , NDEEP , NPL , IH , JH , HEAT , COEF , IBT, UPR 
+T , LWRT , RT , LT , FT , BT , ANSN ) 

END 



SUBROUTINE INTRO 



TITLE: 

SUBROUT: 

AUTHOR: 

DATE: 

COMPILER: 

LINKER: 



MODEL BUILDER 
INTRO 

LT STEVE GLASER 
09 JUL 1991 

MICROSOFT VERSION 4. 01 
MICROSOFT VERSION 3. 55 



THIS SUBROUTINE PROVIDES THE USER WITH AN OVERVIEW OF THE THERMAL 
ANALYZER MODEL BUILDER. 



COMMON IBT , UPRT , LWRT , RT , LT , FT , BT , UT , ANSWER 



DEFINE ONE VARIABLE CHARACTER VALUES 
CHARACTER*! ANSWER 



PROVIDE THE USER WITH A PROGRAM OVERVIEW 
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7002 FORMAT (Al) 

CALL CLS 
WRITE(* , 7003) 

7003 FORMAT( //////, ' ********* OVERVIEW ********** 

+' THIS PROGRAM PERFORMS A NODAL ANALYSIS OF A PRINTED 
+' CIRCUIT BOARD CONTAINING UP TO FOUR COPPER LAYERS 

+' (WITH EPOXY LAYERS IN BETWEEN). THE OUTPUT CONSISTS 
+' OF UP TO 720 COEFFICIENTS THAT CONTRIBUTE TO THE 
+' DETERMINATION OF THE TEMPERATURE DISTRIBUTION 
+' OF THE PRINTED CIRCUIT BOARD (PCB) WHEN FED INTO THE 
+' THERMAL ANALYZER. ' ,//, 

+' THE FOLLOWING IS AN OUTLINE OF THE MAJOR ' ,/, 

+' SECTIONS OF THIS PROGRAM AND WHAT ENTRIES ARE 
+' REQUIRED OF THE USER. 

+ 1 

+' PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 
+ ' 

+' PRESS <ENTER> TO CONTINUE 

RE AD(*, 7002) ANSWER 
C 

CALL CLS 
C 

WRITE(*, 7004) 

7004 FORMAT( //////, ' ********* OVERVIEW ********** 

+' A. DATA OUTPUT FILE: 

+ ' 

+' THIS PROGRAM GENERATES AN OUTPUT DATA FILE WHICH IS 
+' TO BE THE INPUT TO THE THERMAL ANALYZER. 

+' DURING THE COURSE OF THE PROGRAM, THE USER WILL BE 
+' ASKED TO PROVIDE A NAME FOR THE OUTPUT FILE. 

+ ' 

+' WHEN PROMPTED PLEASE ENTER THE NAME OF THE OUTPUT 
+' DATA FILE. THE DATA FILE NAME SHOULD BE NO LONGER 
+' THAN SIX LETTERS, AND MAY NOT HAVE ANY SPACES. 

+ 1 

+' PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 
+ ' 

+' PRESS <ENTER> TO CONTINUE 

RE AD(*, 7 002) ANSWER 
C 



f 




,// 



y 






r 



t 

> 



2X, ) 



' >/, 

' ,/, 

’,2X, ) 



7005 



CALL CLS 
WRITE (*,7005) 

FORMAT( //////,' ********* OVERVIEW ********** 

+' B. STRUCTURE PHYSICAL CHARACTERISTICS 
+ ' 

+' 1. YOU WILL BE ASKED TO SELECT UNIT TYPE 

+' (SI OR ENGLISH). 

+' 2. THE PRINTED CIRCUIT BOARD IS MAINLY COMPOSED OF 

+' ALTERNATING COPPER AND EPOXY LAYERS. THE PROGRAM 

+' IS DESIGNED TO PROVIDE THE USER WITH A RANGE 

+' OF 1 TO 4 COPPER LAYERS. IT IS ASSUMED THAT EACH 

+' COPPER LAYER LIES BETWEEN EPOXY LAYERS. 

+' DEPENDING ON THE NUMBER OF COPPER LAYERS SELECTED 

+' THE PROGRAM WILL ALLOW THE USER TO CHOOSE FROM 



f 



,/ 

,/ 

J 

J 



/, 

/, 

/, 

/, 




,// 



y 
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7006 



7007 



7008 



+' FOUR ALTERNATIVE NODAL ASPECT RATIOS. 

+ ' 

+' PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 
+ ' 

+' PRESS <ENTER> TO CONTINUE 

READ(*,7002)ANSWER 
CALL CLS 
WRITE(*, 7006) 

FORMAT ( //////, ? ********* OVERVIEW ********** 



■;/: 



,2X, ) 



+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 



B. 

3. 



STRUCTURE PHYSICAL CHARACTERISTICS (CONTINUED) 



y / y 
/, 



r 



AFTER SELECTING THE DESIRED NUMBER OF COPPER 
LAYERS AND NODAL ASPECT RATIO, THE PROGRAM WILL 
THEN ASK THE USER TO PROVIDE LAYER CHARACTERISTICS. ',/, 
EPOXY AND COPPER LENGTH, THICKNESS, AND WIDTH, ',/, 
AS WELL AS CONSTANTS CONSTITUTE THE QUERIES. ',/, 

THE PROGRAM WILL ALSO ASK THE USER TO PROVIDE ' ,/, 

THE PERCENT COVERAGE FOR EACH COPPER LAYER. ' , / . 

,/> 

PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 



PRESS <ENTER> TO CONTINUE 
READ(* , 7002) ANSWER 
CALL CLS 
WRITE(*, 7007) 

FORMATC //////, ' ********* OVERVIEW ********** 

+' C. INTITIAL AND AMBIENT TEMPERATURES 
+ ' 

+ ' 1 . 

+ ' 

+ ' 

+ ' 

+ ' 

+ ' 

+ ' 

+ ' 

+ ' 

+ ' 



V’ 



AFTER SELECTING THE DESIRED NUMBER OF COPPER 
LAYERS AND NODAL ASPECT RATIO, THE PROGRAM WILL 



,2X, ) 



• 

y / y 



THEN ASK THE USER TO PROVIDE TEMPERATURES FOR THE . ’,/, 
BOARD. ',/, 

PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS ',/, 

\2X, 



PRESS <ENTER> TO CONTINUE 
READ(*, 7 002) ANSWER 
CALL CLS 
WRITE(*,7008) 

FORMAT(//////,' ********* OVERVIEW ********** 



+ ’ 
+ ' 
+ ' 
+ ' 
+ • 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ' 
+ ’ 
+ ’ 
+ ' 



D. HEAT INPUT. 

1. HEAT INJECTION OCCURS ONLY ON THE UPPER COPPER 
LAYER. THIS PROGRAM SUPPLIES THE USER FOR 
ALTERNATIVE METHODS FOR ENTERING HEAT. 



A. ) TOTAL HEAT OVER SURFACE. 

B. ) AVERAGE HEAT PER UNIT AREA 

C. ) INPUT HEAT NODE BY NODE 

D. ) NO HEAT INPUT 

PLEASE NOTE: ENTRIES MUST BE IN UPPER CASE LETTERS 
****THIS CONCLUDES THE PROGRAM OVERVIEW**** 



' >/, 



\ll 

' J, 

' ,/, 

f/ / 
y / y 

' >/, 



',//, 
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' ,2X, ) 



C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



+ ’ PRESS <ENTER> TO CONTINUE 

READ(*, 7002) ANSWER 
CALL CLS 
END 

SUBROUTINE PCBS1 

********************************************************************* 



TITLE: MODEL BUILDER 

SUBROUT: PCBS1 

AUTHOR: LT STEVE GLASER 

DATE: 09 MAY 1991 

COMPILER: MICROSOFT VERSION 4.01 

LINKER: MICROSOFT VERSION 3. 55 

THIS SUBROUTINE HANDLES THE TEMPERATURE INPUTS TO THE PRINTED CIRCUIT 
BOARD. 



COMMON I BT , UPRT , LWRT , RT , LT , FT , BT , UT , ANSWER 
DEFINE REAL VARIABLES 



REAL I BT , UPRT , LWRT , RT , LT , FT , BT 
DEFINE ONE CHARACTER VARIABLES 



CHARACTER* 1 ANSTEM , ANSTEB , ANSTEU , ANSTEL , ANSTER , ANSLFT , ANSFRT , ANSBC 
+K,UT 
C 

31180 CALL CLS 

WRITE(*, 60000) 

60000 FORMATC III ,' ****************************************************** 
******************** ' , 

+’ *******************^bient temperature inputs******************** 

4 .******** 1 j 

+ 1 ***************************************************************** 

+******** ’,//,) 
c 
c 

WRITE(*, 60001) UT 

60001 FORMATC/,’ ENTER THE INITIAL BOARD TEMPERATURE ( ’ , A1 , ’ ) : ’ ,2X, 

+) 

READ *,IBT 
C 
C 

WRITEC*, 60002) UT 

60002 FORMATC/,’ ENTER THE UPPER SURFACE AMBIENT TEMPERATURE (',A1,') 
+: \2X, ) 

READ * , UPRT 
C 
C 
C 
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WRITE(*, 60003) UT 

60003 FORMAT( / , ' ENTER THE LOWER SURFACE AMBIENT TEMPERATURE ( ' ,A1, ’ ) 
+: ' , 2X, ) 

READ * , LWRT 
C 
C 

WRITEC*, 60004) UT 

60004 FORMAT( / , ' ENTER THE RIGHT SURFACE AMBIENT TEMPERATURE C'.Al,’) 
+: ' ,2X, ) 

READ *,RT 
C 
C 

c 



WRITEC*, 60005) UT 

60005 FORMAT( / , ' ENTER THE LEFT SURFACE AMBIENT TEMPERATURE ( ’ ,A1, ’ ): 
+ \2X, ) 

READ *,LT 
C 
C 

WRITEC*, 60006) UT 

60006 FORMAT( / , ' ENTER THE FRONT SURFACE AMBIENT TEMPERATURE C',A1,') 
+: ’ > 2X, ) 

READ *,FT 
C 
C 

c 



WRITEC*, 60007) UT 

60007 FORMAT( / , ' ENTER THE REAR SURFACE AMBIENT TEMPERATURE (',A1,'): 
+ ' j 2X, ) 

READ *,BT 



REVIEW THE TEMPERATURE ENTRIES 



C 
C 
C 
C 

60052 CALL CLS 

WRITEC*, 60008) IBT,UT 

60008 FORMATC////,’ YOU HAVE MADE THE FOLLOWING AMBIENT TEMPERATURE E 

+NTRIES: ',////,’ 1. INITIAL BOARD TEMPERATURE: ' , 1X,F9. 3, 1X,A1) 

WRITE(*, 60009) UPRT.UT 

60009 FORMAT( ’ 2. UPPER AMBIENT TEMPERATURE: 

WRITE(*, 60010) LWRT , UT 

60010 FORMATC' 3. LOWER AMBIENT TEMPERATURE: 

WRITEC*, 60011) RT.UT 

60011 FORMATC' RIGHT AMBIENT TEMPERATURE: 

WRITEC*, 60012) LT,UT 

60012 FORMATC' 5. LEFT AMBIENT TEMPERATURE: ' , IX, F9. 3, IX, Al) 

WRITEC*, 60013) FT,UT 

60013 FORMATC’ 6. FRONT AMBIENT TEMPERATURE: ' , IX, F9. 3 , IX, Al) 

WRITEC*, 60014) BT,UT 

60014 FORMATC' 7. REAR AMBIENT TEMPERATURE: ' , IX, F9. 3 , IX, Al , //) 

C 
C 
C 



’ ,1X,F9. 3,1X,A1) 
' ,1X,F9. 3,1X,A1) 
' ,1X,F9. 3,1X,A1) 



WRITEC *,600 15) 

60015 FORMATC/,' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES A 
+ND N FOR NO: ' ,2X, ) 
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READC*, 600 16) ANSTEM 

60016 FORMAT(Al) 

C 

C 

C MAKE CORRECTIONS BOARD TEMPERATURES 

C 

c 

IFC ANSTEM. EQ. 'Y' ) THEN 

60022 CALL CLS 

WRITEC*, 60017) IBT,UT 

60017 FORMATC /// , ' THE INITIAL BOARD TEMPERATURE IS: ' ,F9. 3 , IX, Al) 

WRITEC*, 60018) 

60018 FORMATC/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X, ) 

READ(*, 60019) ANSTEB 

60019 FORMATC Al) 

IFC ANSTEB. EQ. 'Y') THEN 
WRITEC*, 60020) UT 

60020 FORMATC/,' ENTER THE NEW VALUE C ' ,A1, ' ): ',2X, ) 

READ *, IBT 

ELSEIFC ANSTEB. EQ. 'N' ) THEN 
GOTO 60021 
ELSE 

GOTO 60022 
ENDIF 
C 

60021 CALL CLS 

WRITEC*, 70022) UPRT,UT 

70022 FORMATC///,' THE UPPER AMBIENT TEMPERATURE IS: ' ,F9. 3 , IX, Al) 
WRITEC*, 60023) 

60023 FORMATC/,’ WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 

+X, ) 

READC*, 60024) ANSTEU ' 

60024 FORMATC Al) 

IF( ANSTEU. EQ. 'Y' ) THEN 
WRITEC*, 60025) UT 

60025 FORMATC/,’ ENTER THE NEW VALUE C ', A 1 ,’) : ',2X, ) 

READ *, UPRT 

ELSEIFC ANSTEU. EQ. 'N') THEN 
GOTO 60026 
ELSE 

GOTO 60021 
ENDIF 
C 
C 

60026 CALL CLS 

WRITEC*, 60027) LWRT,UT 

60027 FORMATC///,' THE LOWER AMBIENT TEMPERATURE IS: ’ ,F9. 3 , IX, Al) 
WRITEC*, 60028) 

60028 FORMATC/,’ WOULD YOU LIKE TO CHANGE THIS VALUE? CY OR N): ’,2 

+X, ) 

READC*, 60029) ANSTEL 

60029 FORMATC Al) 

IFC ANSTEL. EQ. 'Y') THEN 
WRITEC*, 60030) UT 

60030 FORMATC/,' ENTER THE NEW VALUE C ’ ,A1 ,’ ): ' ,2X, ) 
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READ *,LWRT 

ELSEIF( ANSTEL. EQ. ’ N ' ) THEN 
GOTO 60031 
ELSE 

GOTO 60026 
END IF 
C 
C 
C 

60031 CALL CLS 

WRITEC*, 60032) RT,UT 

60032 FORMATC///,' THE RIGHT AMBIENT TEMPERATURE IS: ' ,F9. 3 , 1X,A1) 
WRITEC*, 60033) 

60033 FORMATC/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 
+X, ) 

RE ADC*, 60034) ANSTER 

60034 FORMAT(Al) 

IF( ANSTER. EQ. ' Y* ) THEN 
WRITEC*, 60035) UT 

60035 FORMAT(/ , 1 ENTER THE NEW VALUE (’ ,A1, ' ): ' ,2X, ) 

READ *,RT 

ELSEIF( ANSTER. EQ. 'N') THEN 
GOTO 60036 
ELSE 

GOTO 60031 
ENDIF 
C 
C 
C 

c 

60036 CALL CLS 

WRITEC*, 60037) LT,UT 

60037 FORMATC///,' THE LEFT AMBIENT TEMPERATURE IS: ’ ,F9. 3,1X,A1) 
WRITEC*, 60038) 

60038 FORMATC/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 

+X, ) 

READC*, 60039) ANSLFT 

60039 FORMAT(Al) 

IF( ANSLFT. EQ. ’ Y* ) THEN 
WRITEC*, 60040) UT 

60040 FORMATC/,' ENTER THE NEW VALUE (' ,A1, ' ): ',2X, ) 

READ * , LT 

ELSEIF(ANSLFT. EQ. 'N') THEN 
GOTO 60041 
ELSE 

GOTO 60036 
ENDIF 
C 
C 

c 

c 

60041 CALL CLS 

WRITE(*, 60042) FT,UT 

60042 FORMATC///,' THE FRONT AMBIENT TEMPERATURE IS: ' ,F9. 3 , IX, Al) 
WRITE(*, 60043) 

60043 FORMATC/,' WOULD YOU LIKE TO CHANGE THIS VALUE? CY OR N): ',2 
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+x, ) 

READ(*, 60044) ANSFRT 

60044 FORMATC Al) 

IF( ANSFRT. EQ. 'Y' ) THEN 
WRITE( ,V , 60045) UT 

60045 FORMAT( / , ' ENTER THE NEW VALUE (’, A 1 , ') : ' ,2X, ) 
READ *,FT 

ELSEIF( ANSFRT. EQ. 'N' ) THEN 
GOTO 60046 
ELSE 

GOTO 60041 
ENDIF 



C 

C 

c 

c 

60046 CALL CLS 

WRITEC*, 60047) BT,UT 

60047 FORMATC///,' THE REAR AMBIENT TEMPERATURE IS: ' ,F9. 3 , IX, Al) 
WRITEC*, 60048) 

60048 FORMATC/,' WOULD YOU LIKE TO CHANGE THIS VALUE? (Y OR N): ',2 

+X, ) 

READ(*, 60049) ANSBCK 

60049 FORMATC Al) 

IFCANSBCK. EQ. 'Y' ) THEN 
WRITEC*, 60050) UT 

60050 FORMATC/,' ENTER THE NEW VALUE (' ,A1 ,' ): ' ,2X, ) 

READ *,BT 

ELSEIFC ANSBCK. EQ. 'N') THEN 
CONTINUE 
ELSE 

GOTO 60046 
ENDIF 

ELSEIFC ANSTEM. EQ. 'N') THEN 
GOTO 60051 
ELSE 

GOTO 60052 
ENDIF 

60051 END 

SUBROUTINE PCBS2CTHEAT,THPN,AHEAT,NHEAT,NPL,NWIDE,NDEEP,HEAT, IH, JH 
+,ANSN,EL,EW) 



C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



TITLE: MODEL BUILDER 

SUBROUT: PCBS2 -- HEAT INPUT SUBROUTINE 

AUTHOR: LT STEVE GLASER 

DATE: 09 MAY 1991 

COMPILER: MICROSOFT VERSION 4.01 

LINKER: MICROSOFT VERSION 3. 55 

THIS SUBROUTINE HANDLES THE HEAT INPUTS TO THE PRINTED CIRCUIT BOARD 
DEFINE REAL VARIABLES 



REAL THE AT , THPN , AHE AT , NHE AT , S L , SW 
C 
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c 

C DEFINE INTEGERS 

C 

INTEGER NWIDE ,NDEEP,NPL,NUM,NUMA,CH,H,TOTNOD ,NC,NN, I , J, IM 
C 
C 

C DEFINE ONE CHARACTER VARIABLES 

C 

CHARACTERS SELH,ANSH,ANSN,ATH,ANSHA,AHN, DUMMY 
C 

C DEFINE CHARACTER VARIABLES 

C 

CHARACTER UH*13,UAH*6 
C 

C DEFINE REAL MATRICES 

C 

REAL HEAT( 100 ,100) ,COEFF( 740 ,9) 

C 

C DEFINE INTEGER VARIABLES 

C 

INTEGER IH( 1000), JH( 1000), UMMY( 1000) 

C 

C PROVIDE A CORRELATION BETWEEN NODE NUMBERS AND MATRIX LOCATION 

C 

C 

51 NUM=1 

DO 60 1=1 ,NPL/NWIDE 
DO 61 J=l, NWIDE 
JH(NUM)=J 
NUM=NUM+1 

61 CONTINUE 

60 CONTINUE 

C 

NUMA=1 

CH=0 

H=1 

DO 62 1=1 ,NPL 
IH(NUMA)=H 
CH=CH+1 

IF(CH.EQ. NWIDE) THEN 
H=H+1 
CH=0 
ELSE 

CONTINUE 
END IF 

NUMA=NUMA+1 

62 CONTINUE 
C 

c 

c 

1241 CALL CLS 

WRITE (*,8000) 

8000 FORMAT( /// , 1 ****************************************************** 
******************** 1 5 5 

+ ' *******************************he AT I NPUTS^*********^**’'”'^^** 

********* 1 } j 
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+ * ***************************************************************** 

+******** 1 , / / , 

+' HEAT’ INPUT TO THE PCB OCCURS ONLY ON THE UPPER COPPER',/, 

+’ LAYER. HEAT INPUT IS ACCOMPLISHED BY ONE OF THE ',/, 

+' FOLLOWING METHODS: ',//, 

+' 1. ENTER AS TOTAL HEAT APPLIED TO THE PCB',/, 

+' 2. ENTER AS AVERAGE HEAT PER UNIT AREA ',/, 

+' 3. ENTER HEAT NODE BY NODE ' ,/, 

+ ' 4. NO HEAT INPUT' ,//, 

+' PLEASE SELECT NUMBER 1 THROUGH 4: ' ,2X, ) 

C 

READ(*,5799) SELH 
5799 FORMAT(Al) 

IF( SELH. EQ. ' 1 ' . OR. SELH. EQ. ' 2 ' . OR. SELH. EQ. ' 3 ' . OR. SELH. EQ. ' 4 ' ) THEN 
444 WRITE(*,291) SELH 

291 FORMATC//,' YOU HAVE SELECTED NUMER ' ,A1,' OF FOUR ALTERNATIVES 
+.',/,' IS THIS THE DESIRED SELECTION? (Y OR N): ' ,2X, ) 

ELSE 

GOTO 1241 
ENDIF 

READ(*,5409)ANSH 
5409 FORMATC Al) 

C 

C 

IF(ANSH. EQ. 'Y') THEN 
5441 GOTO 63 

ELSEIF( ANSH. EQ. 'N') THEN 
GOTO 51 
ELSE 

CALL CLS 
WRITE(* , 8000) 

GOTO 444 
63 ENDIF 

C 
C 

C DETERMINE UNIT FOR HEAT INPUT 
C 

IF ((ANSN.EQ. 'E'). AND. (SELH.EQ. '2')) THEN 
UH = 'Btu/(hr*in 2)' 

UAH = ' Btu/hr ' 

ELSEIF ((ANSN.EQ. ' E '). AND. ( SELH. EQ. ' 1 ' . OR. SELH. EQ. '3')) THEN 
UH = 'Btu/hr' 

ELSEIF ((ANSN.EQ. 'S' ). AND. (SELH.EQ. '2' )) THEN 
UH = 'Watts/(cm 2)' 

UAH = 'Watts* 

ELSEIF ((ANSN.EQ. ’S'). AND. (SELH.EQ. '1'. OR. SELH.EQ. '3')) THEN 
UH = 'Watts' 

ENDIF 

C 

C 

C ALLOW FOR RE-SELECTION OF HEAT INPUT METHOD OR CONTINUE WITH 
C INITIAL SELECTION. 

C 

C CHOICE # 1 

c 

c 
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CALL CLS 

IF (SELH. EQ. ' l'.AND. ANSN. EQ. ' E ’ ) THEN 

989 WRITE (*,988) UH 

988 FORMAT(///,' YOU HAVE SELECTED TO INPUT HEAT AS A TOTAL HEAT 

APPLIED TO THE SURFACE. 

+' ENTER TOTAL HEAT APPLIED TO THE SURFACE (',A6,'): ',2X, ) 

READ * , THEAT 

990 WRITE(*,991) 

991 FORMAT( / , ' IS THIS THE CORRECT ENTRY? (Y OR N): ' ,2X, ) 

READ(*,992)ATH 

992 FORMAT(Al) 

C 

C MAKE ENTRY AND ALLOW FOR CORRECTIONS 
C 

IF (ATH. EQ. ' Y’ ) THEN 
THPN=THEAT/NPL 
WRITE(*,993) THPN,UH 

993 FORMAT( / , ' TOTAL HEAT PER NODE IS: ' ,F9. 4 , IX, A13) 

WRITE(* , 9323) 

9323 FORMAT(//,' <PRESS ENTER TO CONTINUE^) 

READ(*,9324)DUMMY 

9324 FORMAT(Al) 

C 

C 

C FILL HEAT MATRIX WITH DESIRED VALUES 

C 

C 

DO 994 1=1 ,NDEEP 
DO 995 J=1 ,NWIDE 
HEAT( I , J)=THPN 

995 CONTINUE 

994 CONTINUE 
C 

C 

ELSEIF(ATH. EQ. ’ N’ ) THEN 
CALL CLS 
GOTO 989 
ELSE 

CALL CLS 

WRITE(*,996) THEAT, UH 

996 FORMAT(////,' TOTAL HEAT APPLIED TO THE SURFACE IS: ’ ,F9 

+. 4, IX , A13 , ) 

GOTO 990 

ENDIF 

C 

C 

C CHOICE #2 

C 

C 

ELSEIFC SELH. EQ. ' 2 ' ) THEN 
998 WRITE(* , 997 ) UH 

997 FORMAT(////, ’ YOU HAVE SELECTED TO ENTER THE AVERAGE HEA 

+T OVER THE',/,' UPPER PCB SURFACE. ENTER THE DESIR 

+ED HEAT INPUT: (’,A13,'): ' ,2X, ) 

READ *,AHEAT 
*C 
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C MAKE ENTRY AND ALLOW FOR CORRECTION 
C 

1000 WRITE(*,999) 

999 FORMAT( / , ' IS THIS THE CORRECT ENTRY? (Y OR N): ' ,2X, ) 

READ(*,1001)ANSHA 

1001 FORMAT(Al) 

IF (ANSHA. EQ. * Y' ) THEN 
THPN=AHEAT*EL*EW/NPL 
WRITE(*, 1002) THPN.UAH 

1002 FORMATC///,' TOTAL HEAT PER NODE IS: ' ,F9. 4, IX, A6) 

WRITE(* ,4323) 

4323 FORMATC//,' <PRESS ENTER TO CONTINUE>') 

RE AD (*,43 24) DUMMY 

4324 FORMAT(Al) 

C 

C FILL HEAT MATRIX WITH DESIRED VALUES 

C 

C 

DO 1003 1=1 ,NDEEP 
DO 1004 J=1 ,NWIDE 
HEAT( I , J)=THPN 
1004 CONTINUE 

1003 CONTINUE 
C 

C 

ELSEIF( ANSHA. EQ. ' N * ) THEN 
CALL CLS 
GOTO 998 
ELSE 

CALL CLS 

WRITE(* , 1006) AHEAT.UH 

1006 FORMATC////, ' AVERAGE HEAT OVER PCB SURFACE IS: ',F9.4,1 

+X, A13 , ) 

GOTO 1000 

ENDIF 

C 

C CHOICE it 3 

C 

C 

ELSEIF( SELH. EQ. ' 3 ' ) THEN 

1200 WRITE(*, 1201) 

1201 FORMATC////,' YOU HAVE SELECTED TO ENTER THE HEAT NODALLY’ 

+,//,' ENTER THE TOTAL NUMBER OF NODES DESIGNATED FOR HEAT INP 

+UT: ' ,2X, ) 

READ * ,TOTNOD 
C 
C 

C THIS IS DONE NODE BY NODE. GET NUMBER OF ENTRIES AND THEN LOOP UNTIL 
C ALL ENTRIES HAVE BEEN MADE. 

C 

C 

C TELL USER MAXIMUM NUMBER OF ENTRIES POSSIBLE 
C 

IF (TOTNOD. GT. NPL) THEN 
WRITE(* , 1202) NPL 

1202 FORMATC///,' THE MAXIMUM ENTRY IS: ',14) 
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WRITE(*,2239) 

2239 FORMAT( // / , ' PLEASE PRESS <ENTER> TO CONTINUE ') 

READ(*, 2240) DUMMY 

2240 FORMAT(Al) 

CALL CLS 
GOTO 1200 
ENDIF 

C 

C MAKE ENTRIES 

C 

C 



DO 1203 1=1 , TOTNOD 
NC=I 



1204 

1205 
5345 



CALL CLS 

WRITE(*,1205) NC, TOTNOD 

FORMAT(///,' THIS IS NUMBER ',13,' OF ',13,' ENTRIES') 

WRITE(*,5345) 

FORMAT(//,’ ENTER THE NODE NUMBER FOR HEAT INPUT: ',2 



READ (*,8032) NN 
8032 FORMAT( 14) 

UMMY( I)=NN 

IF (NN. EQ. O.OR. NN. GT. NPL) THEN 
GOTO 1204 
ENDIF 

WRITE(*, 1206) UH 

1206 FORMAT( / , ' ENTER THE HEAT INPUT (',A6,'): ',2X, ) 

RFAU * NHFAT 

HEAT( IH(NN) , JH(NN) ) = NHEAT 
1203 CONTINUE 
C 

C PROVIDE OPPORTUNITY TO MAKE CORRECTIONS 

C 

C 

1301 CALL CLS 

WRITE(* , 1302) TOTNOD 

1302 FORMAT( /// , ' YOU HAVE MADE ',13,' NODAL ENTRIES: ' ,/) 

DO 11203 IM=1, TOTNOD 

WRITE( * ,11204) UMMY( IM) , HEAT( IH( UMMY( IM) ) , JH( UMMY( IM) ) ) , UH 
11204 FORMAT( IX, ' NODE NUMBER ’ , 14, ' : ' , IX, F9. 4, IX, A6) 

11203 CONTINUE 

WRITE(*, 1303) 

1303 FORMAT( / , ' DO YOU WISH TO MAKE ANY MORE ENTRIES OR CORRECTION 

+S? (Y OR N): ' ,2X, ) 

READ(*, 1304) AHN 

1304 FORMAT(Al) 

IF( AHN. EQ. 'Y' ) THEN 
CALL CLS 
GOTO 1200 

ELSEIF(AHN. EQ. 'N' ) THEN 
GOTO 1305 
ELSE 

GOTO 1301 

1305 ENDIF 

ELSEIF(SELH. EQ. '4' ) THEN 
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nnooon nnonnnnooooonononnn 



GOTO 1306 
1306 END IF 
END 

SUBROUTINE PCB S3 ( XEW , XEL , XCW , XCL , XUKE , XDELE , XUKC , SELECT , X 1 L 1 , X2L1 , 
+X2L2,X3L1 ,X3L2 ,X3L3 ,X4L1 ,X4L2 ,X4L3 ,X4L4 ,XWIDE ,XDEEP,XPL,XIH,XJH,X 
+EAT , XOEF , XIBT , XPRT , XWRT , XRT , XLT , XFT , XBT , XNITS ) 

C******************* ************************** ************************** 



TITLE: 

SUBROUTINE: 

AUTHOR: 

DATE: 

COMPILER: 

LINKER: 



MODEL BUILDER 
PCBS3 

LT STEVE GLASER 
09 JUL 1991 

MICROSOFT VERSION 4. 01 
MICROSOFT VERSION 3. 55 



THIS SUBROUTINE IS THE "HEART" OF THE MODEL BUILDER. THE 
NODE COEFFICIENTS ARE CALCULATED IN THIS SUBROUTINE AND 
PLACED IN MATRICES. PCBS3 IN TURN ALSO CALLS OTHER 
SUBROUTINES: SI, S2, S3, AND S4. THESE SUBROUTINES CALLED 
BY PCBS3 GENERATE THE OUTPUT DATA FILE FOR THE THERMAL 
ANALYZER. 



DEFINE REAL VARIABLES 



REAL XEL, XEW, XDELE, XUKE, XCL, XCW, X1L1 ,X2L1 ,X2L2 ,X3L1 ,X3L2 ,X3L3 ,X4L1 
+,X4L2,X4L3,X4L4,XUKC,EZB 

REAL EYLR,EFXB ,EYY,EZC11 ,EZC21 ,EZC22,EZC31 ,EZC32,EZC33,EZC41 ,EZC42 
+.EZC43 ,EZC44,EXX,DELX,DELY 

REAL CYLR1 1 , CYLR2 1 , CYLR22 , CYLR3 1 , CYLR32 , CYLR33 , CYLR4 1 , CYLR42 , CYLR4 
+3 ,CYLR44 

REAL CFXB 1 1 , CFXB21 , CFXB22 , CFXB 3 1 , CFXB 3 2 , CFXB 33 , CFXB4 1 , CFXB42 , CFXB4 
+3 ,CFXB44 

REAL CYY11 , CYY2 1 , CYY22 , CYY3 1 , CYY32 , CYY33 , CYY4 1 , CYY42 , CYY43 , CYY44 
REAL CXX1 1 , CXX2 1 , CXX22 , CXX3 1 , CXX32 , CXX33 , CXX41 , CXX42 , CXX43 , CXX44 
REAL CZE11 ,CZE21 ,CZE22 ,CZE31 ,CZE32 ,CZE33 ,CZE41 ,CZE42 ,CZE43 ,CZE44 
REAL PC11 ,PC21 ,PC22,PC31 ,PC32 ,PC33 ,PC41 ,PC42 ,PC43 ,PC44 
REAL XEAT( 100,100) ,XOEF( 1000,12) 



DEFINE CHARACTER VARIABLES 



CHARACTER* 1 SELECT, ANS, XNITS 
DEFINE INTEGER VARIABLES 
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c 

c 

INTEGER XWIDE , XDEEP , XPL , NUMC , I , N , IB , USEL 

INTEGER XIH( 1000), XJH( 1000) 

C 

c*********************************************************************** 
c***********************determine COEFFICIENTS ************************** 

Qirk'tc kkkkkkick^kkkkkkkkkkkkkkickititki'kkkk-k-kk'k'k'kic’k-k'k'k'k'k-k'kk'kk-kk'kiT'k-k'k'k-k'k-kk-kkk-k'k 

c 

c 

C DETERMINE INCREMENTAL MEASUREMENTS IN THE X AND Y DIRECTIONS 

C 

C 

C XNITS IS A CODE USED IN THE OUTPUT DATA FILE. 1 MEANS ENGLISH UNITS, 
C AND 2 MEANS SI UNITS 
C 

c 

IF (XNITS. EQ. *E* ) THEN 
USEL=1 
ELSE 

USEL=2 

ENDIF 

C 

C DETERMINE INCREMENTS IN THE ARRAY 
C 

DELX = XEL/XDEEP 
DELY = XEW/XWIDE 
1111 CONTINUE 
C 

c 

C COPPER COVERAGE, ONE COPPER LAYER CASE 

C 

C 

IF( SELECT. EQ. ’ 1 ' ) THEN 
931 CALL CLS 
C 

C ENTER PERCENT COPPER COVERAGE 

C 

C 

WRITE(* , 3941) 

3941 FORMAT( ////,' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR 

+ LAYER 1: ' ,2X, ) 

READ *,PC11 

IF (PC11. GT. 100) THEN 

GOTO 931 

ENDIF 

PC11=PC11/100 

c 

c 

C COPPER COVERAGE, TWO COPPER LAYERS CASE 

C 

C 

C 

C 
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ELSEIF( SELECT. EQ. '2' ) THEN 

932 CALL CLS 
WRITE(* , 3942) 

3942 FORMATC////, ' ENTER THE PERCENT COPPER COVERAGE ( i. e. 67) FOR 

+ LAYER 1: ' ,2X, ) 

READ *,PC21 
IF (PC21.GT. 100) THEN 
GOTO 932 
END IF 

PC21=PC21/100 

933 WRITE(*,3943) 

3943 FORMATC/,' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR LA 

+YER 2: ' ,2X, ) 

READ *,PC22 
IF (PC22. GT. 100) THEN 
GOTO 933 
END IF 

PC22=PC22/ 100 
C 
C 

C COPPER COVERAGE, THREE COPPER LAYERS CASE 

C 

C 

ELSEIF( SELECT. EQ. ' 3 ' ) THEN 
942 CALL CLS 

WRITE(*,4942) 

4942 FORMATC////,' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR 

+ LAYER 1: ' ,2X, ) 

READ *,PC31 
IF (PC31.GT. 100) THEN 
GOTO 942 
ENDIF 

PC31=PC31/100 
433 WRITE(*,4943) 

4943 FORMATC/,' ENTER THE PERCENT COPPER COVERAGE Ci. e. 67) FOR LA 

+YER 2: ' ,2X, ) 

READ * , PC32 

IF CPC32. GT. 100) THEN 

GOTO 433 

ENDIF 

PC32=PC32/100 
633 WRITEC*,8943) 

8943 FORMATC/,' ENTER THE PERCENT COPPER COVERAGE Ci- e. 67) FOR LA 

+YER 3: ' ,2X, ) 

READ *,PC33 
IF CPC33.GT. 100) THEN 
GOTO 633 
ENDIF 

PC33=PC33/100 

C 

C 

C COPPER COVERAGE, FOUR COPPER LAYERS CASE 

C 

C 

ELSEIFC SELECT. EQ. '4') THEN 
242 CALL CLS 
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WRITE (*, 2942) 

2942 FORMAT( ////,' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR 

+ LAYER 1: ' ,2X, ) 

READ *,PC41 

IF (PC41. GT. 100) THEN 

GOTO 242 

ENDIF 

PC4l=PC41/100 
233 WRITE(*,2943) 

2943 FORMAT( / , ' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR LA 

+YER 2: ' ,2X, ) 

READ * ,PC42 

IF (PC42. GT. 100) THEN 

GOTO 233 

ENDIF 

PC42=PC42/ 100 
133 WRITE(*, 1943) 

1943 FORMAT(/ , 1 ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR LA 

+YER 3: ' ,2X, ) 

RFAT) * PC43 

IF (PC43. GT. 100) THEN 

GOTO 133 

ENDIF 

PC43=PC43/100 
33 WRITE(*,43) 

43 FORMAT(/,' ENTER THE PERCENT COPPER COVERAGE (i. e. 67) FOR LA 

+YER 4: ' ,2X, ) 

READ *,PC44 

IF (PC44. GT. 100) THEN 

GOTO 33 

ENDIF 

PC44=PC44/100 
ENDIF 
CALL CLS 
C 
C 
C 

C VERIFY THAT THE SELECTED ENTRIES ARE CORRECT 
C 
C 
C 

IF( SELECT. EQ. ' 1 ’ ) THEN 
305 WRITE(* , 300) PC11*100 

300 FORMAT( /////' YOU HAVE SELECTED ' ,F6. 2 , '% COVERAGE. ' 

+ IS THIS YOUR DESIRED ENTRY? (Y OR N): ' ,1X, ) 

READ (* , 301)ANS 

301 FORMAT(Al) 

IF(ANS. EQ. ' Y’ ) THEN 
GOTO 333 

ELSEIF( ANS. EQ. 'N') THEN 

GOTO 1111 

ELSE 

CALL CLS 
GOTO 305 
ENDIF 
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ELSE IF( SELECT. EQ. ' 2' ) THEN 

306 WRITE(*,307) PC21*100 , PC22*100 

307 F0RMAT( /////' YOU HAVE SELECTED: ' ,F6. 2 , ' % COVERAGE, LAYER l' 

+ ,/,25X,F6. 2, '% COVERAGE, LAYER 2 IS THIS YOUR DESIRED 
+ ENTRY? (Y OR N): ' ,1X, ) 

READ (*,308)ANS 

308 FORMAT(Al) 

IF(ANS. EQ. 'Y') THEN 
GOTO 333 

ELSEIF( ANS. EQ. ’ N ’ ) THEN 

GOTO 1111 

ELSE 

CALL CLS 
GOTO 306 
ENDIF 

ELSEIF( SELECT. EQ. '3') THEN 

310 WRITE(*,311) PC31*100,PC32*100,PC33*100 

311 FORMAT( /////' YOU HAVE SELECTED: ' ,F6. 2 ,'% COVERAGE, LAYER l' 

+ ,/,25X,F6. 2 , '% COVERAGE, LAYER 2' ,/,25X,F6. 2 , '% COVERAGE, LAYER 3. 
+',//,' IS THIS YOUR DESIRED ENTRY? (Y OR N): ',1X, ) 

READ (*, 312)ANS 

312 FORMAT(Al) 

IF(ANS. EQ. 'Y' ) THEN 
GOTO 333 

ELSEIF(ANS.EQ. 'N') THEN 

GOTO 1111 

ELSE 

CALL CLS 
GOTO 310 
ENDIF 

ELSEIF( SELECT. EQ. ' 4 ' ) THEN 

320 WRITE(*,321) PC41*100,PC42*100,PC43*100,PC44*100 

321 FORMAT( /////' YOU HAVE SELECTED: ' ,F6. 2 , '% COVERAGE, LAYER 1* 

+ ,/ ,25X,F6. 2, '% COVERAGE, LAYER 2 ' , / , 25X,F6. 2 , ' % COVERAGE, LAYER 3' 
+ ,/ ,25X,F6. 2, '% COVERAGE, LAYER 4 IS THIS YOUR DESIRED 
+ ENTRY? (Y OR N): ' ,1X, ) 

READ (*,322)ANS 

322 FORMAT(Al) 

IF( ANS. EQ. 'Y') THEN 
GOTO 333 

ELSEIF( ANS. EQ. 'N' ) THEN 

GOTO 1111 

ELSE 

CALL CLS 
GOTO 320 
ENDIF 
ENDIF 

333 CONTINUE 

C 

C 

c 

C GENERATE CONSTANTS FOR THE EPOXY LAYERS 
C 

C LEFT OF RIGHT EDGE TO OUTSIDE 

C 

C 
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EYLR = 2*XUKE*DELX*XDELE/DELY 



FRONT OR BACK TO OUTSIDE 
EFXB = 2*XUKE*DELY*XDELE/DELX 

INNER MATRIX MOVEMENT IN THE Y DIRECTION 
EYY = XUKE*DELX*XDELE/DELY 

INNER MATRIX MOVEMENT IN THE X DIRECTION 
EXX = XUKE*XDELE*DELY/DELX 

EPOXY TO COPPER 
ONE COPPER LAYER CASE 

IF( SELECT. EQ. ' 1 ’ ) THEN 

EZC11= 2*DELX*DELY/( (XDELE/XUKE)+(X1L1*PC11/XUKC) ) 
TWO COPPER LAYERS CASE 

ELSEIF( SELECT. EQ. '2' ) THEN 

EZC21 = 2*DELX*DELY/( (XDELE/XUKE)+(X2L1*PC21/XUKC) ) 
EZC22 = 2*DELX ,V DELY / ( ( XDELE/XUKE )+( X2L2*PC22/XUKC ) ) 

THREE COPPER LAYERS CASE 
ELSEIF( SELECT. EQ. '3' ) THEN 
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EZC31 = 2*DELX*DELY/((XDELE/XUKE)+(X3L1*PC31/XUKC)) 
EZC32 = 2*DELX*DELY/((XDELE/XUKE)+(X3L2*PC32/XUKC)) 
EZC33 = 2*DELX*DELY/( (XDELE/XUKE)+(X3L3*PC33/XUKC) ) 



FOUR COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. '4') THEN 

EZC41 = 2*DELX*DELY/( (XDELE/XUKE)+(X4L1*PC41/XUKC) ) 
EZC42 = 2*DELX*DELY/( (XDELE/XUKE)+(X4L2*PC42/XUKC) ) 
EZC43 = 2*DELX*DELY/( (XDELE/XUKE)+(X4L3*PC43/XUKC) ) 
EZC44 = 2*DELX*DELY/((XDELE/XUKE)+(X4L4*PC44/XUKC)) 



ELSE 

CONTINUE 

ENDIF 



EPOXY TO BOTTOM OUTER EDGE 
EZB = 2*XUKE*DELX*DELY/XDELE 



GENERATE CONSTANTS FOR THE COPPER LAYERS 



LEFT OF RIGHT EDGE TO OUTSIDE 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. ' 1 ’ ) THEN 
CYLR11= 2*XUKC*X1L1*PC11/DELY 



TWO COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. 1 2 ' ) THEN 
CYLR21= 2*XUKC*X2L1*PC21/DELY 
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CYLR22= 2*XUKC*X2L2*PC22/DELY 



THREE COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. 1 3 ' ) THEN 
CYLR31= 2*XUKC*X3L1*PC31/DELY 
CYLR32= 2*XUKC*X3L2*PC32/DELY 
CYLR33= 2*XUKC*X3L3*PC33/DELY 



FOUR COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. '4' ) THEN 
CYLR4l= 2*XUKC*X4L1*PC41/DELY 
CYLR42= 2*XUKC*X4L2*PC42/DELY 
CYLR43= 2*XUKC*X4L3*PC43/DELY 
CYLR44= 2*XUKC*X4L4*PC44/DELY 

ELSE 

CONTINUE 

ENDIF 



FRONT OR BACK TO OUTSIDE 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. ' 1' ) THEN 

CFXB11= 2*XUKC*DELY*X1L1*PC11/DELX 



TWO COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. 1 2 1 ) THEN 
CFXB21= 2*XUKC*DELY*X2L1*PC21/DELX 
CFXB22= 2*XUKC*DELY*X2L2*PC22/DELX 



THREE COPPER LAYERS CASE 



SO 



ooooooo ooooooo oooooooo ooooooo noon 



ELSEIF( SELECT. EQ. '3') THEN 
CFXB31= 2*XUKC*DE LY*X 3 L 1* PC 3 1 / DELX 
CFXB32= 2*XUKC*DELY*X3L2*PC32/DELX 
CFXB33= 2*XUKC*DELY*X3L3*PC33/DELX 



FOUR COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. ’4’ ) THEN 
CFXB41= 2*XUKC*DELY*X4L1*PC41/DELX 
CFXB42= 2*XUKC*DELY*X4L2*PC42/DELX 
CFXB43= 2*XUKC*DELY*X4L3*PC43/DELX 
CFXB44= 2*XUKC*DELY*X4L4*PC44/DELX 

ELSE 

CONTINUE 

ENDIF 

INNER MATRIX MOVEMENT IN THE Y DIRECTION 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. * 1 ’ ) THEN 
CYY11= XUKC*DELX*X 1L1*PC 1 1/DELY 



TWO COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. '2') THEN 
CYY21= XUKC*DELX*X2L1*PC2 1/DELY 
CYY22= XUKC*DELX*X2L2*PC22/DELY 



THREE COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. ' 3 ' ) THEN 
CYY31= XUKC*DELX*X3L1*PC 3 1/DELY 
CYY32= XUKC*DELX*X3L2*PC32/DELY 
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CYY33= XUKC*DELX*X3L3*PC33/DELY 



FOUR COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. 1 4 ' ) THEN 
CYY41= XUKC*DELX*X4L1*PC41/DELY 
CYY42= XUKC*DELX*X4L2*PC42/DELY 
CYY43= XUKC*DELX*X4L3*PC43/DELY 
CYY44= XUKC*DELX*X4L4*PC44/DELY 



ELSE 

CONTINUE 
END IF 



INNER MATRIX MOVEMENT IN THE X DIRECTION 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. ' 1' ) THEN 
CXX11= XUKC*X1L1*PC11*DELY/DELX 



TWO COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. ' 2 ' ) THEN 
CXX21= XUKC*X2L1*PC21*DELY/DELX 
CXX22= XUKC*X2L2*PC22*DELY/DELX 



THREE COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. '3' ) THEN 
CXX31= XUKC*X3L1*PC31*DELY/DELX 
CXX32= XUKC*X3L2*PC32*DELY/DELX 
CXX33= XUKC*X3L3*PC33*DELY/DELX 



FOUR COPPER LAYERS CASE 
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ELSEIFC SELECT. EQ. '4') THEN 
CXX4l= XUKC*X4L1*PC41*DELY/DELX 
CXX42= XUKC*X4L2*PC42*DELY/DELX 
CXX43= XUKC*X4L3*PC43*DELY/DELX 
CXX44= XUKC*X4L4*PC44*DELY/DELX 



ELSE 

CONTINUE 

ENDIF 



COPPER TO EPOXY (OR AIR) 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. 'l') THEN 
C ZE 1 1= 2*XUKC*DE LX*DE LY/ ( X 1 L 1*PC 1 1 ) 



TWO COPPER LAYERS CASE 



ELSEIF( SELECT. EQ. '2' ) THEN 

CZE21 = 2*DELX*DELY/( (XDELE/XUKE)+(X2L1*PC21/XUKC) ) 
CZE22 = 2*XUKC*DELX*DELY/(X2L2*PC22) 



THREE COPPER LAYER CASE 



ELSEIF( SELECT. EQ. ' 3 ' ) THEN 

CZE31 = 2*DELX*DELY/((XDELE/XUKE)+(X3L1*PC31/XUKC)) 
CZE32 = 2*DE LX*DELY / ( ( XDE LE / XUKE )+(X3L2*PC32/ XUKC ) ) 
CZE33 = 2*XUKC*DELX*DELY/(PC33*X3L3) 



FOUR COPPER LAYER CASE 



ELSEIFC SELECT. EQ. '4') THEN 

CZE41 = 2*DELX*DELY/((XDELE/XUKE)+(X4L1*PC41/XUKC)) 
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CZE42 = 2*DELX*DELY/((XDELE/XUKE)+(X4L2*PC42/XUKC)) 
CZE43 = 2*DELX*DELY/((XDELE/XUKE)+(X4L3*PC43/XUKC)) 
CZE44 = 2*XUKC*DELX*DELY/(PC44*X4L4) 

ELSE 

CONTINUE 

ENDIF 

CALL THE OUTPUT DATA FILE SUBROUTINES 



ONE COPPER LAYER CASE 



IF( SELECT. EQ. * 1 ' ) THEN 

CALL S 1 ( EYLR , EFXB , EYY ,EXX , EZC 1 1 , EZB , CYLR1 1 , CFXB 1 1 , CYY1 1 , CXX1 1 , CZE 1 
+ 1 , XPL , XIH , X JH , XEAT , XWIDE , XDEEP , XOEF , XIBT , XPRT , XWRT , XRT , XLT , XFT , XBT 
+,USEL) 



TWO COPPER LAYER CASE 



ELSEIF( SELECT. EQ. '2') THEN 

CALL S2(EYLR, EFXB, EYY, EXX,EZC21,EZC22, EZB, CYLR21,CYLR22,CFXB21,CFX 
+B22,CYY21 ,CYY2 2 ,CXX21 ,CXX22, CZE2 1,CZE22, XPL, XIH, XJH, XEAT, XWIDE, XDE 
+EP , XOEF , XIBT , XPRT , XWRT , XRT , XLT , XFT , XBT , USEL) 



THREE COPPER LAYER CASE 



ELSEIF( SELECT. EQ. '3') THEN 

CALL S3( EYLR, EFXB, EYY, EXX.EZC31 ,EZC32 ,EZC33 , EZB, CYLR31,CYLR32,CYLR 
+33 ,CFXB31 ,CFXB32 ,CFXB33 ,CYY31 ,CYY32 ,CYY33 ,CXX31 ,CXX32 ,CXX33 ,CZE31 , 
+CZE32,CZE33, XPL, XIH, XJH, XEAT, XWIDE, XDEEP, XOEF, XIBT, XPRT, XWRT, XRT, X 
+LT, XFT, XBT, USEL) 



FOUR COPPER LAYER CASE 



ELSEIF( SELECT. EQ. ' 4 * ) THEN 

CALL S4(EYLR, EFXB, EYY, EXX,EZC41 ,EZC42 ,EZC43 ,EZC44, EZB, CYLR41,CYLR4 
+2 , CYLR43 , CYLR44 , CFXB4 1 , CFXB42 , CFXB43 , CFXB44 , CYY4 1 , CYY42 , CYY43 , CYY4 
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+ , XWIDE , XDEEP , XOEF , XIBT , XPRT , XWRT , XRT , XLT , XFT , XBT , USEL) 

ELSE 

CONTINUE 
END IF 
END 

SUBROUTINE COPPER( EW, EL, ANSN, UK, SELECT, T1L1,T2L1,T2L2 ,T3L1 ,T3L2 ,T3 
+L3 , T4L 1 , T4L2 , T4L3 , T4L4 , SPEVAL , UKC , CL , CW ) 



C COMPILER: MICROSOFT VERSION 4. 01 

C LINKER: MICROSOFT VERSION 3. 55 

C 

C DEFINE REAL VARIABLES 

C 

C SUBROUTINE COPPER GATHERS ALL THE INFORMATION REQUIRED FOR THE COPPER 
C LAYERS. 

REAL EL,EW,UN,T1L1 ,T2L1 ,T2L2 ,T3L1 ,T3L2 ,T3L3 ,T4L1 ,T4L2 ,T4L3 ,T4L4,V0 
+L 1 , V0L2 , V0L3 , VOL4 , T 1W 1 , T2W 1 , T2W2 , T3W 1 , T3W2 , T3W3 , T4W 1 , T4W2 , T4W3 , T4W 
+4, UKC 

INTEGER NWIDE,NDEEP,NPL 
C DEFINE ALL ONE CHARACTER VARIABLES 

C 

CHARACTER* 1 ANSN , SELECT, WORL , SPEVAL , ANSE 1 , ANST 1 , ANSK 1 

C CHARACTER VARIABLES OF MORE THAN ONE POSITION 

CHARACTER LOCO*2 ,UK*10 ,UT*1 
C 

C DEFINE MATRICES 



C 

C 

C 

C 

C 



C 

c 

5458 CALL CLS 
31171 WRITE(*, 18000) 

18000 FORMAT( / // , ' ****************************************************** 
******************** * , / , 

+ ' ******************coPPER LAYER CHARACTERISTICS******************* 
********* ' / 

+ 1 ***************************************************************** 

********* ' , / , ) 

C 

C PROVIDE CORRECT UNIT ABBREVIATIONS 
C 



c 

C TITLE: 

C AUTHOR: 

C DATE: 



MODEL BUILDER 
LT STEVE GLASER 
09 JUL 1991 



C 



COPPER LAYER CHARACTERISTICS 



AREA=EW*EL 
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IF( ANSN. EQ. 'S' ) THEN 

UN=12*2. 54*12*2. 54*. 00134374*2. 54 
WRITEC* .15400) 

15400 FORMATS ALL ENTRIES ARE IN SI NOTATION.',/) 

LOCO=' cm' 

C UK=' Watts/cm/C' 

C UT='C' 

ELSEIF( ANSN. EQ. ' E ’ ) THEN 
UN=12*12*0. 00134374 
WRITEC*, 15401) 

15401 FORMATC ALL ENTRIES ARE IN ENGLISH NOTATION’,/) 

LOCO=' in' 

C UK=' Btu/hr/F ' 

C UT='F’ 

ENDIF 

C 

C 

C 

WRITEC*, 15402) EL, LOCO 

15402 FORMATC/,' THE COPPER LAYER LENGTH IS THE SAME AS THE EPOXY LAYER 
+: ' ,F9. 4, 1X,A2,2X, ) 

CL=EL 

C 

WRITEC*, 15403) EW.LOCO 

15403 FORMATC//,' THE COPPER LAYER WIDTH IS THE SAME AS THE EPOXY LAYER 
+: ’ ,F9.4,1X,A2,2X, ) 

CW=EW 

14337 WRITEC*, 16403) SELECT 

16403 FORMATC//,' YOU SELECTED ’ , A1 , ’ COPPER LAYERCS) FOR THE PCB. ',/,’ 
+ YOU NOW WILL BE ASKED TO ENTER THE COPPER LAYER THICKNESS' ,2X, ) 
99999 WRITEC*, 31112) 

31112 FORMATC//,' DO YOU WANT TO SPECIFY THICKNESS BY LENGTH OR WEIGHT? 
+ ’ ,/,' ENTER L IF YOU WANT TO ENTER LENGTH, W IF YOU WANT TO ENTER 
+ WEIGHT: ' ,2X, ) 

READC*, 31113) WORL 

31113 FORMATC Al) 

IF( WORL. EQ. ' L' . AND. SELECT. EQ. ’ 1 ' ) THEN 
CALL CLS 

WRITEC*, 31114) LOCO 

31114 FORMATC////,' ENTER THE THICKNESS FOR LAYER 1 C',A2,'): ' ,2X, ) 
READ * ,T1L1 

GOTO 31170 

ELSEIF( WORL. EQ. ' L' . AND. SELECT. EQ. ' 2' )THEN 
CALL CLS 

WRITEC*, 31115) LOCO 

31115 FORMATC////,' ENTER THE THICKNESS FOR LAYER 1 ( ' , A2 , ' ): ' ,2X, ) 
READ * ,T2L1 

WRITEC*, 31116) LOCO 

31116 FORMATC/,' ENTER THE THICKNESS FOR LAYER 2 ( ' , A2, ' ) : ',2X, ) 
READ * ,T2L2 

GOTO 31170 

ELSEIFCWORL.EQ. ' L* . AND. SELECT. EQ. '3')THEN 
CALL CLS 

WRITEC*, 31117) LOCO 

31117 FORMATC////,' ENTER THE THICKNESS FOR LAYER 1 ( ' ,A2, ' ): ',2X, ) 
READ * ,T3L1 
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WRITE(*, 31118) LOCO 

31118 FORMAT( / , ' ENTER THE THICKNESS FOR LAYER 2 ( ’ , A2 , ' ' ,2X, ) 
READ *,T3L2 

WRITE(*, 31119) LOCO 

31119 FORMAT( / , ' ENTER THE THICKNESS FOR LAYER 3 (’ ,A2 ,’) : ',2X, ) 
READ *,T3L3 

GOTO 31170 

ELSEIF(WORL. EQ. ' L' . AND. SELECT. EQ. '4')THEN 
CALL CLS 

WRITE(*, 31120) LOCO 

31120 FORMAT( //// , ' ENTER THE THICKNESS FOR LAYER 1 ( ’ , A2 , ’ ) : \2X, ) 
READ *,T4L1 

WRITE(*,31121) LOCO 

31121 FORMAT( / , ' ENTER THE THICKNESS FOR LAYER 2 ( ’ , A2, ' ) : ’ ,2X, ) 
READ * ,T4L2 

WRITE(*,31122) LOCO 

31122 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 3 ( ’ , A2 , ’ ) : ' ,2X, ) 
READ *,T4L3 

WRITE(*, 31123) LOCO 

31123 FORMAT(/,' ENTER THE THICKNESS FOR LAYER 4 (',A2,'): ',2X, ) 
READ *,T4L4 

GOTO 31170 

ELSEIF( WORL. EQ. ' W' . AND. SELECT. EQ. ' 1 ’ ) THEN 
CALL CLS 
WRITE(*, 31142) 

31142 FORMAT( ////,’ ENTER THE WEIGHT FOR LAYER 1 (oz): ’,2X, ) 

IF(UW. EQ. 'gm' ) THEN 

READ * ,T1W1 
T1W1=T1W1*62. 5 
ELSE 

READ *,T1W1 
ENDIF 

V0L1=UN*T1W1 
T1L1=V0L1 /AREA 
GOTO 31170 

ELSEIF(WORL. EQ. 'W'.AND. SELECT. EQ. 1 2' )THEN 
CALL CLS 
WRITE(*, 31143) 

31143 FORMAT( //// , ' ENTER THE WEIGHT FOR LAYER 1 (oz): ' ,2X, ) 

IF(UW. EQ. 'gm' ) THEN 

READ * ,T2W1 
T2W1=T2W1*62. 5 
ELSE 

READ *,T2W1 
ENDIF 

V0L1=UN*T2W1 
T2L1=V0L1 /AREA 
WRITE(*, 31144) 

31144 FORMAT(/,' ENTER THE WEIGHT FOR LAYER 2 (oz): ',2X, ) 

IF(UW. EQ. 'gra' ) THEN 

READ * ,T2W2 
T2W2=T2W2*62. 5 
ELSE 

READ * ,T2W2 
ENDIF 

V0L2=UN*T2W2 
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T2L2=V0L2 / AREA 
GOTO 31170 

ELSEIF(WORL.EQ. 'W'.AND. SELECT. EQ. ' 3 1 )THEN 
CALL CLS 
WRITE(*, 31145) 

31145 FORMAT(////,’ ENTER THE WEIGHT FOR LAYER 1 (oz): ' ,2X, ) 
IF(UW. EQ. 'gm' ) THEN 

READ * ,T3W1 
T3W1=T3W1*62. 5 
ELSE 

READ * ,T3W1 
END IF 

V0L1=UN*T3W1 
T3L1=V0L1/AREA 
WRITE(*, 31146) 

31146 FORMAT ( / , ' ENTER THE WEIGHT FOR LAYER 2 (oz): ' ,2X, ) 
IF(UW. EQ. 'gm' ) THEN 

READ * ,T3W2 
T3W2=T3W2*62. 5 
ELSE 

READ *,T3W2 
ENDIF 

VOL2=UN*T3W2 
T3L2=VOL2 / AREA 
WRITE(*, 31147) 

31147 F0RMAT(/,' ENTER THE WEIGHT FOR LAYER 3 (oz): ' ,2X, ) 
IF(UW. EQ. 'gm' ) THEN 

READ * ,T3W3 
T3W3=T3W3*62. 5 
ELSE 

READ *,T3W3 
ENDIF 

VOL3=UN*T3W3 
T3L3=VOL3/AREA 
GOTO 31170 

ELSEIF(WORL. EQ. ' W’.AND. SELECT. EQ. '4')THEN 
CALL CLS 
WRITEC*, 31148) 

31148 FORMAT( ////,’ ENTER THE WEIGHT FOR LAYER 1 (oz): ' ,2X, ) 
IF(UW. EQ. 'gm' ) THEN 

READ *,T4W1 
T4W1=T4W1*62. 5 
ELSE 

READ *,T4W1 
ENDIF 

V0L1=UN*T4W1 
T4L1=V0L1/AREA 
WRITE(*, 31149) 

31149 FORMAT( / , ' ENTER THE WEIGHT FOR LAYER 2 (oz): ’ ,2X, ) 
IF(UW. EQ. 'gm' ) THEN 

READ * ,T4W2 
T4W2=T4W2*62. 5 
ELSE 

READ *,T4W2 
ENDIF 

V0L2=UN*T4W2 
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T4L2=V0L2/AREA 
WRITEC*, 31150) 

31150 FORMATC/,' ENTER THE WEIGHT FOR LAYER 3 (oz): ’,2X, ) 

IF(UW. EQ. ’gm’ ) THEN 

READ *,T4W3 
T4W3=T4W3*62. 5 
ELSE 

READ *,T4W3 
ENDIF 

VOL3=UN*T4W3 
T4L3=VOL3 / AREA 
WRITE(*, 31151) 

31151 FORMATC/,' ENTER THE WEIGHT FOR LAYER 4 (oz): ' ,2X, ) 

IF(UW. EQ. 'gm' ) THEN 

READ * ,T4W4 
T4W4=T4W4*62. 5 
ELSE 

READ *,T4W4 
ENDIF 

V 0L4=UN*T 4W 4 
T4L4=V0L4/AREA 
GOTO 31170 
CALL CLS 
GOTO 31170 
ELSE 
CALL CLS 
GOTO 31171 
ENDIF 

31170 IF(SPEVAL. EQ. ’ B ' ) THEN 
GOTO 41186 
ELSE 
CONTINUE 
ENDIF 

WRITE(*, 31172) UK 

31172 FORMATC/,' ENTER COPPER LAYER THERMAL CONDUCTIVITY ( ’ , A10 , ’ ) : 

+2X, ) 

READ * ,UKC 
C 

C MAKE CHANGES OR CORRECTIONS TO COPPER LAYER ENTRIES 
C 

31193 IF( SELECT. EQ. ' 1 ’ )THEN 
CALL CLS 

WRITEC*, 31173) CL, LOCO 

31173 FORMATC///,' YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
+R LAYER(S). ',///, 

+' 1.) LENGTH: ’ ,F9. 4, 1X,A2,2X, ) 

WRITEC*, 73) CW,LOCO 

73 FORMATC/,’ 2. ) WIDTH: ’ ,F9. 4, 1X,A2,2X, ) 

WRITEC*, 74) T1L1 , LOCO 

74 FORMATC/,' 3.) THICKNESS LAYER 1: ’ ,F9. 4, IX, A2 ,2X, ) 

WRITEC *,75) UKC ,UK 

75 FORMATC/,’ 4.) k: ’ ,F9. 4, 1X,A10,2X, ) 

ELSEIFC SELECT. EQ. ’2')THEN 

CALL CLS 

WRITEC*, 311) CL, LOCO 

311 FORMATC///,’ YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE COPPE 
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+R LAYER(S). ' ,1/1, 

+ ' 1.) LENGTH: ' ,F9. 4, IX, A2 ,2X, ) 

WRITEC*, 312) CW , LOCO 



312 


FORMAT( / , ' 
WRITEC*, 314) 
FORMATC / , ' 


2. ) WIDTH: 
T2L1.L0C0 




i 


' ,F9. 4, 1X,A2,2X, 


) 


314 


3. ) THICKNESS 
T2L2 , LOCO 


LAYER 


1: 1 


' ,F9. 4, 1X,A2,2X, 


) 


WRITEC*,315) 
FORMATC / , ' 


315 


THICKNESS 


LAYER 


2: ’ 


' ,F9. 4, 1X,A2,2X, 


) 




WRITEC*, 316) 
FORMATC / , ' 


UKC ,UK 










316 


4. ) k: 




1 


' ,F9. 4, 1X,A10,2X, 


) 




ELSEIFC SELECT. EQ. ' 3 ' )THEN 












CALL CLS 
WRITEC*, 317) 


CL, LOCO 










317 


FORMATC///, ' 


YOU HAVE MADE THE 


FOLLOWING ENTRIES FOR THE COPPE 




+R LAYER(S). ' 


,///, 












+' 1. ) LENGTH: 


' ,F9. 4, IX 


,A2 , 2X, ) 






WRITEC*, 318) 
FORMATC / , ' 


CW , LOCO 










318 


2. ) WIDTH: 
T3L1 , LOCO 






' ,F9. 4, 1X,A2,2X, 


) 


WRITEC*, 319) 
FORMATC / , ' 


319 


3. ) THICKNESS 
T3L2,LOCO 


LAYER 


1: 


' ,F9. 4, 1X,A2,2X, 


) 


WRITEC*, 320) 
FORMATC / , ' 


320 


THICKNESS 


LAYER 


2: 


' ,F9. 4, 1X,A2, 2X, 


) 




WRITEC *,321) 
FORMATC / , ' 


T3L3 ,LOCO 










321 


THICKNESS 


LAYER 


3: 


' ,F9. 4 , IX, A2 ,2X, 


) 




WRITEC*, 322) 
FORMATC / , ' 


UKC, UK 










322 


4. ) k: 






' ,F9. 4 , IX, A10 , 2X, 


) 




ELSEIFC SELECT. EQ. 4 )THEN 












CALL CLS 
WRITEC*, 323) 


CL, LOCO 










323 


FORMATC///, ' 


YOU HAVE MADE THE 


FOLLOWING ENTRIES FOR THE COPPE 




+R LAYER(S). ' 


,///, 












+ ' 1. ) LENGTH: 


' ,F9. 4 , IX 


,A2 , 2X, ) 






WRITEC*, 324) 
FORMATC / , ' 


CW , LOCO 










324 


2. ) WIDTH: 






' ,F9. 4, 1X,A2,2X, 


) 




WRITEC *,325) 


T4L1 ,LOCO 










325 


FORMATC / , ' 
WRITEC *,326) 
FORMATC / , ' 


3. ) THICKNESS 
T4L2.LOCO 


LAYER 


1: 


' ,F9. 4, 1X,A2,2X, 


) 


326 


THICKNESS 


LAYER 


2: 


' ,F9. 4, 1X,A2,2X, 


) 




WRITEC *,327) 
FORMATC / , ' 


T4L3 , LOCO 










327 


THICKNESS 


LAYER 


3: 


' ,F9. 4, 1X,A2,2X, 


) 




WRITEC*, 328) 
FORMATC / , ' 


T4L4,L0C0 










328 


THICKNESS 


LAYER 


4: 


' ,F9. 4 , IX, A2 , 2X, 


) 




WRITEC *,329) 
FORMATC / , ' 


UKC, UK 










329 


4. ) k: 






' ,F9. 4, 1X,A10, 2X, 


) 




ENDIF 













C 

31194 WRITEC*, 71180) 

71180 FORMAT( /// , ' DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR Y 

+ES AND N FOR NO: ' ,2X, ) 

READ(*,31181)ANSE1 
31181 FORMAT(Al) 

IF(ANSE1. EQ. * Y* ) THEN 
56342 WRITEC*, 31184) 

31184 FORMATC/,' WOULD YOU LIKE TO CHANGE THE THICKNESS? (Y OR N 
+): \2X, ) 
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READ(*,31185)ANST1 
31185 FORMAT(Al) 

PRINT * 

IF( ANST1. EQ. 'Y') THEN 
SPEVAL=’B' 

CALL CLS 
WRITEC*, 99998) 

99998 FORMATC////) 

GOTO 99999 

ELSEIF(ANST1. EQ. 'N' ) THEN 

CONTINUE 

ELSE 

CALL CLS 
WRITE( ,V , 56343) 

56343 FORMAT( ////) 

GOTO 56342 
ENDIF 

C 

41186 CALL CLS 

WRITEC*, 31187) 

31187 FORMATC////) 

PRINT *, ' THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS ' ,UK 
+C 1 1 UK 

32186 ’ WRITE(*, 31188) 

31188 FORMATC/,' WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCTIVIT 

+Y? (Y OR N): ' ,2X, ) 

READ(*,31189)ANSK1 

31189 FORMAT( Al) 

PRINT * 

IFCANSKl. EQ. 'Y' ) THEN 
WRITEC*, 31190) UK 

31190 FORMATC/,' ENTER THE THERMAL CONDUCTIVITY (' ,A10 ,' ): ' ,2X, ) 

READ *,UKC 

ELSEIFCANSK1.EQ. 'N') THEN 

GOTO 31191 

ELSE 

CALL CLS 
WRITEC*, 56344) 

56344 FORMATC////) 

GOTO 32186 
ENDIF 

31191 GOTO 31193 
C 

C 

ELSEIF( ANSE1. EQ. 'N') THEN 

GOTO 31180 

ELSE 

CALL CLS 
GOTO 31194 
ENDIF 
31180 END 

Q ********************************************************************** 

c 

C TITLE: MODEL BUILDER 

C SUBROUTINE: SI 

C DATE: 09 JUL 91 
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AUTHOR: LT STEVE GLASER 

COMPILER: MICROSOFT VERSION 4. 01 

LINKER: MICROSOFT VERSION 3.55 

CALLED FROM SUBROUTINE PCBS3. THIS SUBROUTINE GENERATES THE OUTPUT 

DATA FILE FOR THE ONE COPPER LAYER CASE 



SUBROUTINE SI ( GYLR , GFXB , GYY , GXX , GZC 1 1 , GZB , GYLR1 1 , GFXB 1 1 , GYY 1 1 , GXX 
+11, GZE 1 1 , GPL , GIH , G JH , GEAT , GWIDE , GDEEP , GOEF , GIBT , GPRT , GWRT , GRT , GLT , 
+GFT,GBT,GSEL) 

REAL GYLR , GFXB , GYY , GXX , GZC 1 1 , GZB , GYLR1 1 , GFXB 1 1 , GYY 1 1 , GXX1 1 , GZE 1 1 , A 
+CC , DAMP , CONFAC , GIBT , GPRT , GWRT , GRT , GLT , GFT , GBT 
REAL GEAT( 100,100) ,GOEF( 1000,12) 



CHARACTER* 1 SELECT, ANS 
CHARACTER DATAF*79 ,NAME*6 

INTEGER GPL , I , N , IB , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX ,TMAX , HT 
+RS ,D1 ,D2 ,D3 ,D4 ,D5 ,D6 ,D7 ,MAXIT,LOCVAR 

INTEGER GIH( 1000) ,GJH( 1000) ,GCON( 1000,12) 

CONTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
C0NFAC=0. 8 
COUNT=2*GPL 
LOCVAR=COUNT 



898 CALL CLS 

WRITE(* , 803) 

803 FORMAT(///,' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE' ,/, 

+' EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES',/, 

+' NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE',/, 

+' THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ' ,/, 

+' PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND',/, 

+' SHOULD NOT HAVE ANY SPACES. ’,///, 
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+' PLEASE ENTER THE DESIRED DATA FILE NAME: ',2X ) 





READ(* ,804) NAME 






804 


FORMAT( A6 ) 






552 


WRITE(*, 910) NAME 






910 


FORMAT( /////,’ YOU SELECTED ' 

+’,//) 


,A6 , ' FOR YOUR DATA 


FILE NAME 


911 


WRITE(*. 811) 

FORMATS IS THIS THE DESIRED 






811 


SELECTION? ENTER Y 


FOR YES AND N 




+FOR NO: ' ,2X, ) 








READ(* , 812) ANS 






812 


FORMAT(Al) 







C 

IF(ANS. EQ. 'N') THEN 
GOTO 898 
ELSE 

CONTINUE 

ENDIF 

IF(ANS. EQ. 'Y') THEN 
GOTO 897 
ELSE 

CALL CLS 
GOTO 552 
897 ENDIF 

C ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 

833 CALL CLS 

WRITE(*,805) 

805 FORMAT( //// , ' ENTER THE DESIRED TITLE TO BE PLACED ON LINE 

NUMBER ONE OF THE OUTPUT DATA FILE: ’,///, 

+' ' >2X, ) 





READ(*,806) DATAF 






806 


FORMAT( A79) 






835 


WRITE(* ,831) 
FORMAT( ////,' 






831 


DO YOU WISH TO CHANGE THE 


TITLE OF YOUR OUTPU 




+T DATA FILE?' ’ 


ENTER Y FOR YES AND N 


FOR NO: ’ ,2X, ) 




READ(*,832) ANS 




832 

r 


FORMAT(Al) 






u 


IF( ANS. EQ. ’ Y’ ) THEN 







GOTO 833 
ELSE 



CONTINUE 

ENDIF 

IF(ANS. EQ. 'N') THEN 
GOTO 834 
ELSE 

GOTO 835 
834 ENDIF 



C****** ******************^ 'k'k'k 'k'k'k'k'k'k'k'k'k'k'kjc'fr'kjt'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k'k 

C**************cOEFFICIENTS FOR EPOXY AND COPPER LAYERS***************** 
q******** *************************************************************** 



C 



DO 90 1=1 , GPL 
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N=1 

IB=GPL+I 



IF ((GIH(I).EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH(I).EQ. 1. OR. GJH(I). EQ 
+.GWIDE)) THEN 

DETERMINE COEFFICIENTS FOR TOP LAYER 
IF ( GEAT(GIH( I) ,GJH( I) ) . EQ. 0.0) THEN 
GC0N(I,N) = 6 

ELSE 

GC0N(I,N) = 7 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 



LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
LEFT EDGE 

IF (GJH(I). EQ. 1) THEN 
LEFT COEFFICIENT 
GOEF(I.N) = GYLR11 
GOEF( IB ,N) = GYLR 
N=N+1 

GCON(I.N) = 7551 
GCON( IB ,N) = 7551 

RIGHT COEFFICIENT 
GOEF(I,N) = GYY11 
G0EF( IB ,N) = GYY 
N=N+1 

GC0N(I,N) = 10 * ( 1+1) + 1 
GCON( IB ,N) = TO* ( IB+1) + 1 

RIGHT EDGE 

ELSEIF (GJH(I). EQ. GWIDE) THEN 

LEFT COEFFICIENT 

GOEF(I,N) = GYY11 
GOEF( IB ,N) = GYY 
N=N+1 

GCON(I,N) = 10*(I-1)+1 
GC0N( IB ,N) = 10*( IB- 1)+1 

RIGHT COEFFICIENT 

G0EF(I,N) = GYLR 11 
GOEF( IB ,N) = GYLR 
N=N+1 

GCON(I,N) = 7541 
GCON( IB ,N) = 7541 
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END IF 



FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
FRONT COEFFICIENT 
GOEF( I ,N)=GFXB11 
GOEF ( I B , N ) =GFXB 
N=N+1 

GCON( I ,N) =7521 
GC0N( IB ,N) = 7521 

BACK COEFFICIENT 
GOEF( I ,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 

GCON( I ,N) =10 ,v ( I+GWIDE)+1 
GCON( IB ,N) = 10*( IB+GWIDE)+1 

BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
FRONT COEFFICIENT 
GOEF( I ,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 

BACK COEFFICIENT 
GOEF( I ,N)=GFXB11 
GOEF( IB ,N)=GFXB 
N=N+1 

GCON( I ,N) = 7511 
GCON( IB ,N) = 7511 
END IF 

TOP COEFFICIENT 
G0EF(I,N)=GZE11 
GOEF( IB ,N)=GZC11 
N=N+1 

GCON(I,N) = 7511 
GCON( IB ,N) = 10*1+1 

BOTTOM COEFFICIENT 

G0EF(I,N) = GZC11 
G0EF( IB ,N) = GZB 
N=N+1 

GC0N(I,N) =10*( I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 

HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) , G JH( I ) ) 

N=N+1 
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GC0N( I ,N)=9991 
ENDIF 



C 

ELSEIF( ( GIH( I) . EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH(I).NE. 1. OR.GJH(I) 
+. NE.GWIDE)) THEN 
C 

C DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH(I),GJH(I)).EQ. 0. 0) THEN 
GC0N(I,N) = 6 

ELSE 

GC0N(I,N) = 7 
ENDIF 
C 

C CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 



LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 
G0EF(I,N) = GYY11 
GOEF(IB.N) = GYY 
N=N+1 

GCON( I ,N) = 10*(I-1)+1 
GCON( IB ,N) = 10*(IB-1)+1 

RIGHT COEFFICIENT 
GOEF( I ,N) = GYY 11 
GOEF( IB,N) = GYY 
N=N+1 

GCON(I,N) = 10 * (1+1) + 1 
GCON(IB.N) = 10* (IB+1) + 1 



FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I). EQ. 1) THEN 
C FRONT COEFFICIENT 
GOEF( I ,N)=GFXB11 
GOEF (IB, N) =GFXB 
N=N+1 

GCON( I ,N) =7521 
GCON( IB ,N) = 7521 

C BACK COEFFICIENT 

GOEF( I ,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 

GC0N(I,N) =10*( I+GWIDE)+1 
GCON( IB ,N) = 10*( IB+GWIDE)+1 

C BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

G0EF( I ,N)=GXX11 
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GOEF( IB ,N)=GXX 
N=N+1 

GCON(I,N) =10*(I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 

C BACK COEFFICIENT 
G0EF(I,N)=GFXB11 
GOEF( IB ,N)=GFXB 
N=N+1 

GCON( I ,N) = 7511 
GCON( IB ,N) = 7511 
END IF 

TOP COEFFICIENT 
GOEF( I ,N)=GZE11 
GOEF( IB ,N)=GZC1 1 
N=N+1 

GCON(I.N) = 7511 
GCON( IB ,N) = 10*1+1 

BOTTOM COEFFICIENT 

G0EF(I,N) = GZC11 
G0EF( IB ,N) = GZB 
N=N+1 

GC0N( I ,N) =10*( I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 

HEAT INPUT 

IF(GEAT(GIH( I) ,GJH(I)). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) , G JH( I ) ) 
N=N+1 

GCON(I,N)=9991 
ENDIF 



LEFT AND RIGHT EDGES EXCLUDING CORNERS 

ELSEIF( (GIH( I). EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH(I). EQ. 1. OR. GJH(I) 
+. EQ. GWIDE)) THEN 

DETERMINE COEFFICIENTS FOR TOP LAYER 
IF ( GEAT(GIH( I) ,GJH( I) ) . NE. 0. 0) THEN 
GC0N( I ,N) = 7 

ELSE 

GC0N( I ,N) = 6 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 



LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
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c 

C LEFT EDGE 

IF (GJH(I). EQ. 1) THEN 
C LEFT COEFFICIENT 

GOEF( I ,N) = GYLR11 
GOEF(IB.N) = GYLR 
N=N+1 

GCON( I ,N) = 7551 
GCON( IB ,N) = 7551 

RIGHT COEFFICIENT 
GOEF( I ,N) = GYY11 
GOEF(IB,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * (1+1) + 1 
GCON(IB.N) = 10* ( IB+1) + 1 

RIGHT EDGE 

ELSEIF (GJH(I). EQ. GVIDE) THEN 

LEFT COEFFICIENT 

G0EF(I,N) = GYY 11 
GOEF(IB.N) = GYY 
N=N+1 

GCON( I ,N) = 10*(I-1)+1 
GCON(IB,N) = 10*(IB-1)+1 

RIGHT COEFFICIENT 

GOEF(I,N) = GYLR11 
GOEF( IB ,N) = GYLR 
N=N+1 

GCON( I ,N) = 7541 
GCON(IB,N) = 7541 
END IF 



FRONT COEFFICIENT 
GOEF( I ,N)=GXX11 
G0EF( IB ,N)=GXX 
N=N+1 

GC0N( I ,N) =10*( I -GWIDE )+l 
GCON( IB ,N) = 10*(IB-GWIDE)+1 

C BACK COEFFICIENT 
GOEF( I ,N)=GXX11 
G0EF( IB ,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I+GWIDE)+1 
GCON( IB ,N) = 10*(IB+GWIDE)+1 



C 

C TOP COEFFICIENT 
G0EF( I ,N)=GZE11 
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GOEF( IB ,N)=GZC11 
N=N+1 

GCON( I ,N) =7511 
GCON( IB ,N) = 10*1+1 

BOTTOM COEFFICIENT 

G0EF(I,N) = GZC11 
G0EF( IB,N) = GZB 
N=N+1 

GC0N( I ,N) =10*( I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 

HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) ,GJH( I ) ) 
N=N+1 

GC0N( I ,N)=9991 
ENDIF 



****DETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE******* 



ELSEIF( ( GIH( I ) . NE. 1. OR. GIH( I ) . NE. GDEEP) . AND. (GJH( I ) . NE. 1. OR. GJH( I ) 
+. NE.GWIDE)) THEN 

DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I) ,GJH(I)).NE. 0. 0) THEN 
GC0N( I ,N) = 7 

ELSE 

GC0N(I,N) = 6 
ENDIF 



CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 



LEFT COEFFICIENT 
GOEF( I ,N) = GYY11 
GOEF( IB ,N) = GYY 
N=N+1 

GC0N( I ,N) = 10*(I-1)+1 
GC0N( IB ,N) = 10*(IB-1)+1 



RIGHT COEFFICIENT 
G0EF( I ,N) = GYY 11 
G0EF(IB,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* ( IB+1) + 1 
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c 

C FRONT COEFFICIENT 
GOEF(I,N)=GXXll 
G0EF( IB ,N)=GXX 
N=N+1 

GCON(I,N) =10*(I-GWIDE)+1 
GCON(IB.N) = 10*( IB -GWIDE)+1 

C BACK COEFFICIENT 

GOEF( I ,N)=GXX11 
GOEF( IB ,N)=GXX 
N=N+1 

GCON(I,N) =10*( I+GWIDE)+1 
GCON(IB.N) = 10*(IB+GWIDE)+1 

C 

C TOP COEFFICIENT 

GOEF( I ,N)=GZE11 
GOEF( IB ,N)=GZC11 
N=N+1 

GCON(I,N) = 7511 
GCON( IB ,N) = 10*1+1 
C 

C BOTTOM COEFFICIENT 

C 

GOEF(I,N) = GZC11 
GOEF( IB ,N) = GZB 
N=N+1 

GCON(I,N) =10*( I+GPL)+1 
GCON(IB.N) = 10*( IB+GPL)+1 
C 

C HEAT INPUT 

IF( GEAT( GIH( I ) , GJH( I ) ) . NE. 0. 0 ) THEN 
GOEF(I,N) = GEAT( GIH( I ) , GJH( I ) ) 

N=N+1 

GCON(I,N)=9991 

ENDIF 

ENDIF 

90 CONTINUE 

C 

C GENERATE DATA FILE 

C 

OPEN ( 3, F I LE=NAME,FORM=' FORMATTED' ,ACCESS=' DIRECT' , RECL=108 , STATUS 
+= ' NEW ' ) 

WRITE( 3 , 909) DATAF 
909 FORMAT( 1X,A79) 

WRITE ( 3 , 908 ) COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908 FORMAT( 2X, 9( 13 ,5X) ) 

WRITE( 3 ,907) ZER, ZER, ZER 
907 FORMAT( 2X, 3( 13 ,5X) ) 

WRITE( 3,9081) NMAX , TMAX , HTRS ,D1 ,D2 ,D3 ,D4,D5 ,D6 ,D7 
9081 FORMAT( 2X,9(I3,5X)) 

WRITE (3, 905) ACC .DAMP, MAXIT.CONFAC.GIBT 

905 FORMAT( 1X,2(F9. 7 , IX) , 14X, 12 , 1X.F9. 7.1X.F9.5) 

WRITE( 3,906) GPRT , GFT , GBT , GRT , GLT , GWRT 

906 FORMAT( 1X,6(F12. 4, IX)) 
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DO 112 I=l,LOCVAR 

WRITE( 3,9100) GCON( 1,1) ,GCON( 1,2) ,GCON( 1,3) ,GCON( 1,4) ,GCON( I ,5 ) ,GC 
+ON( 1,6), GCON( 1,7) ,GCON( 1,8) 

9100 FORMAT( I4,3X,7(I5,7X)) 

WRITE(3,9110) GOEF( 1,1) ,GOEF( 1,2) ,GOEF( 1,3) ,GOEF( 1,4) ,GOEF( I , 5) ,GO 
+EF( 1,6), GOEF( 1,7) 

9110 FORMAT(7(F9. 3,3X)) 

112 CONTINUE 
CLOSE (3) 

CALL CLS 

WRITE(*,999) NAME 

999 FORMAT( ////,' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 

+NAMED ' ,A6, /////,' <PRESS ENTER TO CONTINUE^) 

READ(*,5912) ANS 
5912 FORMAT(Al) 



END 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



TITLE: 

SUBROUTINE: 

DATE: 

AUTHOR: 

COMPILER: 

LINKER: 



MODEL BUILDER 
S2 

09 JUL 91 
LT STEVE GLASER 
MICROSOFT VERSION 4.01 
MICROSOFT VERSION 3.55 



S2 IS CALLED FROM PCBS3. S2 IS THE OUTPUT DATA FILE GENERATOR FOR THE 
THE COPPER LAYER CASE. 



SUBROUTINE S2 (GYLR,GFXB,GYY,GXX,GZC21,GZC22,GZB,GYLR21,GYLR22,GFX 
+B2 1 , GFXB22 , GYY2 1 , GYY22 , GXX2 1 , GXX22 , GZE2 1 , GZE22 , GPL , GIH , GJH , GEAT, GW 
+ I DE , GDEEP , GOEF , G I BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GSEL ) 



REAL GYLR , GFXB , GYY , GXX , GZC2 1 , GZB , GYLR2 1 , GFXB21 , GYY2 1 , GXX21 , GZE2 1 , A 
+CC, DAMP, CONF AC, GIBT, GPRT, GWRT, GRT, GLT, GFT, GBT, GYLR22,GFXB 22, GYY22, 
+GXX22,GZE22,GZC22 



REAL GEAT( 100, 100), GOEF( 1000,12) 



CHARACTER* 1 SELECT, ANS 
CHARACTER DATAF*79 ,NAME*6 

INTEGER GPL , I , N , I B , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX ,TMAX , HT 
+RS ,D1 ,D2 ,D3 ,D4 ,D5 ,D6 ,D7 ,MAXIT, LOCVAR, IC , ID 

INTEGER GIH( 1000) ,GJH( 1000) ,GCON( 1000,12) 

CONTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 
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D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
CONFAC=0. 8 
C0UNT=4*GPL 
LOCVAR=COUNT 



C 

C 

898 CALL CLS 

WRITE(*,803) 

803 FORMAT(///,' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE' ,/, 

+' EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES',/, 

+' NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE',/, 

+’ THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ’,/, 

+' PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND',/, 

+’ SHOULD NOT HAVE ANY SPACES. ',///, 

+' PLEASE ENTER THE DESIRED DATA FILE NAME: ',2X ) 

READ(*,804) NAME 

804 FORMAT(A6) 

552 WRITE (* ,9 10) NAME 

910 FORMAT( /////,' YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 

+',//) 

911 WRITE(* 811) 

811 FORMAT( ' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N 

+FOR NO: ' ,2X, ) 

READ(* , 812) ANS 

812 FORMAT(Al) 

C 

IF( ANS. EQ. 'N' ) THEN 
GOTO 898 
ELSE 

CONTINUE 

ENDIF 

IF( ANS. EQ. ’ Y' ) THEN 
GOTO 897 
ELSE 

CALL CLS 
GOTO 552 
897 ENDIF 

C ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 

833 CALL CLS 

WRITE(*,805) 

805 FORMATC////,’ ENTER THE DESIRED TITLE TO BE PLACED ON LINE 

+ ' ,/,' NUMBER ONE OF THE OUTPUT DATA FILE: ' ,///, 

+' ' , 2X ) 

READ(*,806)’DATAF 

806 FORMAT( A79) 
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835 WRITE(*,831) 

831 FORMATC////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 

+T DATA FILE?',/,' ENTER Y FOR YES AND N FOR NO: ',2X, ) 

READ(* , 832) ANS 

832 FORMAT( Al) 

C 

IF( ANS. EQ. ’Y’ ) THEN 
GOTO 833 
ELSE 

CONTINUE 

ENDIF 

IF(ANS.EQ. V) THEN 
GOTO 834 
ELSE 

GOTO 835 
834 ENDIF 



C**************C0EFFICIENTS FOR EPOXY AND COPPER LAYERS***************** 

Q***************************** ********************************* ********* 



c 

DO 90 1=1, GPL 
N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

C 

C*********************************CORNERS******************************** 

c 

IF ((GIH(I).EQ. 1. OR. GIH( I ) . EQ. GDEEP) . AND. (GJH(I).EQ. 1. OR. GJH( I ) . EQ 
+. GWIDE)) THEN 
C 

C DETERMINE COEFFICIENTS FOR TOP COPPER 
IF ( GEAT( GIH( I ) ,GJH( I ) ) . EQ. 0.0) THEN 
GCON( I ,N) = 6 

ELSE 

GCON( I ,N) = 7 
ENDIF 
C 

C CONNECTIONS FOR EPOXY LAYER 

GCON(IB.N) = 6 
GCON( ID ,N) = 6 
GCON( IC,N) = 6 
C 
C 

C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 

C 

C LEFT EDGE 

IF (GJH(I).EQ. 1) THEN 
C LEFT COEFFICIENT 

GOEF( I ,N) = GYLR21 
GOEF(IC,N) = GYLR22 
GOEF( IB ,N) = GYLR 
GOEF(ID.N) = GYLR 
N=N+1 
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GCON(I,N) = 7551 
GCON(IB.N) = 7551 
GC0N( IC ,N) = 7551 
GCON( ID ,N) = 7551 

RIGHT COEFFICIENT 
GOEF( I ,N) = GYY21 
GOEF(IC.N) = GYY22 
G0EF( IB ,N) = GYY 
GOEF(ID.N) = GYY 
N=N+1 

GCON(I.N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* (IB+1) + 1 
GC0N(IC,N) = 10 * (IC+1) + 1 

GC0N( ID,N) = 10 * (ID+1) + 1 

RIGHT EDGE 

ELSEIF (GJH(I).EQ. GWIDE) THEN 

LEFT COEFFICIENT 

G0EF(I,N) = GYY21 
G0EF( IC ,N) = GYY22 

GOEF(IB.N) = GYY 

GOEF(ID.N) = GYY 

N=N+1 

GC0N( I ,N) = 10*(I-1)+1 
GCON(IB.N) = 10*( IB-1)+1 

GCON(IC.N) = 10*(IC-1)+1 

GC0N(ID,N) = 10*( ID- 1)+1 



RIGHT COEFFICIENT 

GOEF( I ,N) = GYLR21 
G0EF( IC ,N) = GYLR22 
G0EF( IB ,N) = GYLR 
GOEF(ID.N) = GYLR 
N=N+1 

GC0N( I ,N) = 7541 
GCON(IB.N) = 7541 
GC0N( IC ,N) = 7541 
GCON(ID.N) = 7541 
ENDIF 

FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
C FRONT COEFFICIENT 
G0EF( I ,N)=GFXB21 
G0EF( IC ,N)=GFXB22 
GOEF( IB ,N)=GFXB 
GOEF(ID,N)=GFXB 
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N=N+1 

GCONC I ,N) =7521 
GCON( IB ,N) = 7521 
GCONC IC,N) = 7521 
GCON( ID, N) = 7521 

C BACK COEFFICIENT 
GOEFC I ,N)=GXX21 
GOEFC IC ,N)=GXX22 
GOEFC IB ,N)=GXX 
GOEFC ID, N)=GXX 
N=N+1 

GCONC I, N) =10*C I+GWIDE)+1 
GCON( IB ,N) = 10*C IB+GWIDE )+l 
GCONC IC,N) = 10*(IC+GWIDE)+1 
GCONC ID, N) = 10*(ID+GWIDE)+1 



C BACK EDGE 

ELSEIF C GIH( I ) . EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEFC I ,N)=GXX21 
GOEFC IC,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEFC ID, N)=GXX 
N=N+1 

GCONC I, N) =10*CI-GWIDE)+1 
GCONC IB, N) = 10*(IB-GWIDE)+1 
GCONC IC,N) = 10*C IC-GWIDE)+1 
GCONC ID, N) = 10*(ID-GWIDE)+1 

C BACK COEFFICIENT 

GOEFC I ,N)=GFXB21 
GOEFC IC,N)=GFXB22 
GOEFC IB ,N)=GFXB 
GOEF C I D , N ) =GFXB 
N=N+1 

GCONC I, N) = 7511 
GCON( IB ,N) = 7511 
GCONC IC ,N) = 7511 
GCON( ID,N) = 7511 

ENDIF 

C 

C TOP COEFFICIENT 
GOEFC I ,N)=GZE21 
GOEFC IC,N)=GZE22 
GOEF( IB, N)=GZC21 
GOEFC ID, N)=GZC22 
N=N+1 

GCON( I, N) = 7511 
GCON( IC ,N) = 7511 
GCONC IB, N) = 10*1+1 
GCONC ID, N) = 10*IC+1 
C 
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BOTTOM COEFFICIENT 

GOEF(I,N) = GZC21 
GOEF(IC,N) = GZC22 
GOEF( IB,N) = GZE22 
GOEF(ID,N) = GZB 
N=N+1 

GCON(I.N) =10*( I+GPL)+1 
GCON( IB ,N) = 10*(IB+GPL)+1 
GCON(IC,N) = 10*( IC+GPL)+1 
GCON( ID ,N) = 10*( ID+GPL)+1 



HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT(GIH( I) ,GJH( I) ) 
N=N+1 

GC0N( I ,N)=9991 
ENDIF 



ELSEIF( (GIH( I ). EQ. 1. OR. GIH(I). EQ. GDEEP) . AND. (GJH(I).NE. l.OR.GJH(I) 
+. NE.GWIDE)) THEN 

DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH( I ) , GJH( I ) ) . EQ. 0. 0) THEN 
GCON( I ,N) = 6 

ELSE 

GCON( I ,N) = 7 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GCON( IB ,N) = 6 
GCON(ID.N) = 6 
GCON( IC,N) = 6 



LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 
GOEF(I,N) = GYY21 
GOEF(IC.N) = GYY22 
GOEF(IB.N) = GYY 
GOEF(ID.N) = GYY 
N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON(IB.N) = 10*(IB-1)+1 
GCON(IC.N) = 10*(IC-1)+1 
GCON(ID.N) = 10*(ID-1)+1 



RIGHT COEFFICIENT 
GOEF( I ,N) = GYY21 
GOEF(IC.N) = GYY 2 2 
GOEF(IB.N) = GYY 
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GOEF( ID,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * (1+1) + 1 
GCON( IB ,N) = 10* (IB+1) + 1 
GC0N( IC ,N) = 10 * (IC+1) + 1 
GC0N( ID ,N) = 10 * (ID+1) + 1 



FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I). EQ. 1) THEN 
C FRONT COEFFICIENT 
G0EF(I,N)=GFXB21 
GOEF( IC,N)=GFXB22 
GOEF(IB,N)=GFXB 
GOEF ( ID , N ) =GFXB 
N=N+1 

GCON( I ,N) =7521 
GCON( IB ,N) = 7521 
GCON( IC ,N) = 7521 
GCON( ID,N) = 7521 

C BACK COEFFICIENT 

GOEF( I ,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I+GWIDE)+1 
GCON( IB ,N) = 10*( IB+GWIDE)+1 
GCON( IC ,N) = 10*( IC+GWIDE)+1 
GCON( ID ,N) = 10*(ID+GWIDE)+1 

C BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEF( I ,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 
GCON(IC,N) = 10*( IC-GWIDE)+1 
GCON(ID,N) = 10*( ID-GWIDE)+1 

C BACK COEFFICIENT 

GOEF( I , N)=GFXB2 1 
GOEF( IC ,N)=GFXB22 
GOEF( IB ,N)=GFXB 
GOEF( ID,N)=GFXB 
N=N+1 

GCON( I ,N) = 7511 
GCON( IB , N) = 7511 
GCON( IC ,N) = 7511 
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GCON(ID.N) = 7511 
ENDIF 

TOP COEFFICIENT 
GOEF( I ,N)=GZE21 
GOEF(IC,N)=GZE22 
GOEF( IB ,N)=GZC21 
GOEF( ID,N)=GZC22 
N=N+1 

GCON( I ,N) = 7511 
GCON(IC.N) = 7511 
GCON( IB,N) = 10*1+1 
GCON(ID.N) = 10*IC+1 

BOTTOM COEFFICIENT 

G0EF( I ,N) = GZC21 
G0EF( IC,N) = GZC22 
GOEF(IB.N) = GZB 
G0EF( ID,N) = GZB 
N=N+1 

GC0N( I ,N) =10*(I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 
GC0N(IC,N) = 10*(IC+GPL)+1 
GC0N( ID,N) = 10*( ID+GPL)+1 

HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) ,GJH( I ) ) 
N=N+1 

GC0N( I ,N)=999 1 
ENDIF 



LEFT AND RIGHT EDGES EXCLUDING CORNERS 

ELSEIF( (GIH( I). EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH(I). EQ. 1. OR. GJH(I) 
+. EQ. GWIDE)) THEN 

DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH( I ) ,GJH( I ) ). NE. 0. 0) THEN 
GCON( I ,N) = 7 

ELSE 

GCON( I ,N) = 6 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GCON( IB ,N) = 6 
GCON(IC.N) = 6 
GCON(ID,N) = 6 
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C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 

C 

C LEFT EDGE 

IF (GJH(I). EQ. 1) THEN 
C LEFT COEFFICIENT 

GOEF( I ,N) = GYLR21 
GOEF(IC.N) = GYLR22 
GOEF(IB.N) = GYLR 
GOEF(ID.N) = GYLR 
N=N+1 

GCON( I ,N) = 7551 
GCON( IB ,N) = 7551 
GCON(IC,N) = 7551 
GCON(ID,N) = 7551 
C 

C RIGHT COEFFICIENT 

GOEF(I,N) = GYY21 
GOEF(IC,N) = GYY22 
GOEF( IB ,N) = GYY 
GOEF( ID,N) = GYY 
N=N+1 

GCON(I,N) = 10 * (1+1) + 1 
GCON( IB ,N) = 10* ( IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 
GCON(ID,N) = 10 * ( ID+1) + 1 
C 

C RIGHT EDGE 

ELSEIF (GJH(I). EQ. GWIDE) THEN 
C 

C LEFT COEFFICIENT 

GOEF(I,N) = GYY21 
GOEF( IC,N) = GYY22 
GOEF( IB,N) = GYY 
GOEF( ID,N) = GYY 

N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON( IB ,N) = 10*( IB-1)+1 
GC0N(IC,N) = 10*(IC-1)+1 
GCON(ID,N) = 10*(ID-1)+1 
C 

C RIGHT COEFFICIENT 

C 

GOEF( I ,N) = GYLR21 
GOEF(IC,N) = GYLR22 
GOEF( IB ,N) = GYLR 
GOEF(ID,N) = GYLR 
N=N+1 

GCON( I ,N) = 7541 
GCON( IB ,N) = 7541 
GCON(IC,N) = 7541 
GCON( ID,N) = 7541 
ENDIF 
C 
C 

C FRONT COEFFICIENT 
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GOEF( I ,N)=GXX21 
G0EF( IC ,N)=GXX22 
GOEF( IB ,N)=GXX 
G0EF( ID,N)=GXX 
N=N+1 

GCON( I ,N) =10*( I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 
GCON( IC ,N) = 10*(IC-GWIDE)+1 
GCON( ID,N) = 10*(ID-GWIDE)+1 

C BACK COEFFICIENT 
G0EF( I ,N)=GXX21 
GOEF( IC ,N)=GXX22 
GOEF(IB,N)=GXX 
GOEF( ID ,N)=GXX 
N=N+1 

GC0N(I,N) =10*(I+GWIDE)+1 
GCON( IB,N) = 10*( IB+GWIDE)+1 
GC0N( IC ,N) = 10*(IC+GWIDE)+1 
GCON( ID,N) = 10*( ID+GWIDE)+1 



TOP COEFFICIENT 
GOEF( I ,N)=GZE21 
G0EF( IC ,N)=GZE22 
GOEF( IB ,N)=GZC21 
G0EF(ID,N)=GZC22 
N=N+1 

GCON( I ,N) = 7511 
GCON(IC.N) = 7511 
GCON( IB ,N) = 10*1+1 
GC0N(ID,N) = 10*IC+1 

BOTTOM COEFFICIENT 

G0EF( I ,N) = GZC21 
G0EF( IC ,N) = GZC22 
G0EF( IB ,N) = GZE22 
G0EF( ID,N) = GZB 
N=N+1 

GC0N( I ,N) =10*(I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 
GC0N( IC ,N) = 10*( IC+GPL)+1 
GC0N( ID,N) = 10*( ID+GPL)+1 

HEAT INPUT 

IF(GEAT(GIH( I) ,GJH(I)). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) , GJH( I ) ) 
N=N+1 

GC0N( I ,N)=9991 
ENDIF 



c**** DE termine COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE******* 
C 
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ELSEIF( (GIH( I ) . NE. 1. OR. GIH( I). NE. GDEEP). AND. (GJH(I).NE. 1. OR. GJH(I) 
+. NE.GWIDE)) THEN 

DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I) ,GJH( I) ). NE. 0. 0) THEN 
GCON( I ,N) = 7 

ELSE 

GC0N( I ,N) = 6 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GCON( IB ,N) = 6 
GC0N(IC,N) = 6 
GCON( ID ,N) = 6 



LEFT COEFFICIENT 
GOEF(I,N) = GYY21 
GOEF(IC,N) = GYY22 

G0EF( IB,N) = GYY 

GOEF( ID,N) = GYY 

N=N+1 

GCON( I ,N) = 10 * (1-1) + 1 
GCON( IB ,N) = 10* (IB-1) + 1 
GCON(IC,N) = 10 * (IC-1) + 1 

GCON( ID ,N) = 10 * (ID-1) + 1 



RIGHT COEFFICIENT 
G0EF( I ,N) = GYY21 
GOEF(IC,N) = GYY 2 2 
GOEF( IB ,N) = GYY 
GOEF( ID,N) = GYY 
N=N+1 

GCON(I,N) = 10 * (1+1) + 1 
GCON( IB ,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 
GCON( ID ,N) = 10 * (ID+1) + 1 



FRONT COEFFICIENT 
GOEF(I,N)=GXX21 
GOEF(IC,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GC0N(I,N) =10*( I-GWIDE)+1 
GC0N( IB ,N) = 10*(IB-GWIDE)+1 
GC0N(IC,N) = 10*(IC-GWIDE)+1 
GC0N(ID,N) = 10*( ID-GWIDE)+1 



C BACK COEFFICIENT 
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GOEF( I ,N)=GXX21 
GOEF(IC,N)=GXX22 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I+GWIDE)+1 
GCON( IB ,N) = 10*( IB+GWIDE)+1 
GCON(IC,N) = 10*(IC+GWIDE)+1 
GCON(ID.N) = 10*( ID+GWIDE)+1 



TOP COEFFICIENT 
GOEF( I ,N)=GZE21 
GOEF( IC ,N)=GZE22 
GOEF( IB ,N)=GZC21 
GOEF(ID,N)=GZC22 
N=N+1 

GCON( I ,N) =7511 
GCON( IB ,N) = 10*1+1 
GCON(IC,N) = 7511 
GCON( ID,N) = 10*IC+1 



BOTTOM COEFFICIENT 

GOEF( I ,N) = GZC21 
GOEF( IC ,N) = GZC22 
GOEF(IB,N) = GZB 
GOEF( ID ,N) = GZE22 
N=N+1 

GCON( I ,N) =10*( I+GPL)+1 
GCON(IB.N) = 10*( IB+GPL)+1 
GCON( IC ,N) = 10*(IC+GPL)+1 
GCON( ID ,N) = 10*( ID+GPL)+1 



HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I)). NE. 0.0) THEN 
GOEF( I ,N) = GEAT( GIH( I ) , GJH( I ) ) 

N=N+1 

GCON( I ,N)=999 1 
ENDIF 
END IF 
CONTINUE 

GENERATE DATA FILE 

OPEN ( 3, F I LE=NAME,FORM=' FORMATTED' ,ACCESS='DIRECT’ ,RECL=108, STATUS 
+= ' NEW ' ) 

WRITE( 3 ,909) DATAF 
909 FORMAT( 1X,A79) 

WRITE( 3 , 9 08 ) COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908 FORMAT( 2X,9( 13 ,5X) ) 

WRITE( 3 ,907 ) ZER, ZER, ZER 
907 FORMAT( 2X, 3( 13 ,5X) ) 
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WRITE( 3 , 9081 ) NMAX , TMAX , HTRS ,D1 ,D2 ,D3 ,D4 ,D5 ,D6 ,D7 
9081 F0RMAT( 2X, 9( 13 , 5X) ) 

WRITE( 3,905) ACC , DAMP , MAXIT , CONFAC , GIBT 

905 FORMAT( IX, 2(F9. 7 , IX) , 14X, 12 , 1X,F9. 7,1X,F9.5) 

WRITE (3, 906) GPRT, GFT , GBT , GRT , GLT , GWRT 

906 FORMAT( 1X,6(F12. 3, IX)) 

DO 112 1=1 ,LOCVAR 

WRITE( 3,9100) GCON( 1,1), GCON( 1,2), GCON( 1,3), GCON( 1,4), GCON( I , 5 ) , GC 
+0N( 1,6) ,GCON( 1,7) ,GC0N( 1,8) 

9100 FORMAT( 14 , 3X, 7( 15 , 7X) ) 

WRITE( 3 , 9110) GOEF( 1,1) ,GOEF( 1,2) ,GOEF( 1,3) ,GOEF( 1,4) ,GOEF( I , 5) , GO 
+EF (1,6), GOEF (1,7) 

9110 FORMAT( 7(F9. 3,3X)) 

112 CONTINUE 
CLOSE (3) 

CALL CLS 

WRITE(* , 999 ) NAME 

999 FORMAT(////, ' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 

+NAMED ' ,A6, /////,' <PRESS ENTER TO CONTINUE^) 

READ(*,5912) ANS 
5912 FORMAT(Al) 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



END 



TITLE: 

SUBROUTINE: 

DATE: 

AUTHOR: 

COMPILER: 

LINKER: 



MODEL BUILDER 
S3 

09 JUL 91 
LT STEVE GLASER 
MICROSOFT VERSION 4.01 
MICROSOFT VERSION 3. 55 



SUBROUTINE S3 IS CALLED FROM SUBROUTINE PCBS3. S3 GENERATES THE 
OUTPUT DATA FILE FOR THE THREE COPPER LAYER CASE. 



SUBROUTINE S3 ( GYLR , GFXB , GYY , GXX , GZC3 1 , GZC32 , GZC33 , GZB , GYLR3 1 , GYLR 
+32 , GYLR33 , GFXB31 , GFXB32 , GFXB 3 3 , GYY3 1 , GYY32 , GYY33 , GXX3 1 , GXX32 , GXX33 
+ , GZE3 1 , GZE32 , GZE33 , GPL , GIH , G JH , GEAT , GWIDE , GDEEP , GOEF , GIBT , GPRT , GWR 
+T , GRT , GLT , GFT , GBT , GSEL) 



REAL GYLR , GFXB , GYY , GXX , GZC3 1 , GZB , GYLR3 1 , GFXB 31 , GYY 3 1 , GXX3 1 , GZE3 1 , A 
+CC , DAMP , CONFAC , G I BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GYLR32 , GFXB 32 , GYY3 2 , 
+GXX32 , GZE32 , GZC32 , GYLR33 , GFXB 3 3 , GYY33 , GXX33 , GZE33 , GZC33 



REAL GEAT( 100,100) ,GOEF( 1000 ,12) 



CHARACTER* 1 SELECT, ANS 
CHARACTER DATAF*79 .NAME* 6 

INTEGER GPL , I , N , IB , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS,D1 ,D2,D3,D4,D5 ,D6,D7 ,MAXIT,LOCVAR, IC, ID, IE, IF 

INTEGER GIH( 1000) ,GJH( 1000) ,GCON( 1000 , 12) 

CONTEMP=6 

ZER=0 
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NMAX=750 

TMAX=5 0 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
C0NFAC=0. 8 
C0UNT=6*GPL 
L0CVAR=C0UNT 



C 

C 

898 CALL CLS 

WRITE(*,803) 

803 FORMATC///,' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE’ ,/, 

+' EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES’,/, 

+' NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE',/, 

+' THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ',/, 

+' PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND',/, 

+' SHOULD NOT HAVE ANY SPACES. ',///, 

+' PLEASE ENTER THE DESIRED DATA FILE NAME: ' ,2X ) 

READ(* , 804) NAME 

804 FORMAT( A6) 

552 WRITE(*,910)NAME 

910 FORMATC/////,' YOU SELECTED ',A6,' FOR YOUR DATA FILE NAME 

+',//) 

911 WRITE(* 811) 

811 FORMATC ' IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND N 

+FOR NO: ' ,2X, ) 

READ(*,812) ANS 

812 FORMAT(Al) 

C 

IF( ANS. EQ. 'N' ) THEN 
GOTO 898 
ELSE 

CONTINUE 

ENDIF 

IF( ANS. EQ. 'Y') THEN 
GOTO 897 
ELSE 

CALL CLS 
GOTO 552 
897 ENDIF 

C ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 
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833 CALL CLS 
WRITE(*,805) 

805 FORMAT( ////,' ENTER THE DESIRED TITLE TO BE PLACED ON LINE 

NUMBER ONE OF THE OUTPUT DATA FILE: ',///, 

+ T T , 2X, ) 

READ(*,806) *DATAF 

806 FORMAT( A79) 

835 WRITE(*,831) 

831 FORMAT( ////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 

+T DATA FILE?',// ENTER Y FOR YES AND N FOR NO: ' ,2X, ) 

READ(*,832) ANS 

832 FORMAT(Al) 

C 

IF(ANS. EQ. ’Y’) THEN 
GOTO 833 
ELSE 

CONTINUE 
END IF 

IF(ANS. EQ. ’ N' ) THEN 
GOTO 834 
ELSE 

GOTO 835 

834 END IF 
C 



C**************C0EFFICIENTS FOR EPOXY AND COPPER LAYERS***************** 



c 

DO 90 1=1, GPL 
N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

IE=4*GPL+I 

IF=5*GPL+I 



C 



C**7v********7v**rr*****************7TC0RNERS ******* ************************* 

c 

IF ((GIH(I).EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH( I). EQ. 1. OR. GJH( I). EQ 
+. GWIDE)) THEN 
C 

C DETERMINE COEFFICIENTS FOR TOP COPPER 
IF ( GEAT(GIH( I ) ,GJH( I ) ) . EQ. 0.0) THEN 
GCON( I ,N) = 6 

ELSE 



C 

C 



GCON( I ,N) = 7 
ENDIF 



CONNECTIONS FOR 
GCON( IB ,N) = 6 
GCON( ID,N) = 6 
GCON( IC,N) = 6 
GCON( ID,N) = 6 



EPOXY LAYER 
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GCON(IE,N) = 6 
GCON( IF ,N) = 6 



LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
LEFT EDGE 

IF (GJH(I).EQ. 1) THEN 
LEFT COEFFICIENT 
GOEF( I ,N) = GYLR31 
GOEF(IC.N) = GYLR32 
GOEF( IE ,N) = GYLR33 
G0EF( IB ,N) = GYLR 
GOEF( ID,N) = GYLR 
G0EF( IF,N) = GYLR 
N=N+1 

GCON(I.N) = 7551 
GCON( IB ,N) = 7551 
GC0N( IC ,N) = 7551 
GCON( ID ,N) = 7551 
GCON( IE,N) = 7551 
GCON( IF ,N) = 7551 



RIGHT COEFFICIENT 
GOEF( I ,N) = GYY31 
GOEF( IC ,N) = GYY32 
G0EF( IE ,N) = GYY33 
GOEF( IB ,N) = GYY 
GOEF(ID,N) = GYY 

G0EF( IF,N) = GYY 

N=N+1 

GCON(I,N) = 10 * (1+1) + 1 
GCON( IB ,N) = 10* ( IB+1) + 1 
GCON( IC ,N) = 10 * (IC+1) + 1 

GCON( ID ,N) = 10 * (ID+1) + 1 

GCON( IE ,N) = 10 * ( IE+1) + 1 

GCON(IF.N) = 10 * (IF+1) + 1 



RIGHT EDGE 

ELSEIF (GJH(I).EQ. GWIDE) THEN 

LEFT COEFFICIENT 

GOEF( I ,N) = GYY31 
GOEF(IC.N) = GYY 3 2 
G0EF( IE ,N) = GYY33 
GOEF(IB.N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 

GCON( I ,N) = 10*(I-1)+1 
GCON( IB ,N) = 10*(IB-1)+1 
GCON( IC ,N) = 10*(IC-1)+1 
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GCON( ID,N) = 10*(ID-1)+1 
GCON(IE.N) = 10*( IE-1)+1 
GCON(IF.N) = 10*( IF- 1)+1 

RIGHT COEFFICIENT 

GOEF(I,N) = GYLR31 
GOEF(IC,N) = GYLR32 
GOEF(IE.N) = GYLR33 
GOEF(IB.N) = GYLR 
GOEF(ID,N) = GYLR 
GOEF(IF.N) = GYLR 
N=N+1 

GCON(I,N) = 7541 
GCON( IB ,N) = 7541 
GCON(IC.N) = 7541 
GCON(ID.N) = 7541 
GCON(IE.N) = 7541 
GCON(IF.N) = 7541 

ENDIF 

FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
C FRONT COEFFICIENT 
GOEF( I ,N)=GFXB31 
GOEF(IC,N)=GFXB32 
GOEF( IE ,N)=GFXB33 
GOEF ( I B , N ) =GFXB 
GOEF( ID,N)=GFXB 
GOEF(EF,N)=GFXB 
N=N+1 

GCON(I,N) = 7521 
GCON( IB ,N) = 7521 
GCON(IC,N) = 7521 
GCON( ID ,N) = 7521 
GCON( IE ,N) = 7521 
GCON( IF,N) = 7521 
C BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
GOEF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
N=N+1 

GCON( I ,N) =10*( I+GWIDE)+1 
GCON( IB,N) = 10*( IB+GWIDE)+1 
GCON(IC.N) = 10*(IC+GWIDE)+1 
GCON(ID.N) = 10*( ID+GWIDE)+1 
GCON( IE ,N) = 10*(IE+GWIDE)+1 
GCON(IF.N) = 10*( IF+GWIDE)+1 
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C BACK EDGE 

ELSEIF (GIH(I). EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEF( I ,N)=GXX31 
GOEF(IC,N)=GXX32 
GOEF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEFC ID, N)=GXX 
GOEF(IF,N)=GXX 
N=N+1 

GCON( I ,N) = 10*(I-GWIDE)+1 

GCON( IB ,N) = 10*(IB-GWIDE)+1 
GCON( IC,N) = 10*(IC-GWIDE)+1 
GCON(ID,N) = 10*(ID-GWIDE)+1 
GCON(IE.N) = 10*( IE-GWIDE)+1 
GCON(IF,N) = 10*( IF-GWIDE)+1 

C BACK COEFFICIENT 

GOEF( I ,N)=GFXB31 
GOEF( IC ,N)=GFXB32 
GOEF( IE ,N)=GFXB33 
GOEFC IB, N)=GFXB 
GOEF(ID,N)=GFXB 
GOEFC IF,N)=GFXB 
N=N+1 

GCON( I ,N) = 7511 
GCON(IB.N) = 7511 
GCON(IC.N) = 7511 
GCON(ID.N) = 7511 
GCON( IE ,N) = 7511 
GCON(IF,N) = 7511 

ENDIF 

TOP COEFFICIENT 
GOEF(I,N)=GZE31 
GOEF( IC,N)=GZE32 
GOEF( IE ,N)=GZE33 
GOEF( IB ,N)=GZC31 
GOEF( ID,N)=GZC32 
GOEF(IF,N)=GZC33 
N=N+1 

GCON(I,N) = 7511 
GCON( IC,N) = 7511 
GCON( IE ,N) = 7511 
GCON( IB ,N) = 10*1+1 
GC0N( ID,N) = 10*IC+1 
GCON( IF ,N) = 10*IE+1 

BOTTOM COEFFICIENT 

G0EF( I, N) = GZE31 
GOEFC IC,N) = GZE32 
G0EF( IE, N) = GZE33 
G0EF( IB, N) = GZC31 
GOEFC ID, N) = GZC32 
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GOEF( IF,N) = GZB 
N=N+1 

GC0N(I,N) =10*(I+GPL)+1 
GCON( IB ,N) = 10*(IB+GPL)+1 
GCON( IC ,N) = 10*(IC+GPL)+1 
GCON(ID,N) = 10*( ID+GPL)+1 



HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) , GJH( I ) ) 
N=N+1 

GC0N( I ,N)=9991 
ENDIF 



ELSEIF( (GIH( I) . EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH(I).NE. 1. OR. GJH(I) 
+.NE.GWIDE)) THEN 

DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH( I ) ,GJH( I) ). EQ. 0. 0) THEN 
GC0N( I ,N) = 6 

ELSE 

GC0N(I,N) = 7 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GCON(IB.N) = 6 
GCON(ID.N) = 6 
GCON( IC ,N) = 6 
GC0N( IE ,N) = 6 
GCON( IF ,N) = 6 



LEFT AND RIGHT COEFFICIENTS 
LEFT COEFFICIENT 
G0EF( I ,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF( IB ,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 

GCON( I ,N) = 10*(I-1)+1 
GC0N( IB ,N) = 10*( IB-1)+1 
GCON( IC ,N) = 10*(IC-1)+1 
GCON( ID,N) = 10*(ID-1)+1 
GCON( IE ,N) = 10*(IE-1)+1 
GC0N( IF,N) = 10*(IF-1)+1 
C 

C RIGHT COEFFICIENT 
GOEF( I ,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF( IB ,N) = GYY 
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GOEF(ID,N) = GYY 
GOEF( IF,N) = GYY 
N=N+1 

GCON(I,N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* (IB+1) + 1 
GC0N(IC,N) = 10 * (IC+1) + 1 

GC0N(ID,N) = 10 * (ID+1) + 1 

GC0N( IE ,N) = 10 * ( IE+1) + 1 

GC0N( IF ,N) = 10 * (IF+1) + 1 

FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 

FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
C FRONT COEFFICIENT 
G0EF( I ,N)=GFXB31 
G0EF( IC ,N)=GFXB32 
G0EF( IE ,N)=GFXB33 
G0EF( IB ,N)=GFXB 
G0EF( ID,N)=GFXB 
G0EF( IF ,N)=GFXB 
N=N+1 

GC0N( I ,N) =7521 
GC0N( IB ,N) = 7521 
GC0N( IC ,N) = 7521 
GCON( ID,N) = 7521 
GCON( IE ,N) = 7521 
GC0N( IF ,N) = 7521 

C BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
G0EF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF,N)=GXX 
N=N+1 

GCON(I,N) =10*(I+GWIDE)+1 
GCON( IB ,N) = 10*(IB+GWIDE)+1 
GCON(IC,N) = 10*( IC+GWIDE)+1 
GCON(ID,N) = 10*( ID+GWIDE)+1 
GCON( IE ,N) = 10*(IE+GWIDE)+1 
GCON(IF,N) = 10*( IF+GWIDE)+1 

C BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEF( I ,N)=GXX31 
GOEF( IC ,N)=GXX32 
G0EF( IE ,N)=GXX33 
GOEF( IB ,N)=GXX 
GOEF( ID,N)=GXX 
G0EF( IF ,N)=GXX 
N=N+1 

GCON(I,N) =10*(I-GWIDE)+1 
GC0N( IB ,N) = 10*(IB-GWIDE)+1 
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GCON(IC.N) = 10*(IC-GWIDE)+1 
GCON( ID ,N) = 10*(ID-GWIDE)+1 
GCON( IE ,N) = 10*(IE-GWIDE)+1 
GCON( IF,N) = 10*(IF-GWIDE)+1 

C BACK COEFFICIENT 
G0EF( I ,N)=GFXB31 
G0EF( IC , N)=GFXB32 
G0EF( IE , N)=GFXB33 
G0EF(IB,N) =GFXB 
G0EF( ID,N)=GFXB 
GOEF ( I F , N ) =GFXB 
N=N+1 

GCON(I,N) = 7511 
GCON( IB,N) = 7511 
GCON( IC ,N) = 7511 
GCON(ID,N) = 7511 
GCON( IE ,N) = 7511 
GCON( IF,N) = 7511 
ENDIF 



TOP COEFFICIENT 
G0EF( I ,N)=GZE31 
G0EF( IC , N)=GZE32 
G0EF( IE,N)=GZE33 
G0EF( IB ,N)=GZC31 
G0EF( ID,N)=GZC32 
GOEF( IF ,N)=GZC33 



N=N+1 

GCON( I ,N) = 7511 
GCON( IC,N) = 7511 
GCON(IB,N) = 10*1+1 
GC0N( ID,N) = 10*IC+1 
GC0N( IE ,N) = 7511 
GC0N( IF ,N) = 10*IE+1 



BOTTOM COEFFICIENT 

G0EF(I,N) = GZE31 
GOEF(IC,N) = GZE32 
GOEFC IE ,N) = GZE33 
GOEF(IB,N) = GZC31 
GOEF( ID, N) = GZC32 
GOEFC IF, N) = GZB 

N=N+1 

GCONC I ,N) =10*(I+GPL)+1 
GCONC IB ,N) = 10*C IB+GPL)+1 
GCONC IC,N) = 10*( IC+GPL)+1 
GCONC ID, N) = 10*C ID+GPL)+1 
GCONC IE, N) = 10*( IE+GPL)+1 
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C 

C HEAT INPUT 

IF(GEAT(GIH( I ) ,GJH( I ) ) . NE. 0. 0) THEN 
G0EF( I ,N) = GEAT( GIH( I ) , GJH( I ) ) 

N=N+1 

GCON( I ,N)=9991 
ENDIF 

CXCXCXC 

CXCXCXC 

C 

C 

C LEFT AND RIGHT EDGES EXCLUDING CORNERS 

C 

ELSEIF( (GIH( I). EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH(I).EQ. l.OR. GJH(I) 
+. EQ. GWIDE)) THEN 
C 

C DETERMINE COEFFICIENTS FOR TOP LAYER 

IF (GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
GCON( I , N) = 7 

ELSE 

GC0N( I ,N) = 6 
ENDIF 
C 

C CONNECTIONS FOR EPOXY LAYER 

GCON( IB ,N) = 6 
GCON( IC ,N) = 6 
GCON( ID ,N) = 6 
GCON( IE,N) = 6 
GCON( IF ,N) = 6 
C 
C 

C LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 

C 

C LEFT EDGE 

IF (GJH(I).EQ. 1) THEN 
C LEFT COEFFICIENT 

GOEF( I ,N) = GYLR31 
GOEF(IC.N) = GYLR32 
GOEF(IE.N) = GYLR33 
GOEF( IB ,N) = GYLR 
GOEF(ID.N) = GYLR 
GOEF(IF.N) = GYLR 
N=N+1 

GCON(I,N) = 7551 
GCON(IB.N) = 7551 
GCON(IC,N) = 7551 
GCON( ID ,N) = 7551 
GCON( IE ,N) = 7551 
GCON(IF.N) = 7551 
C 

C RIGHT COEFFICIENT 

GOEF( I ,N) = GYY31 
GOEF(IC.N) = GYY32 
GOEF( IE ,N) = GYY33 
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c 

c 



GOEF( IB ,N) = GYY 
G0EF( ID,N) = GYY 
GOEF( IF,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * 
GC0N( IB ,N) 

GC0N( IC ,N) 

GC0N( ID,N) 

GC0N( IE,N) 

GC0N( IF ,N) 



d+1) + 1 

10* (IB+1) + 1 
10 * (IC+1) + 1 
10 * (ID+1) + 1 
10 * (IE+1) + 1 
10 * ( IF+1) + 1 



RIGHT EDGE 

ELSEIF (GJH(I).EQ. GWIDE) THEN 

LEFT COEFFICIENT 
G0EF( I ,N) = GYY31 
G0EF( IC ,N) = GYY32 
G0EF( IE ,N) = GYY33 
G0EF( IB,N) = GYY 
G0EF( ID ,N) = GYY 
G0EF(IF,N) = GYY 
N=N+1 

GC0N( I ,N) = 10*(I-1)+1 



GCON(IB.N) 
GC0N( IC ,N) 
GC0N( ID ,N) 
GC0N( IE ,N) 
GC0N( IF ,N) 



10*(IB-1)+1 
10*(IC-1)+1 
10*(ID-1)+1 
10*(IE-1)+1 
10*( IF- 1)+1 



RIGHT COEFFICIENT 



G0EF(I,N) 
G0EF( IC ,N) 
GOEF( IE ,N) 
G0EF(IB,N) 
G0EF( ID ,N) 
G0EF( IF,N) 
N=N+1 
GC0N(I,N) 
GC0N( IB ,N) 
GC0N( IC ,N) 
GC0N( ID ,N) 
GC0N( IE ,N) 
GC0N( IF ,N) 
ENDIF 



GYLR31 

GYLR32 

GYLR33 

GYLR 

GYLR 

GYLR 

7541 

7541 

7541 

7541 

7541 

7541 



FRONT COEFFICIENT 
GOEF(I,N) =GXX3 1 
GOEF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
G0EF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
G0EF( IF,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I-GWIDE)+1 
GCON( IB,N) = 10*(IB-GWIDE)+1 
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GCON'( IC ,N) = 10*(IC-GWIDE)+1 
GCON(ID,N) = 10*(ID-GWIDE)+1 
GCON( IE ,N) = 10*(IE-GWIDE)+1 
GCON( IF ,N) = 10*(IF-GWIDE)+1 
BACK COEFFICIENT 
GOEF( I ,N)=GXX31 
G0EF( IC ,N)=GXX32 
GOEF( IE ,N)=GXX33 
G0EF( IB ,N)=GXX 
G0EF( ID,N)=GXX 
G0EF( IF ,N)=GXX 
N=N+1 

GCON(I.N) =10*( I+GWIDE)+1 



GCON( IB ,N) 
GCON( IC ,N) 
GCON( ID,N) 
GCON( IE ,N) 
GCON( IF ,N) 



10*( IB+GWIDE)+1 
10*( IC+GWIDE)+1 
10*( ID+GWIDE)+1 
10*( IE+GWIDE)+1 
10*( IF+GWIDE)+1 



TOP COEFFICIENT 
GOEF( I ,N)=GZE31 
GOEF( IC ,N)=GZE32 
G0EF(ID,N)=GZE33 
GOEF( IB ,N)=GZC31 
G0EF(ID,N)=GZC32 
GOEF( IF ,N)=GZC33 
N=N+1 

GCON(I,N) = 7511 



GCON( IC ,N) 
GCON( IB ,N) 
GCON( ID ,N) 
GCON( IE ,N) 
GCON( IF ,N) 



7511 

10 * 1+1 

10*IC+1 

7511 

10*IE+1 



BOTTOM COEFFICIENT 



G0EF(I,N) 
G0EF( IC ,N) 
G0EF( IE ,N) 



GZE31 

GZE32 

GZE33 



G0EF( IB ,N) 
G0EF( ID ,N) 
G0EF( IF ,N) 
N=N+1 
GCON(I,N) 
GCON( IB ,N) 
GC0N( IC ,N) 
GC0N(ID,N) 
GC0N( IE ,N) 
GC0N( IF ,N) 



GZC31 

GZC32 

EZB 

10*(I+GPL)+1 
10*(IB+GPL)+1 
10*( IC+GPL)+1 
10*(ID+GPL)+1 
10*(IE+GPL)+1 
10*( IF+GPL)+1 



HEAT INPUT 

IF(GEAT(GIH(I),GJH(I)).NE. 0.0) THEN 
G0EF(I,N) = GEAT( GIH( I) ,GJH( I) ) 
N=N+1 
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GCON( I ,N)=9991 
ENDIF 



****DETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN edge******* 



ELSEIF( ( GIH( I ) . NE. 1. OR. GIH( I ) . NE. GDEEP) . AND. ( GJH( I ) . NE. 1. OR. GJH( I ) 
+. NE. GWIDE)) THEN 

DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I),GJH(I)).NE. 0.0) THEN 
GCON(I,N) = 7 

ELSE 

GCON(I,N) = 6 
ENDIF 



CONNECTIONS FOR EPOXY LAYER 
GCON( IB,N) = 6 
GCON( IC,N) = 6 
GCON(ID.N) = 6 
GC0N( IE ,N) = 6 
GCON(IF,N) = 6 



LEFT COEFFICIENT 
GOEF( I ,N) = GYY31 
GOEF(IC,N) = GYY32 
GOEF(IE,N) = GYY33 
GOEF( IB ,N) = GYY 
GOEF(ID,N) = GYY 
GOEF(IF,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * (1-1) + 1 
GC0N( IB ,N) = 10* (IB-1) + 1 
GCON(IC,N) = 10 * (IC-1) + 1 

GCON(ID,N) = 10 * (ID-1) + 1 

GCON( IE ,N) = 10 * (IE-1) + 1 

GCON( IF ,N) = 10 * ( IF- 1) + 1 

C RIGHT COEFFICIENT 
G0EF( I ,N) = GYY31 
GOEF(IC.N) = GYY32 
G0EF( IE ,N) = GYY33 
GOEF(IB.N) = GYY 
G0EF(ID,N) = GYY 
GOEF( IF,N) = GYY 
N=N+1 

GC0N( I ,N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 

GCON(ID,N) = 10 * ( ID+1) + 1 

GCON( IE ,N) = 10 * ( IE+1) + 1 
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GCONC IF,N) = 10 * (IF+1) + 1 

FRONT COEFFICIENT 
G0EF(I,N)=GXX31 
G0EF( IC ,N)=GXX32 
G0EF( IE , N)=GXX33 
G0EF( IB ,N)=GXX 
G0EF( ID,N)=GXX 
GOEFC IF, N)=GXX 
N=N+1 

GC0N( I , N) =10*(I-GWIDE)+1 
GC0N( IB ,N) = 10*(IB-GWIDE)+1 
GC0N( IC,N) = 10*(IC-GWIDE)+1 
GC0N( ID, N) = 10*(ID-GWIDE)+1 
GCONC IE, N) = 10*C IE-GWIDE)+1 
GCONC IF, N) = 10*C IF-GWIDE)+1 

BACK COEFFICIENT 
G0EF( I ,N)=GXX31 
GOEFCIC,N)=GXX32 
GOEFC IE ,N)=GXX33 
GOEFCIB,N)=GXX 
G0EF( ID, N)=GXX 
GOEFC IF, N)=GXX 
N=N+1 

GCONC I, N) =10*CI+GWIDE)+1 
GCONC IB, N) = 10*CIB+GWIDE)+1 
GCONC IC,N) = 10*( IC+GWIDE)+1 
GCONC ID, N) = 10*CID+GWIDE)+1 
GCONC IE, N) = 10*C IE+GWIDE)+1 
GCONC IF, N) = 10*CIF+GWIDE)+1 

TOP COEFFICIENT 
GOEFC I ,N)=GZE31 
GOEF( IC,N)=GZE32 
GOEFC IE, N)=GZE33 
GOEFC IB, N)=GZC31 
GOEFC ID, N)=GZC32 
GOEFC IF, N)=GZC33 
N=N+1 

GCONC I ,N) = 7511 
GCONC IB, N) = 10*1+1 
GCONC IC,N) = 7511 
GCONC ID, N) = 10*IC+1 
GCONC IE,N) = 7511 
GCONC IF, N) = 10*IE+1 



BOTTOM COEFFICIENT 

GOEFC I ,N) = GZE31 
GOEFC IC,N) = GZE32 
GOEFC IE, N) = GZE33 
GOEFC IB, N) = GZC31 
GOEFC ID, N) = GZC32 
GOEFC IF, N) = GZC33 
N=N+1 

GCONC I, N) =10*CI+GPL)+1 
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GCON( IB ,N) = 10*(IB+GPL)+1 
GCON( IC,N) = 10*(IC+GPL)+1 
GCON( ID,N) = 10*(ID+GPL)+1 
GCON(IE.N) = 10*(IE+GPL)+1 
GCON(IF,N) = 10*( IF+GPL)+1 
C 

C HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
GOEF( I ,N) = GEAT( GIH( I ) , GJH( I ) ) 

N=N+1 

GCON( I ,N)=9991 

ENDIF 

ENDIF 

90 CONTINUE 

C 

C GENERATE DATA FILE 

C 

OPEN ( 3 , F I LE=NAME , FORM= ' FORMATTED ' , ACCESS=' DIRECT' ,RECL=108 , STATUS 
+= ' NEW ' ) 

WRITE (3, 9 09) DATAF 
909 FORMAT( 1X,A79) 

WRITE ( 3 , 9 08 ) COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908 F0RMAT(2X,9(I3,5X)) 

WRITE(3,907) ZER, ZER, ZER 
907 FORMAT( 2X,3( 13 ,5X) ) 

WRITE( 3 , 9081) NMAX , TMAX , HTRS ,D1 ,D2 ,D3 ,D4 ,D5 ,D6 ,D7 
9081 FORMAT( 2X, 9( 13 , 5X) ) 

WRITE ( 3 , 9 05 ) ACC , DAMP , MAXIT , CONFAC ,GIBT 

905 FORMAT( 1X,2(F9. 7 , IX) , 14X, 12 , IX, F9. 7,1X,F9.5) 

WRITE( 3 , 9 06 ) GPRT , GFT , GBT , GRT , GLT , GWRT 

906 F0RMAT(1X,6(F12. 3, IX)) 

DO 112 1=1 ,LOCVAR 

WRITE ( 3 , 9 10 0 ) GCON( 1,1), GCON( 1,2), GCON( 1,3), GCON( 1,4), GCON( 1 , 5 ) , GC 
+ON( I , 6 ) , GCON( 1,7), GCON( 1,8) 

9100 FORMAT( 14 , 3X, 7( 15 , 7X) ) 

WRITE( 3 , 9 1 10 ) GOEF( 1,1), GOEF( 1,2), GOEFC 1,3), G0EF( 1,4), GOEF( I , 5 ) , GO 
+EF( 1,6) ,GOEF( 1,7) 

9110 FORMAT( 7(F9. 3,3X)) 

112 CONTINUE 
CLOSE (3) 

CALL CLS 

WRITE(*,999) NAME 

999 FORMAT( //// , ' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 

+NAMED ' ,A6 , /////, ' <PRESS ENTER TO CONTINUE>') 

READ(* ,5912) ANS 
5912 FORMAT(Al) 

END 



C 

C 

C 

C 

C 

C 

c 

c 



TITLE: MODEL BUILDER 

SUBROUTINE: S4 

DATE: 09 JUL 91 

AUTHOR: LT STEVE GLASER 

COMPILER: MICROSOFT VERSION 4.01 

LINKER: MICROSOFT VERSION 3. 55 

SUBROUTINE S4 IS CALLED FROM SUBROUTINE PCBS3. S4 GENERATES THE 
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OUTPUT DATA FILE FOR THE FOUR COPPER LAYER CASE. 



SUBROUTINE S4 ( GYLR , GFXB , GYY , GXX , GZC4 1 , GZC42 , GZC43 , GZC44 , GZB , GYLR4 
+1 , GYLR42 , GYLR43 , GYLR44 , GFXB41 , GFXB42 , GFXB 4 3 , GFXB44 , GYY41 , GYY42 , GYY 
+43 , GYY44 , GXX4 1 , GXX42 , GXX43 , GXX44 , GZE41 , GZE42 , GZE43 , GZE44 , GPL , GIH , G 
+JH,GEAT,GWIDE,GDEEP,GOEF,GIBT,GPRT,GWRT,GRT,GLT,GFT,GBT,GSEL) 

REAL GYLR , GFXB , GYY , GXX , GZC4 1 , GZB , GYLR4 1 , GFXB 41 , GYY4 1 , GXX4 1 , GZE4 1 , A 
+CC , DAMP , CONFAC , G I BT , GPRT , GWRT , GRT , GLT , GFT , GBT , GYLR42 , GFXB 42 , GYY42 , 
+GXX42 , GZE42 , GZC42 , GYLR43 , GFXB43 , GYY43 , GXX43 , GZE43 , GZC43 , GYLR44 , GFX 
+B44 , GYY44 , GXX44 , GZE44 , GZC44 

REAL GEAT( 100,100) ,GOEF( 1000, 12) 



CHARACTER* 1 SELECT, ANS 
CHARACTER DATAF* 7 9 , NAME*6 

INTEGER GPL , I , N , I B , GWIDE , GDEEP , COUNT , CONTEMP , ZER , GSEL , NMAX , TMAX , HT 
+RS ,D1 ,D2 ,D3 ,D4,D5 ,D6 ,D7 ,MAXIT,LOCVAR, IC , ID , IE , IF, IG, IH 

INTEGER GIH( 1000), GJH( 1000), GCON( 1000, 12) 

C0NTEMP=6 

ZER=0 

NMAX=750 

TMAX=50 

HTRS=6 

Dl=2 

D2=4 

D3=6 

D4=0 

D5=0 

D6=0 

D7=0 

ACC=0. 05 
DAMP=0. 666667 
MAXIT=12 
C0NFAC=0. 8 
C0UNT=8*GPL 
LOCVAR=COUNT 



898 CALL CLS 

WRITE(*,803) 

803 FORMAT(///,' THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE ’ , / , 

+' EXISTING THERMAL ANALYZER, FURTHERMORE, THIS PROGRAM DOES 1 ,/, 

+' NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE',/, 

+' THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ',/, 

+' PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS, AND',/, 

+' SHOULD NOT HAVE ANY SPACES. ',///> 

+' PLEASE ENTER THE DESIRED DATA FILE NAME: ' ,2X ) 
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READ(*,804) NAME 




804 


FORMATC A6) 




552 


WRITE(*,910)NAME 




910 


FORMATC/////,' YOU SELECTED ',A6,' FOR YOUR DATA 

+’,//) 


FILE NAME 


911 


WRITE(*.811) 

FORMATC IS THIS THE DESIRED SELECTION? ENTER Y 




811 


FOR YES AND N 




+FOR NO: ' ,2X, ) 






READ(*,812) ANS 




812 


FORMATC Al) 





c 

IF( ANS. EQ. 'N' ) THEN 
GOTO 898 



ELSE 

CONTINUE 

ENDIF 

IF(ANS. EQ. 'Y') THEN 
GOTO 897 
ELSE 

CALL CLS 
GOTO 552 
897 ENDIF 

C ALLOW THE USER TO PROVIDE TITLE LINE FOR DATA FILE 

833 CALL CLS 
WRITE (*,805) 

805 FORMATC ////,’ ENTER THE DESIRED TITLE TO BE PLACED ON LINE 

+ ' ,/,' NUMBER ONE OF THE OUTPUT DATA FILE: ' ,///, 

+ 1 1 2X , ) 

READ(*,806)’dATAF 

806 FORMAT( A79) 

835 WRITE(*,831) 

831 FORMAT(////,' DO YOU WISH TO CHANGE THE TITLE OF YOUR OUTPU 

+T DATA FILE?',/,' ENTER Y FOR YES AND N FOR NO: ',2X, ) 

READ(*, 832) ANS 

832 FORMATC Al) 

C 

IFCANS. EQ. 'Y') THEN 
GOTO 833 
ELSE 

CONTINUE 

ENDIF 

IFCANS. EQ. 'N' ) THEN 
GOTO 834 
ELSE 

GOTO 835 

834 ENDIF 
C 



c *********************************************************************** 
C**************C0EFFICIENTS FOR EPOXY AND COPPER LAYERS***************** 



c 



DO 90 1=1, GPL 
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N=1 

IB=GPL+I 

IC=2*GPL+I 

ID=3*GPL+I 

IE=4*GPL+I 

IF=5*GPL+I 

IG=6*GPL+I 

IH=7*GPL+I 

C 

Q*********^***********Vr***********QORNERS************^*******^^ f *^**^^**^* 

c 

IF ((GIH(I).EQ. 1. OR. GIH( I). EQ. GDEEP). AND. (GJH( I). EQ. 1. OR. GJH( I). EQ 
+. GWIDE)) THEN 



DETERMINE COEFFICIENTS FOR TOP COPPER 
IF (GEAT(GIH( I) ,GJH( I ) ). EQ. 0.0) THEN 
GCON( I ,N) = 6 

ELSE 

GCON( I ,N) = 7 
ENDIF 



CONNECTIONS 
GCON(IB.N) 
GCON( ID ,N) 
GCON(IC,N) 
GCON( ID,N) 
GCON(IE.N) 
GCON( IF ,N) 
GCON(IG.N) 
GCON(IH,N) 



FOR EPOXY LAYER 
= 6 
= 6 
= 6 
= 6 
= 6 
= 6 
= 6 
= 6 



LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
LEFT EDGE 

IF (GJH(I). EQ. 1) THEN 
C LEFT COEFFICIENT 

G0EF( I ,N) = GYLR41 

GOEF(IC,N) = GYLR42 
G0EF( IE,N) = GYLR43 
GOEF(IG,N) = GYLR44 
GOEF( IB ,N) = GYLR 
G0EF( ID,N) = GYLR 
G0EF( IF,N) = GYLR 
G0EF( IG,N) = GYLR 
G0EF( IH,N) = GYLR 
N=N+1 

GCON( I ,N) = 7551 
GCON( IB ,N) = 7551 
GCON( IC ,N) = 7551 
GCON( ID ,N) = 7551 
GCON( IE ,N) = 7551 
GCON( IF ,N) = 7551 
GCON( IG,N) = 7551 
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GCON(IH.N) = 7551 



RIGHT COEFFICIENT 
GOEF( I ,N) = GYY41 
GOEF(IC,N) = GYY42 
GOEF(IE.N) = GYY43 
GOEF(IG,N) = GYY44 
GOEF( IB,N) = GYY 
GOEF(ID.N) = GYY 
GOEF( IF,N) = GYY 
GOEF( IH,N) = GYY 
N=N+1 

GCON( I ,N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* (IB+1) + 1 
GC0N(IC,N) = 10 * (IC+1) + 1 

GCON(ID.N) = 10 * ( ID+1) + 1 

GC0N( IE ,N) = 10 * (IE+1) + 1 

GC0N(IF,N) = 10 * (IF+1) + 1 

GC0N(IG,N) = 10 * ( IG+1) + 1 

GC0N( IH,N) = 10 * (IH+1) + 1 



RIGHT EDGE 

ELSEIF (GJH(I).EQ. GWIDE) THEN 
LEFT COEFFICIENT 



GOEF(I.N) = GYY41 
G0EF(IC,N) = GYY42 
G0EF( IE ,N) = GYY43 
GOEF(IG.N) = GYY44 
G0EF( IB ,N) = GYY 
G0EF(ID,N) = GYY 
G0EF(IF,N) = GYY 
GOEF(IH.N) = GYY 
N=N+1 



GCON(I,N) 
GC0N( IB ,N) 
GCON(IC.N) 
GCON(ID,N) 
GCON( IE ,N) 
GCON(IF,N) 
GCON( IG,N) 
GC0N(IH,N) 



10 *( I - 1)+1 

10*(IB-1)+1 

10*(IC-1)+1 

10*(ID-1)+1 

10*(IE-1)+1 

10*(IF-1)+1 

10*(IG-1)+1 

10*(IH-1)+1 



RIGHT COEFFICIENT 



GOEF(I,N) = GYLR41 
GOEF(IC,N) = GYLR42 
GOEF( IE ,N) = GYLR43 
G0EF( IG,N) = GYLR44 
GOEF( IB,N) = GYLR 
GOEF( ID ,N) = GYLR 
GOEF( IF,N) = GYLR 
GOEF(IH,N) = GYLR 
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N=N+1 

GCON(I,N) = 7541 
GCON( IB,N) = 7541 
GCON(IC,N) = 7541 
GCON( ID ,N) = 7541 
GCON( IE ,N) = 7541 
GCON( IF,N) = 7541 
GCON( IG,N) = 7541 
GCON(IH,N) = 7541 
ENDIF 

FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
C FRONT COEFFICIENT 
GOEF ( I , N ) =GFXB4 1 
GOEF( IC ,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
GOEF( IG,N)=GFXB44 
GOEF( IB ,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF ( I F , N ) =GFXB 
GOEF( IH ,N)=GFXB 
N=N+1 

GCON(I,N) = 7521 
GCON( IB , N) = 7521 
GCON( IC ,N) = 7521 
GCON( ID ,N) = 7521 
GCON( IE ,N) = 7521 
GCON( IF ,N) = 7521 
GCON(IG,N) = 7521 
GCON( IH,N) = 7521 
C BACK COEFFICIENT 
GOEF ( I , N ) =GXX4 1 
GOEF( IC ,N)=GXX42 
GOEF ( I E , N ) =GXX43 
GOEF( IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 
N=N+1 

GCON(I,N) =10*(I+GWIDE)+1 
GCON( IB ,N) = 10*( IB+GWIDE)+1 
GCON( IC ,N) = 10*(IC+GWIDE)+1 
GCON( ID ,N) = 10*( ID+GWIDE)+1 
GCON( IE,N) = 10*( IE+GWIDE)+1 
GCON( IF ,N) = 10*(IF+GWIDE)+1 



C BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEF( I ,N)=GXX41 
GOEF( IC ,N)=GXX42 
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GOEF( IE ,N)=GXX43 
G0EF(IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 
N=N+1 

GCON(I,N) = 10*(I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 
GCON(IC,N) = 10*( IC-GWIDE)+1 
GCON( ID,N) = 10*(ID-GWIDE)+1 
GCON( IE ,N) = 10*( IE-GWIDE)+1 
GCON( IF,N) = 10*(IF-GWIDE)+1 
GCON(IG,N) = 10*( IG-GWIDE)+1 
GCON( IH,N) = 10*(IH-GWIDE)+1 

C BACK COEFFICIENT 
GOEF( I ,N)=GFXB41 
G0EF(IC,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
G0EF(IG,N)=GFXB44 
GOEF( IB ,N)=GFXB 
GOEF(ID,N)=GFXB 
GOEF ( I F , N ) =GFXB 
GOEF( IH,N)=GFXB 
N=N+1 

GCON( I ,N) = 7511 
GCON( IB ,N) = 7511 
GCON(IC.N) = 7511 
GCON( ID,N) = 7511 
GCON( IE ,N) = 7511 
GCON( IF ,N) = 7511 
GCON(IG,N) = 7511 
GCON( IH,N) = 7511 
ENDIF 

C 

C TOP COEFFICIENT 
GOEF( I ,N)=GZE41 
GOEF(IC,N)=GZE42 
GOEF( IE ,N)=GZE43 
GOEF( IG,N)=GZE44 
GOEF( IB ,N)=GZC41 
GOEF( ID,N)=GZC42 
GOEF(IF,N)=GZC43 
G0EF(IH,N)=GZC44 
N=N+1 

GCON( I ,N) = 7511 
GCON(IC.N) = 7511 
GCON( IE,N) = 7511 
GCON( IB ,N) = 10*1+1 
GCON(ID.N) = 10*IC+1 
GCON(IF.N) = 10*IE+1 
GCON(IG.N) = 7511 
GCON(IH.N) = 10*IF+1 

C 

C BOTTOM COEFFICIENT 

C 
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GOEF(I,N) = GZE41 
GOEF(IC.N) = GZE42 
GOEF(IE.N) = GZE43 
GOEF(IG,N) = GZE44 
GOEF( IB ,N) = GZC41 
GOEF( ID,N) = GZC42 
GOEF(IF.N) = GZC43 
GOEF( IH,N) = GZB 
N=N+1 

GCON( I ,N) =10*( I+GPL)+1 
GCON( IB ,N) = 10*(IB+GPL)+1 
GCON(IC.N) = 10*(IC+GPL)+1 
GCON( ID ,N) = 10*(ID+GPL)+1 
GCON(IE.N) = 10*(IE+GPL)+1 
GCON(IF.N) = 10*(IF+GPL)+1 
GCON(IG.N) = 10*(IG+GPL)+1 
GCON(IH.N) = 10*(IH+GPL)+1 
C HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
GOEF(I.N) = GEAT( GIH( I ) , G JH( I ) ) 
N=N+1 

GCON( I ,N)=9991 
ENDIF 



C 

ELSEIF((GIH(I). EQ. 1. OR. GIH(I). EQ. GDEEP). AND. (GJH(I).NE. 1. OR. GJH(I) 
+. NE. GWIDE)) THEN 
C 

C DETERMINE NUMBER OF CONNECTIONS FOR COPPER AND EPOXY LAYERS 
IF (GEAT(GIH( I) ,GJH( I) ). EQ. 0. 0) THEN 
GCON(I,N) = 6 

ELSE 

GCON(I,N) = 7 
ENDIF 

CONNECTIONS FOR EPOXY LAYER 
GCON( IB ,N) = 6 
GCON( ID,N) = 6 
GCON(IC,N) = 6 
GCON( IE,N) = 6 
GCON( IF ,N) = 6 
GCON( IG,N) = 6 
GCON( IH,N) = 6 
C 

C LEFT AND RIGHT COEFFICIENTS 

C LEFT COEFFICIENT 

GOEF(I,N) = GYY41 
GOEF( IC , N) = GYY42 
GOEF( IE ,N) = GYY43 
GOEF( IG, N) = GYY44 
GOEF( IB ,N) = GYY 
GOEF( ID,N) = GYY 
GOEF( IF ,N) = GYY 
GOEF( IH,N) = GYY 
N=N+1 

GCON(I,N) = 10*(I-1)+1 
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GCON( IB ,N) = 10*(IB-1)+1 
GCON(IC.N) = 10*(IC-1)+1 
GCON(ID,N) = 10*(ID-1)+1 
GCON( IE,N) = 10*(IE-1)+1 
GCON( IF,N) = 10*(IF-1)+1 
GCON( IG,N) = 10*(IG-1)+1 
GCON(IH,N) = 10*(IH-1)+1 

RIGHT COEFFICIENT 
GOEF( I ,N) = GYY41 

GOEF( IC,N) = GYY42 
GOEF( IE,N) = GYY43 
GOEF( IG,N) = GYY44 
GOEF( IB ,N) = GYY 
GOEF( ID,N) = GYY 
GOEF( IF,N) = GYY 
GOEF(IH.N) = GYY 
N=N+1 

GCON(I,N) = 10 * ( 1+1) + 1 
GCON( IB ,N) = 10* (IB+1) + 1 
GCON( IC,N) = 10 * (IC+1) + 1 
GCON(ID,N) = 10 * (ID+1) + 1 
GCON( IE ,N) = 10 * (IE+1) + 1 
GCON(IF.N) = 10 * (IF+1) + 1 
GCON(IG.N) = 10 * ( IG+1) + 1 
GCON(IH,N) = 10 * (IH+1) + 1 

FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 

FRONT EDGE 

IF (GIH(I).EQ. 1) THEN 
C FRONT COEFFICIENT 
GOEF( I ,N)=GFXB41 
G0EF(IC,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
G0EF(IG,N)=GFXB44 
GOEF ( I B , N ) =GFXB 
GOEF(ID,N)=GFXB 
GOEF(IF,N)=GFXB 
GOEF( IH,N)=GFXB 
N=N+1 

GCON(I,N) =7521 
GCON(IB.N) = 7521 
GCON(IC,N) = 7521 
GCON( ID,N) = 7521 
GCON(IE.N) = 7521 
GCON( IF ,N) = 7521 
GC0N(IG,N) = 7521 
GCON(IH,N) = 7521 

C BACK COEFFICIENT 
GOEF( I ,N)=GXX41 
G0EF(IC,N)=GXX42 
GOEF( IE ,N)=GXX43 
G0EF(IG,N)=GXX44 
GOEF( IB ,N)=GXX 
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GOEF(ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH,N)=GXX 
N=N+1 

GCON( I , N) =10*( I+GWIDE)+1 
GCON( IB ,N) = 10*(IB+GWIDE)+1 
GCON( IC,N) = 10*(IC+GWIDE)+1 
GCON(ID.N) = 10*( ID+GVIDE)+1 
GCON( IE ,N) = 10*(IE+GWIDE)+1 
GCON( IF,N) = 10*(IF+GWIDE)+1 
GCON( IG,N) = 10*(IG+GWIDE)+1 
GCON( IH,N) = 10*(IH+GWIDE)+1 
C BACK EDGE 

ELSEIF (GIH(I).EQ. GDEEP) THEN 
C FRONT COEFFICIENT 

GOEF ( I , N ) =GXX4 1 
GOEF(IC,N)=GXX42 
GOEF(IE,N)=GXX43 
G0EF(IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF ( IH , N ) =GXX 
N=N+1 

GCON(I,N) =10*(I-GWIDE)+1 
GCON(IB.N) = 10*(IB-GWIDE)+1 
GCON(IC,N) = 10*(IC-GWIDE)+1 
GCON( ID ,N) = 10*( ID-GWIDE)+1 
GCON(IE.N) = 10*(IE-GWIDE)+1 
GCON(IF.N) = 10*(IF-GWIDE)+1 
GCON( IG,N) = 10*(IE-GWIDE)+1 
GCON(IH,N) = 10*(IF-GWIDE)+1 
C BACK COEFFICIENT 
GOEF( I ,N)=GFXB41 
GOEF( IC ,N)=GFXB42 
GOEF( IE ,N)=GFXB43 
GOEF( IG,N)=GFXB44 
GOEF( IB ,N)=GFXB 
GOEF( ID ,N)=GFXB 
GOEF( IF,N)=GFXB 
GOEF( IH,N)=GFXB 
N=N+1 

GCON( I ,N) =7511 
GCON( IB,N) = 7511 
GCON(IC.N) = 7511 
GCON( ID , N) =7511 
GCON( IE ,N) = 7511 
GCON( IF,N) = 7511 
GCON(IG.N) = 7511 
GCON( IH,N) = 7511 
ENDIF 

TOP COEFFICIENT 
GOEF( I ,N)=GZE41 
GOEF( IC ,N)=GZE42 
GOEF( IE,N)=GZE43 
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GOEF( IG , N)=GZE44 
G0EF(IB,N)=GZC41 
GOEF( ID ,N)=GZC42 
G0EF(IF,N)=GZC43 
GOEFC IH,N)=GZC44 
N=N+1 

GCON( I ,N) = 7511 
GCON(IC,N) = 7511 
GC0N( IB ,N) = 10*1+1 
GC0N( ID,N) = 10*IC+1 
GC0N( IE ,N) = 7511 
GCONC IF, N) = 10*IE+1 
GCON(IG.N) = 7511 
GC0N(IH,N) = 10*IG+1 

BOTTOM COEFFICIENT 

GOEFC I ,N) = GZE41 
G0EF(IC,N) = GZE42 
G0EF( IE, N) = GZE43 
G0EF( IG,N) = GZE44 
G0EF( IB, N) = GZC41 
G0EF( ID, N) = GZC42 
G0EF( IF, N) = GZC43 
G0EF( IH,N) = GZB 
N=N+1 

GC0N( I ,N) =10*C I+GPL)+1 
GC0N( IB ,N) = 10*( IB+GPL)+1 
GC0N(IC,N) = 10*( IC+GPL)+1 
GCONC ID,N) = 10*C ID+GPL)+1 
GC0N( IE ,N) = 10*( IE+GPL)+1 
GCONC IF, N) = 10*C IF+GPL)+1 
GC0N( IG,N) = 10*C IG+GPL)+1 
GCONC IH,N) = 10*C IH+GPL)+1 

HEAT INPUT 

IFCGEATCGIHCI),GJH(I)).NE. 0. 0) THEN 
GOEFC I, N) = GEAT(GIHCI),GJH(I)) 
N=N+1 

GC0N( I ,N)=999 1 
ENDIF 



LEFT AND RIGHT EDGES EXCLUDING CORNERS 

ELSEIF( C GIH( I). EQ. 1. OR. GIH( I) . EQ. GDEEP) . AND. CGJH(I).EQ. 1. OR. GJH(I) 
+. EQ. GWIDE)) THEN 
C 

C DETERMINE COEFFICIENTS FOR TOP LAYER 
IF (GEAT(GIH(I),GJHCI)).NE. 0. 0) THEN 
GCONC I, N) = 7 

ELSE 

GCONC I, N) = 6 
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ENDIF 



CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 
GCON(IC,N) = 6 
GCON( ID,N) = 6 
GCON(IE,N) = 6 
GC0N( IF,N) = 6 
GCON(IG.N) = 6 
GCON( IH,N) = 6 

LEFT AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
LEFT EDGE 

IF (GJH(I).EQ. 1) THEN 
C LEFT COEFFICIENT 

GOEF( I ,N) = GYLR41 
GOEF(IC.N) = GYLR42 
GOEF(IE.N) = GYLR43 
GOEF( IG,N) = GYLR44 
G0EF( IB ,N) = GYLR 
GOEF( ID ,N) = GYLR 
GOEF(IF.N) = GYLR 
GOEF( IH,N) = GYLR 
N=N+1 

GCON( I ,N) = 7551 
GCON(IB.N) = 7551 
GCON(IC.N) = 7551 
GCON(ID.N) = 7551 
GCON( IE ,N) = 7551 
GCON(IF.N) = 7551 
GCON( IG,N) = 7551 
GCON(IH.N) = 7551 
C RIGHT COEFFICIENT 

GOEF(I,N) = GYY41 
GOEF(IC.N) = GYY42 
GOEF(IE.N) = GYY43 
GOEF(IG,N) = GYY44 
GOEF(IB.N) = GYY 
GOEF(ID.N) = GYY 
GOEF(IF.N) = GYY 
GOEF(IH,N) = GYY 
N=N+1 

GCON(I.N) = 10 * (1+1) + 1 
GC0N( IB ,N) = 10* (IB+1) + 1 
GCON(IC,N) = 10 * (IC+1) + 1 

GC0N( ID,N) = 10 * ( ID+1) + 1 

GCON( IE ,N) = 10 * ( IE+1) + 1 

GCON(IF.N) = 10 * ( IF+1) + 1 

GCON(IG,N) = 10 * (IG+1) + 1 

GC0N( IH,N) = 10 * (IH+1) + 1 

C RIGHT EDGE 

ELSEIF (GJH(I).EQ. GWIDE) THEN 
C 

C LEFT COEFFICIENT 

GOEF( I ,N) = GYY41 
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GOEF(IC,N) = GYY42 
GOEF(IE.N) = GYY43 
GOEF( IG,N) = GYY44 
GOEF( IB ,N) = GYY 
GOEFC ID, N) = GYY 
GOEF( IF,N) = GYY 
GOEF( IH,N) = GYY 
N=N+1 

GCON(I,N) = 10*(I-1)+1 
GCON( IB ,N) = 10*(IB-1)+1 
GCON(IC,N) = 10*(IC-1)+1 
GCON( ID,N) = 10*(ID-1)+1 
GCON( IE ,N) = 10*(IE-1)+1 
GCON( IF,N) = 10*(IF-1)+1 
GCON( IG,N) = 10*(IG-1)+1 
GCON(IH,N) = 10*(IH-1)+1 
RIGHT COEFFICIENT 

GOEF( I ,N) = GYLR41 

GOEFC IC,N) = GYLR42 
GOEF( IE ,N) = GYLR43 
GOEF( IG,N) = GYLR44 
GOEF( IB ,N) = GYLR 
GOEF( ID,N) = GYLR 
GOEF( IF,N) = GYLR 
GOEF( IH,N) = GYLR 
N=N+1 

GCON( I ,N) = 7541 

GCON( IB,N) = 7541 
GCON(IC,N) = 7541 
GCON(ID,N) = 7541 
GCONC IE ,N) = 7541 
GCON( IF,N) = 7541 
GCONC IG,N) = 7541 
GCONC IH,N) = 7541 
ENDIF 

FRONT COEFFICIENT 
GOEFC I, N) =GXX4 1 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEFC ID,N)=GXX 
GOEF(IF,N)=GXX 
GOEFC IH,N)=GXX 
N=N+1 

GCONC I, N) =10*C I-GWIDE)+1 
GCONC IB, N) = 10*(IB-GWIDE)+1 
GCONC IC,N) = 10*(IC-GWIDE)+1 
GCONC ID, N) = 10*(ID-GWIDE)+1 
GCONC IE, N) = 10*(IE-GWIDE)+1 
GCONC IF, N) = 10*(IF-GWIDE)+1 
GCONC IG,N) = 10*CIG-GWIDE)+1 
GCONC IH,N) = 10*C IH-GWIDE)+1 
C BACK COEFFICIENT 
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GOEF(I,N) =GXX4 1 
G0EF(IC,N)=GXX42 
GOEF( IE ,N)=GXX43 
G0EF(IG,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF( ID ,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF(IH,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I+GWIDE)+1 
GCON( IB , N) = 10*(IB+GWIDE)+1 
GCON(IC,N) = 10*(IC+GWIDE)+1 
GCON( ID,N) = 10*( ID+GWIDE )+l 
GCON( IE , N) = 10*( IE+GWIDE )+l 
GCON( IF,N) = 10*(IF+GWIDE)+1 
GCON( IG,N) = 10*( IG+GWIDE)+1 
GCON(IH,N) = 10*(IH+GWIDE)+1 
C TOP COEFFICIENT 
GOEF( I ,N)=GZE41 
GOEF( IC,N)=GZE42 
GOEF( IE , N)=GZE43 
G0EF(IG,N)=GZE44 
G0EF(IB,N)=GZC41 
G0EF(ID,N)=GZC42 
GOEF( IF ,N)=GZC43 
G0EF(IH,N)=GZC44 
N=N+1 

GCON(I,N) = 7511 
GCON(IC,N) = 7511 
GCON( IB ,N) = 10*1+1 
GCON(ID.N) = 10*IC+1 
GCON( IE ,N) = 7511 
GCON( IF ,N) = 10*IE+1 
GCON(IG.N) = 7511 
GCON(IH.N) = 10*IG+1 

BOTTOM COEFFICIENT 



GOEF(I,N) = 
GOEF(IC,N) = 
GOEF(IE.N) = 
GOEF(IG,N) = 
GOEF(IB,N) = 
GOEF( ID,N) = 
GOEF( IF,N) = 
GOEF(IH,N) = 
N=N+1 

GCON(I.N) = 
GCON( IB ,N) = 
GCON(IC,N) = 
GCON( ID,N) = 
GCON( IE ,N) = 
GCON( IF ,N) = 
GCON( IG,N) = 



GZE41 

GZE42 

GZE43 

GZE44 

GZC43 

GZC43 

GZC43 

EZB 

10*( I+GPL)+1 
10*( IB+GPL)+1 
10*( IC+GPL)+1 
10*( ID+GPL)+1 
10*( IE+GPL)+1 
10*( IF+GPL)+1 
10*(IG+GPL)+1 
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GCON( IH,N) = 10*(IH+GPL)+1 
HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ). NE. 0. 0) THEN 
G0EF(I,N) = GEAT( GIH( I ) , GJH( I ) ) 
N=N+1 

GC0N( I ,N)=999 1 
ENDIF 



****DETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE******* 



ELSEIF((GIH( I). NE. 1. OR. GIH(I). NE. GDEEP). AND. (GJH( I). NE. 1. OR. GJH(I) 
+. NE. GWIDE)) THEN 

DETERMINE CONNECTIONS FOR TOP LAYER 
IF (GEAT(GIH(I),GJH(I)).NE. 0. 0) THEN 
GCON( I ,N) = 7 

ELSE 

GC0N(I,N) = 6 
ENDIF 



CONNECTIONS FOR EPOXY LAYER 
GC0N( IB ,N) = 6 
GC0N(IC,N) = 6 
GC0N(ID,N) = 6 
GC0N( IE ,N) = 6 
GC0N( IF ,N) = 6 
GC0N( IG,N) = 6 
GC0N(IH,N) = 6 



C 



LEFT COEFFICIENT 
G0EF(I,N) = GYY41 
G0EF(IC,N) = GYY42 
G0EF( IE,N) = GYY43 
G0EF(IG,N) = GYY44 
G0EF(IB,N) = GYY 
G0EF(ID,N) = GYY 
GOEF(IF.N) = GYY 
G0EF( IH,N) = GYY 
N=N+1 



GC0N( 

GC0N( 

GC0N( 

GC0N( 

GC0N( 

GCON( 

GCON( 

GCON( 

RIGHT 

G0EF( 

G0EF( 

GOEF( 

GOEF( 

GOEF( 

GOEF( 



I »N) 
IB ,N) 
IC,N) 
ID ,N) 
IE,N) 
IF ,N) 

IG, N) 

IH, N) 



10 * ( 1 - 1 ) + 1 
10* (IB-1) + 1 
10 * (IC-1) + 1 
(ID-1) + 1 
(IE-1) + 1 
( IF-1) + 1 
(IG-1) + 1 
( IH-1) + 1 



10 

10 

10 

10 

10 



* 

* 

* 

* 

* 

* 



COEFFICIENT 
I ,N) = GYY41 

IC, N) = GYY42 
IE ,N) = GYY43 
IG,N) = GYY44 
IB ,N) = GYY 

ID , N) = GYY 
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GOEF(IF,N) = GYY 

GOEF( IH,N) = GYY 

N=N+1 

GCON( I ,N) = 10 * ( 1+1) + 1 
GCON( IB ,N) = 10* (IB+1) + 1 
GCON( IC ,N) = 10 * (IC+1) + 1 

GCON( ID ,N) = 10 * ( ID+1) + 1 

GCON( IE ,N) = 10 * ( IE+1) + 1 

GC0N( IF,N) = 10 * ( IF+1) + 1 

GCON(IG,N) = 10 * ( IG+1) + 1 

GC0N(IH,N) = 10 * ( IH+1) + 1 

FRONT COEFFICIENT 
GOEF( I ,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF( I G , N ) =GXX44 
GOEF( IB ,N)=GXX 
GOEF(ID,N)=GXX 
GOEF( IF,N)=GXX 
GOEF(IH,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I-GWIDE)+1 
GCON( IB ,N) = 10*(IB-GWIDE)+1 
GCON( IC ,N) = 10*( IC-GWIDE)+1 
GCON( ID,N) = 10*( ID-GWIDE)+1 
GCON( IE ,N) = 10*(IE-GWIDE)+1 
GCON( IF ,N) = 10*(IF-GWIDE)+1 
GCON(IG.N) = 10*(IG-GVIDE)+1 
GCON( IH,N) = 10*( IH-GWIDE)+1 
C BACK COEFFICIENT 
GOEF( I ,N)=GXX41 
GOEF( IC ,N)=GXX42 
GOEF( IE ,N)=GXX43 
GOEF(IH,N)=GXX44 
GOEF( IB ,N)=GXX 
GOEF( ID,N)=GXX 
GOEF( IF ,N)=GXX 
GOEF( IH ,N)=GXX 
N=N+1 

GCON( I ,N) =10*(I+GWIDE)+1 
GC0N( IB ,N) = 10*( IB+GWIDE)+1 
GCON( IC ,N) = 10*( IC+GWIDE)+1 
GCON( ID ,N) = 10*(ID+GWIDE)+1 
GCON( IE ,N) = 10*( IE+GWIDE)+1 
GCON( IF ,N) = 10*( IF+GWIDE)+1 
GCON( IG,N) = 10*( IG+GWIDE)+1 
GCON( IH,N) = 10*( IH+GWIDE)+1 
C TOP COEFFICIENT 
GOEF( I ,N)=GZE41 
GOEF( IC ,N)=GZE42 
G0EF( IE ,N)=GZE43 
GOEF( IG,N)=GZE44 
GOEF( IB ,N)=GZC41 
GOEF( ID ,N)=GZC42 
GOEF( IF,N)=GZC43 



142 



G0EF(IH,N)=GZC44 

N=N+1 

GCON( I ,N) = 7511 
GCON(IB.N) = 10*1+1 
GCON(IC,N) = 7511 
GCON(ID,N) = 10*IC+1 
GCON(IE.N) = 7511 
GCON( IF,N) = 10*IE+1 
GCON(IG,N) = 7511 
GCON(IH,N) = 10*IG+1 
C 

C BOTTOM COEFFICIENT 

C 

GOEF( I ,N) = GZE41 
GOEF(IC,N) = GZE42 
GOEF(IE,N) = GZE43 
GOEF( IG,N) = GZE44 
GOEF( IB ,N) = GZC41 
GOEF(ID,N) = GZC42 
GOEF( IF,N) = GZC43 
GOEF(IH,N) = EZB 
N=N+1 

GCON( I ,N) =10*(I+GPL)+1 
GCON(IB,N) = 10*(IB+GPL)+1 
GCON(IC.N) = 10*( IC+GPL)+1 
GCON( ID,N) = 10*( ID+GPL)+1 
GCON( IE ,N) = 10*( IE+GPL)+1 
GCON( IF,N) = 10*(IF+GPL)+1 
GCON(IG,N) = 10*(IG+GPL)+1 
GCONCIH.N) = 10*(IH+GPL)+1 
C HEAT INPUT 

IF(GEAT(GIH( I) ,GJH( I) ) . NE. 0. 0) THEN 
GOEFCI ,N) = GEAT( GIH( I ) , G JH( I ) ) 

N=N+1 

GCON(I,N)=9991 

ENDIF 

ENDIF 

90 CONTINUE 

C 

C GENERATE DATA FILE 

C 

OPEN ( 3 , FI LE=NAME , FORM= ' FORMATTED ’ ,ACCESS=' DIRECT’ ,RECL=108, STATUS 
+= 1 NEW 1 ) 

WRITE(3,909) DATAF 
909 F0RMAT(1X,A79) 

WRITEC 3 , 9 08 ) COUNT , CONTEMP , ZER , ZER , ZER , ZER , ZER , ZER , GSEL 
908 FORMAT(2X,9(I3,5X) ) 

WRITE( 3 ,907 ) ZER, ZER, ZER 
907 F0RMAT(2X, 3( 13 ,5X) ) 

WRITEC 3, 9081) NMAX , TMAX , HTRS ,D1 ,D2 ,D3 ,D4,D5 ,D6 ,D7 
9081 FORMATC2X,9CI3,5X)) 

WRITE ( 3 , 905 ) ACC , DAMP , MAXIT , CONFAC , GIBT 

905 FORMATC1X.2CF9. 7 , IX) , 14X, 12 , IX, F9. 7,1X,F9.5) 

WRITEC 3, 906) GPRT , GFT , GBT , GRT , GLT , GWRT 

906 FORMATC IX, 6(F12. 3, IX) ) 

DO 112 1=1 ,LOCVAR 
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VRITE( 3 , 9100) GC0N( 1,1) ,GC0N( 1,2) ,GCON( 1,3) ,GCON( 1,4) ,GCON( I ,5) ,GC 
+0N( 1,6) , GCON( 1,7), GCON( 1,8) 

9100 FORMAT( 14 , 3X , 7 ( 15 , 7X) ) 

WRITE( 3,9110) GOEF( I , 1) ,GOEF( 1,2) ,GOEF( I , 3) , GOEF( I ,4) ,GOEF( I ,5) ,G0 
+EF ( I , 6 ) , GOEF (1,7) 

9110 FORMAT( 7(F9. 3 , 3X) ) 

112 CONTINUE 
CLOSE (3) 

CALL CLS 

WRITE(*,999) NAME 

999 FORMATC////, ' THE OUTPUT DATA HAS BEEN PLACED IN A FILE 

+NAMED ' ,A6 ,/////, ' <PRESS ENTER TO CONTINUE>') 

READ(*,5912) ANS 
5912 FORMATC Al) 

END 
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